Is there a way to make the BLAST results easier to read?
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!
• 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.
Traffic: 2630 users visited in the last hour
Read more here: Source link