Get query coverage using CLI blast

[hezichia]

I'm running blastall locally but can not find an option to include query coverage in the output. How do I add query coverage to the blastall output?

I would even be happy with query length in the output as query coverage would then be easily calculated using alignment length.

This problem seems to be addressed in this thread: http://seqanswers.com/forums/showthread.php?t=9971 however I was unable to find a query length output option in blast_formatter.

Thanks.

[lzsph]

I haven't used blastall before, I just used NCBI BLAST+ 2.2.26.
For example, with the local blastx tool.

$blastx -db db_name -query input_file -out output -evalue evalue -outfmt "6 qlen" -num_threads 4

Only 6, 7, and 10 option of

-outfmt

allows you add additional column. If you want to add extra column(s), delimited using space ,such as

"6 qseqid sseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore qlen"

The first 12 columns are default columns, you can also delete it. The last qlen is what you want to add.

You can dind the detailed documentation http://computing.bio.cam.ac.uk/local/doc/blast+_user_manual.pdf in section 4.2.26 outfmt.