I need help with a blast command, please help
Entering edit mode
4 months ago
neto.pvp • 0

I have this command:

blastn -db /home/BLAST_nt/nt -query otu_greedy_0.97_82.fasta -outfmt 5 -out otu_greedy_0.97_82_blastn.xml -evalue 0.001

And i have to include the following conditions:

percent identity of 95%, e-value of 0.001, minimum query coverage of 100%, best hit score edge of 0.05 and best hit overhang of 0.25.

should I include these conditions in the code above? if yes, how can i do that?

linux blast • 386 views
Entering edit mode
4 months ago
Mensur Dlakic ★ 25k

Your command already has all the conditions that can be specified, which is E-value threshold. The remaining options you want can't be set on command line.

I suggest you specify -outfmt 6 and save the output with an appropriate name. What can be included in the tabular output is specified here. Then you have to post-process this file by sorting according to the percent identity column (which is third from left if you use this format without modification) and removing all lines that have <95% in that column. For the other values you will have to be a bit creative and calculate them based on various quantities from the tabular output, but it should be doable.

This could be useful:


Entering edit mode

Thank you very much!! I'm new to this subject and I believe your answer gave me a great idea of where to go next, thanks again!!


Login before adding your answer.

Traffic: 1467 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6