This is a test version of Biostars. For the public version, visit https://www.biostars.org.
To get number of hits from blastp output file

Hi all,

I have multiple blastp output (format 6) in a directory, I wish to calculate the number of hits with sequence identity of more than 40% for each output file, therefore, I have tried:

for i in *.tsv; do awk '$3>=40' $i | wc -l; done

However, this command only give me a list of number in the terminal without matching it with the blastp output, any modification that I can do so that I know the number belongs to which blastp output file? Thank you.

genome
• 2,848 views • link
updated by AK AK • written by fec2 • fec2

1 answer

Hi fec2,

Try:

for i in *.tsv; do echo -ne "${i}\t" && awk '$3>=40{print $2}' ${i} | sort -u | wc -l; done

It should returns something like:

blast_out1.tsv  217
blast_out2.tsv  172
blast_out3.tsv  215
• 0 views • link
written by AK AK

Thank you very much! Is it possible to get the result in an output file?

• 0 views • link
written by fec2 • fec2

You're welcome. Try this:

(for i in *.tsv; do echo -ne "${i}\t" && awk '$3>=40{print $2}' ${i} | sort -u | wc -l; done) > output.txt
• 0 views • link
written by AK AK

The command working well, thank you again!

• 0 views • link
written by fec2 • fec2

Hi fec2,

If an answer was helpful you should upvote it, if the answer resolved your question you should mark it as accepted.
Upvote|Bookmark|Accept

• 0 views • link
written by lieven.sterck lieven.sterck

Hi thanks for your comment, I have accepted the answer. Thank you.

• 0 views • link
written by fec2 • fec2

Hi,

May I know where can I find the manual for the meaning of all these command? I am new in this field and but have no clue where to find all these information. Really appreciate your help.

• 0 views • link
written by fec2 • fec2

Hello fec2,

You can use man echo, man awk, man sort, and man wc. I'd recommend "Ch3. Remedial Unix Shell" and "Ch7. Unix Data Tools" in the book: Bioinformatics Data Skills by Vince Buffalo. :-)

• 0 views • link
written by AK AK

Log in to answer this question.

More posts like this