1
$\begingroup$

I am trying to reproduce locally on my computer what I get running mirbase on their website using BLAST. The search sequences is: mature miRNAs which I had downloaded on my computer and make it as a BLAST database with command:

./makeblastdb -in /home/marianoavino/Downloads/mature.fa -dbtype 'nucl' -out /home/marianoavino/Downloads/mature

then on mirbase I see they use an e-value of 10, which I leave locally. On mirbase at the end of the analysis they give you these parameter setting:

Search parameters

Search algorithm:
    BLASTN 
Sequence database:
    mature 
Evalue cutoff:
    10 
Max alignments:
    100 
Word size:
    4 
Match score:
    +5 
Mismatch penalty:
    -4 

and this is the command line I use on my computer for BLAST

./blastn -db /home/marianoavino/Downloads/mature -evalue 10 -word_size 4 -query /home/marianoavino/Downloads/testinputblast.fasta -task "blastn" -out /home/marianoavino/Downloads/testBLast.out

The results of the two analysis are different, with mirbase finding much more stuff than local BLAST.

Do you have any idea on which parameters I should use on local blast command line to match those listed mirbase parameters in order to get the same answer?

$\endgroup$
0

1 Answer 1

1
$\begingroup$

The main problem I could identify is that your scoring matrix is not the same; as you said miRbase uses match score of +5 and mismatch penalty of -4.

The BLASTN defaults can be found here. However, I ran the local BLAST using the exact command that you mentioned and the default values that it used were:
reward = +1
penalty = -2

Matrix: blastn matrix 1 -2
Gap Penalties: Existence: 0, Extension: 2.5

The other reasons could be the difference in the database itself (the sequences used) or the version of BLAST used in the two cases.

$\endgroup$

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .