This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Refseq genomic BLAST database

Hi everyone, I'd like to use a preformatted bacterial database on NCBI to run blastn searches and used the update_blastdb script to download the refseq_genomic database. All of the downloaded files are listed as refseq_genomic.## where ## is 04, or 07, 195 etc. I would like to know how I can find which of these files is my bacteria of interest, such as Streptococcus, Shewanella etc? Thanks in advance

ncbi blast database
• 4,024 views • link
updated by GenoMax GenoMax • written by chland • chland

thanks, So I would have to point to the directory as refseq_genomic.*.tar.gz for both blastn searches and extracting sequence using blastdbcmd?

• 1 views • link
written by chland • chland

Please use ADD REPLY/ADD COMMENT when responding to existing posts to keep threads logically organized.

• 1 views • link
written by GenoMax GenoMax

2 answers

I don't think that is even possible.

In any case what you want (or need to do) is to use all parts for your blast searches. They all together form a single DB, you have to use it as -db refseq_genomic in your blastcmdline (so omitting the .## in the name)

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

Use the answer provided by @5heikki in this thread to download the genomes: How to download COMPLETE bacterial genomes from NCBI based on list of names? Then index them and blast away.

You can find the names in this file (I am only looking for reference genomes).

wget ftp://ftp.ncbi.nlm.nih.gov/genomes/refseq/assembly_summary_refseq.txt
grep "reference" assembly_summary_refseq.txt | awk -F '\t' '{print $8}' > names
• 1 views • link
written by GenoMax GenoMax

Log in to answer this question.

More posts like this