blast nr database not working

najoshi · 11-25-2021, 11:55 PM

So I am trying to download the latest nr database and use blast with it on Ubuntu 16.04. I compiled and installed the latest blast and then used update_blastdb.pl to download nr, which worked fine. But when I try to run a blast query using nr I get this error:

BLAST Database error: No alias or index file found for nucleotide database [/share/genomes/ncbi/nr/nr] in search path [/share/biocore/joshi/sandbox::]

And I am definitely using the correct way of specifying nr. Here is my blastn command:

blastn -query test.fa -db /share/genomes/ncbi/nr/nr -num_threads 30

/share/genomes/ncbi/nr/ is a directory containing all of the nr files.

So then I downloaded nt and ran the same command with nt and it worked just fine. It seems that whatever index file is needed by blastn is not included in any of the tar files for nr. Anybody have any ideas?

GenoMax · 11-29-2021, 05:01 AM

nr is a protein database. You are trying to use blastn with that. You will need to use blastx or a translated protein search.

najoshi · 11-29-2021, 01:36 PM

There are both an nr nucleotide and nr protein databases. I'm using the nr nucleotide database.

GenoMax · 11-29-2021, 01:44 PM

If you downloaded `nr` as pre-formatted database from NCBI then it is a protein only database as noted in the blastftp file.

najoshi · 11-29-2021, 02:18 PM

So according to the blastftp file, nt is "Not non-redundant". Where would I download the nucleotide nr? Or are the nt.00.* files actually the nr nucleotide database? Because that is confusing.

najoshi · 11-30-2021, 02:02 AM

So according to that blastftp file:

nt.00.tar.gz first volume of the nucleotide nr database
nt.01.tar.gz second volume of the nucleotide nr database

nt is nucleotide nr! And nr is protein nr! ?!?!

GenoMax · 11-30-2021, 03:36 PM

I don't think there is a non-redundant nt database. You should consider using RefSeq/RefSeq genomes.

Perhaps the only way to put this to rest is to email NCBI help desk and ask. There is some ambiguity for sure since the description of the `nt` FASTA file counterpart says

Quote:

nt.gz nucleotide database from GenBank excluding the
batch division htgs, est, gss,sts, pat divisions,
and wgs entries. Not non-redundant.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
BLAST+ creating custom blast database and using blast+ filtering features	deniz	Bioinformatics	3	07-07-2019 08:04 AM
BLAST Database error: No alias or index file found for nucleotide database [nt]	chris_s	Bioinformatics	14	03-24-2019 07:11 AM
BLAST Database error: No alias or index file found for protein database	gt259	Bioinformatics	2	05-24-2015 11:31 PM
How working with sam format and blast?	Guigra	Bioinformatics	7	04-30-2013 05:59 PM
BLAST database error - when changing to new BLAST+ local program	biobio	Bioinformatics	4	06-15-2011 05:20 AM

11-25-2021, 11:55 PM	#1
najoshi Member Location: Davis, CA Join Date: Feb 2010 Posts: 17	blast nr database not working So I am trying to download the latest nr database and use blast with it on Ubuntu 16.04. I compiled and installed the latest blast and then used update_blastdb.pl to download nr, which worked fine. But when I try to run a blast query using nr I get this error: BLAST Database error: No alias or index file found for nucleotide database [/share/genomes/ncbi/nr/nr] in search path [/share/biocore/joshi/sandbox::] And I am definitely using the correct way of specifying nr. Here is my blastn command: blastn -query test.fa -db /share/genomes/ncbi/nr/nr -num_threads 30 /share/genomes/ncbi/nr/ is a directory containing all of the nr files. So then I downloaded nt and ran the same command with nt and it worked just fine. It seems that whatever index file is needed by blastn is not included in any of the tar files for nr. Anybody have any ideas?

11-29-2021, 01:44 PM	#4
GenoMax Senior Member Location: East Coast USA Join Date: Feb 2008 Posts: 7,143	If you downloaded `nr` as pre-formatted database from NCBI then it is a protein only database as noted in the blastftp file. Last edited by GenoMax; 11-29-2021 at 03:44 PM.

11-29-2021, 02:18 PM	#5
najoshi Member Location: Davis, CA Join Date: Feb 2010 Posts: 17	So according to the blastftp file, nt is "Not non-redundant". Where would I download the nucleotide nr? Or are the nt.00.* files actually the nr nucleotide database? Because that is confusing. Last edited by najoshi; 11-29-2021 at 02:20 PM.

11-29-2021, 05:01 AM	#2
GenoMax Senior Member Location: East Coast USA Join Date: Feb 2008 Posts: 7,143	nr is a protein database. You are trying to use blastn with that. You will need to use blastx or a translated protein search.

11-29-2021, 01:36 PM	#3
najoshi Member Location: Davis, CA Join Date: Feb 2010 Posts: 17	There are both an nr nucleotide and nr protein databases. I'm using the nr nucleotide database.

11-30-2021, 02:02 AM	#6
najoshi Member Location: Davis, CA Join Date: Feb 2010 Posts: 17	So according to that blastftp file: nt.00.tar.gz first volume of the nucleotide nr database nt.01.tar.gz second volume of the nucleotide nr database nt is nucleotide nr! And nr is protein nr! ?!?!