This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Getting scientific name (ssciname) in Blastn output

Hey all,

I would like to get the scientific names of each of my blast hits from blastn run locally. I have downloaded the taxdb.tar.gz taxonomy database from NCBI and have set the $BLASTDB variable to the directory with my blast database and query contigs. In my tabbed output, all I get is an 'N/A' in the spaces where the scientific names should be present. I was wondering if I am missing an important step in the taxonomy database installation?

$ export BLASTDB=$BLASTDB:/path/to/myblastdb
$ blastn -db all_viruses.fasta -query my_contigs.fasta -max_target_seqs 200 -evalue 0.001 -outfmt "6 qseqid sseqid sscinames qcovs pident evalue bitscore length" -out my_blast_results

Thank you in advance, Brie

blast alignment
• 4,055 views • link
written by brianna.flynn • brianna.flynn

2 answers

What blast version are you using? Try ssciname. Not the singular form.

I think you have an older version. Always good practice to post what version num you're using so we can trouble shoot. Also try staxids it's not exactly what you're after but might help you get there.

• 0 views • link
written by Mark Mark

Hi Amar! I finally got it to work for me by downloading the accession2taxid file from the NCBI FTP site, and then creating a new blast database with the taxids using the option -taxid_map following this post: Scientific Names In Blast Output And Databases

• 0 views • link
written by brianna.flynn • brianna.flynn

Log in to answer this question.

More posts like this