Is there a way to make the BLAST results easier to read?

Is there a way to make the BLAST results easier to read?

1

Hi all.

I’m trying to get assembly stat.
As part of that, I would like to evaluate the assembly quality by running BLAST to the fasta file output by Trinity.
I got the result by blastx, but it is difficult to evaluate the result.
Is there any way to evaluate the result of each sequence at once?

I’m referring to table S1 in this paper.
www.ncbi.nlm.nih.gov/pmc/articles/PMC7266049/

I’m very appriciate for your help!


RNA-seq


BLAST


Trinity

• 30 views

To easily generate the count table, you could do the following. Use tabular output in local blast: -outfmt "6 staxid ssciname" then count occurrences by piping into | sort | uniq -c like so:

blastx -db nr -query Trinity.fasta -evalue 1e-6 -outfmt "6 staxid ssciname" `| sort | uniq -c | sort -n > taxon-counts.txt

You may want to adjust parameters, also whether or not staxid and ssciname work depend on how the database was built.


Login
before adding your answer.

Traffic: 2630 users visited in the last hour

Read more here: Source link