Threshold over percentage identity in blastn
2
0
Entering edit mode
8.1 years ago
Tony ▴ 10

Hello folks, I am currently using blastn on commandline trying to blast a complete genome over a database. when i do this i can get only those hits above 90% percent identity. But i would like to see those hits even if they have 40% percent identity. Could someone help me out with this? Thanks in advance...

command i am currently using is:

"blastn -db db.fasta -query query.fasta -out queryout.tsv -outfmt 6"

blast alignment next-gen gene • 13k views
ADD COMMENT
3
Entering edit mode
5.9 years ago
gb ★ 2.2k

A little late (2.2 years later) but if some one come across this. You need to change the default -task:

blastn -db db.fasta -query query.fasta -out queryout.tsv -perc_identity 40 -outfmt 6 -task blastn
ADD COMMENT
0
Entering edit mode

interesting, any idea why that is? == that you have to add the -task blastn while already asking to do blastn ?

ADD REPLY
1
Entering edit mode

blastn does not do an alignment with all sequences, this will take to long. Before the actual alignments there is a pre "filter" step that makes a selection from the reference database based on exact matches. In other words, it takes substrings from your query and it tries to find sequences that have exactly that string. It works like kmers. This is really fast, if a sequence in the reference already not contains that small overlap it is no use to align it with a heavy algoritm. Megablast searches default with substrings(kmers) of 28 bases and blastn with 11. So if u use megablast all the sequences that does not have an exacted match somewhere of 28 bases in a sequence in the database the sequence will be discarded. Megablast is the default task

ADD REPLY
0
Entering edit mode
8.1 years ago
Tony ▴ 10

Hello everyone., After looking keener into help page, i found it is possible to put thresholds over percentage identity in blastn.

command:

blastn -db db.fasta -query query.fasta -out queryout.tsv -perc_identity 40 -outfmt 6

ADD COMMENT

Login before adding your answer.

Traffic: 1553 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6