This is a test version of Biostars. For the public version, visit https://www.biostars.org.
sorting blast results by sseq coverage

my query blastn search returns - 3 sequences tie ( subject sequence) with identical score , e-value and alignment length.

These 3 sequences are of different length. Is there any way to output results so that the hit with highest subject sequence coverage (or in other words, the shortest hit) is displayed first?

blast
• 1,482 views • link
written by nheb1977 • nheb1977

Hi nheb1977,

Try adding -outfmt "6 qseqid sseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore slen" when running BLAST, then sort the results using sort -k1,1 -k12,12nr -k11,11n -k4,4nr -k13,13n.

• 0 views • link
written by AK AK

0 answers

No answers yet.

Log in to answer this question.

More posts like this