getting sscinames in viral blastx query
2
0
Entering edit mode
6.0 years ago

I am running a local blast instance on my computer and I am trying to extract the sscinames from the database.  I downloaded the viral database from ftp.ncbi.nlm.nih.gov/refseq/release/viral and created the blast database (with taxID) via : 

 grep -i "Viral\|virus" /reference/blast_database/refseq/RefSeq-release74.catalog &> viral_taxid.txt

awk 'BEGIN{FS= "\t"}{print $3" " $1}' viral_taxid.txt &> viral_taxid_map.txt

cat viral.2.protein.faa viral.1.protein.faa &> combined.protein.faa

makeblastdb -parse_seqids -in combined.protein.faa -out combined.protein -dbtype prot -taxid_map viral_taxid_map.txt

 

When I blast against my database, e.g.

blastx -num_threads 1 -outfmt '6 qseqid evalue bitscore staxids sgi sacc sscinames stitle' -db combined.protein -query data.fa

it does not return the sscinames.  Is there a way to incorporate this into my blast database?

 

Thanks.

 

refseq blast • 3.4k views
ADD COMMENT
7
Entering edit mode
6.0 years ago

SOLUTION : I "export BLASTDB="/path/to/taxdb" which is the directory where taxdb.btd and taxdb.bti are installed.

ADD COMMENT
0
Entering edit mode
7 weeks ago
Alex Bel ▴ 30

Another stupid but easy way to fix it -- put 'taxdb.bti' and 'taxdb.btd' files into the directory where you start your blastn command

ADD COMMENT

Login before adding your answer.

Traffic: 2255 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6