SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



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 09:04 AM
BLAST Database error: No alias or index file found for nucleotide database [nt] chris_s Bioinformatics 14 03-24-2019 08:11 AM
BLAST Database error: No alias or index file found for protein database gt259 Bioinformatics 2 05-25-2015 12:31 AM
How working with sam format and blast? Guigra Bioinformatics 7 04-30-2013 06:59 PM
BLAST database error - when changing to new BLAST+ local program biobio Bioinformatics 4 06-15-2011 06:20 AM

Reply
 
Thread Tools
Old 11-26-2021, 12:55 AM   #1
najoshi
Member
 
Location: Davis, CA

Join Date: Feb 2010
Posts: 16
Default 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?
najoshi is offline   Reply With Quote
Old 11-29-2021, 06:01 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,138
Default

nr is a protein database. You are trying to use blastn with that. You will need to use blastx or a translated protein search.
GenoMax is offline   Reply With Quote
Old 11-29-2021, 02:36 PM   #3
najoshi
Member
 
Location: Davis, CA

Join Date: Feb 2010
Posts: 16
Default

There are both an nr nucleotide and nr protein databases. I'm using the nr nucleotide database.
najoshi is offline   Reply With Quote
Old 11-29-2021, 02:44 PM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,138
Default

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 04:44 PM.
GenoMax is offline   Reply With Quote
Old 11-29-2021, 03:18 PM   #5
najoshi
Member
 
Location: Davis, CA

Join Date: Feb 2010
Posts: 16
Default

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 03:20 PM.
najoshi is offline   Reply With Quote
Old 11-30-2021, 03:02 AM   #6
najoshi
Member
 
Location: Davis, CA

Join Date: Feb 2010
Posts: 16
Default

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! ?!?!
najoshi is offline   Reply With Quote
Old 11-30-2021, 04:36 PM   #7
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,138
Default

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.
GenoMax is offline   Reply With Quote
Reply

Tags
blast

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 01:29 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2022, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO