Question: BLAST Database error: No alias or index file found for nucleotide database
3
gravatar for arronslacey
4.6 years ago by
arronslacey240
United Kingdom
arronslacey240 wrote:

Hi - I am trying to run blastn, and then also SIFT standalone. I am having database configuration issues however as I am getting the following:

arron@arron-Ideapad-Z570 ~/Phd/programs/sift4.0.3b $ blastn -query test/lacI.fasta -db db/swissprot/
BLAST Database error: No alias or index file found for nucleotide database [db/swissprot/] in search path [/home/arron/Phd/programs/sift4.0.3b:::]

After some advice from other threads, I downloaded a protein database, for example swissprot:

wget ftp://ftp.ebi.ac.uk/pub/databases/fastafiles/uniprot/uniprotkb_swissprot.gz
zcat uniprotkb_swissprot.gz | awk '{if (/^>/) { print ">" $2} else { print $_}}' > swissprot.fa

and then used makeblastdb to create a blast database:

arron@arron-Ideapad-Z570 ~/Phd/programs/sift4.0.3b/db/swissprot $ makeblastdb -in swissprot.fa -dbtype prot

Building a new DB, current time: 10/27/2014 13:18:57
New DB name:   swissprot.fa
New DB title:  swissprot.fa
Sequence type: Protein
Keep Linkouts: T
Keep MBits: T
Maximum file size: 1073741824B
Adding sequences from FASTA; added 546439 sequences in 19.0039 seconds.

yet I am still getting the same problem. What am I doing wrong?

blastn blast sift makeblastdb • 32k views
ADD COMMENTlink modified 10 weeks ago by RamRS22k • written 4.6 years ago by arronslacey240
6
gravatar for RamRS
4.6 years ago by
RamRS22k
Houston, TX
RamRS22k wrote:

Once you create the database using makeblastdb, your database name is swissprot.fa

Your command should be:

blastn -query test/lacI.fasta -db <path_to_db>/swissprot.fa

I'd recommend using a custom name with makeblastdb that you can then supply to the -db parameter - that should avoid confusion with existing file-names.

ADD COMMENTlink modified 10 weeks ago • written 4.6 years ago by RamRS22k
1

hi RamRS - I actually think I should have been using blastp as swissprot as a protein database...yes it works fine with blastp

ADD REPLYlink written 4.6 years ago by arronslacey240
1

It would seem I have run out of coffee once again - should've seen that bit from a mile away. Or probably because it's Monday. 

Anyway, glad you found the solution :-)

ADD REPLYlink written 4.6 years ago by RamRS22k
1

This was actually very helpful for me. Thank you!! (Even though you were missing out on coffee). 

ADD REPLYlink written 4.0 years ago by thorerges60
1

Oh well, another coffee-less Monday, another shot in the dark!

ADD REPLYlink written 4.0 years ago by RamRS22k
0
gravatar for mforthman
3.0 years ago by
mforthman30
mforthman30 wrote:

I'm experiencing a similar issue.

Using a set of genomes/transcriptomes, I invoked makeblastdb as follows:

makeblastdb -in exons_for_blastall.fasta -input_type fasta -dbtype nucl \
-title exons_for_blastall -parse_seqids -out exons_for_blastall

This is what is reported on screen:

Building a new DB, current time: 06/01/2016 21:31:34
New DB name:   exons_for_blastall
New DB title:  exons_for_blastall
Sequence type: Nucleotide
Keep Linkouts: T
Keep MBits: T
Maximum file size: 1000000000B
Adding sequences from FASTA; added 231417 sequences in 52.605 seconds.

Next, I tried to blastn this db with the .fasta file I used to generate it (all vs. all blast):

blastn -db exons_for_blastall -query exons_for_blastall.fasta -perc_identity 80\
 -out exons_for_blastall.blasted -evalue 1e-20

You'll see I used the db name exactly. This resulted in the same error message that started this thread. Interestingly, makeblastdb created a number of files, for which I have no idea what they are for: .nhr, .nin, .nog, .nsd, .nsi, and .nsq.

ADD COMMENTlink modified 3.0 years ago by RamRS22k • written 3.0 years ago by mforthman30

I am not sure why people always choose the same name as their fasta file as the basename for the database.

Are all the ".n" files in the same directory? Could you try to create a new db in a different directory with a different name and see if that helps.

ADD REPLYlink written 3.0 years ago by genomax68k

Figure out the problem. For some reason, blastn doesn't like spaces in paths to a file, even if the appropriate symbol to indicate a space in a directory name is used.

ADD REPLYlink modified 3.0 years ago • written 3.0 years ago by mforthman30
0
gravatar for pazmonvar
10 weeks ago by
pazmonvar0
pazmonvar0 wrote:

Hello! I'm trying to parse the database UNITE to BLAST and I am getting this same error, but actually, not all .n files are generating (just the .nhr, .nin, .nsq).

This is the command I am using:

makeblastdb -in sh_general_release_dynamic_02.02.2019.fasta -input_type fasta -dbtype nucl -parse_seqids -title uniteDB2019 -hash_index -out uniteDB2019o

Is there something missing?

ADD COMMENTlink modified 10 weeks ago by RamRS22k • written 10 weeks ago by pazmonvar0

Please open a new question and reference this post in there. Do not add answers unless you're answering the top-level question.

ADD REPLYlink written 10 weeks ago by RamRS22k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1418 users visited in the last hour