magicblast aganist local nr database
I have downloaded the complete nr datadase to my local server. I have now a list of nr.00* to nr.38*
nr.00.phd
nr.00.phi
nr.00.phr
nr.00.pin
nr.00.pog
nr.00.ppd
nr.00.ppi
nr.00.psq
nr.00.tar.gz.md5
nr.01.phd
...
nr.38.ppd
nr.38.ppi
nr.38.psq
nr.38.tar.gz.md5
nr.pal
nr.pdb
nr.pos
nr.pot
nr.ptf
nr.pto
taxdb.btd
taxdb.bti
I would like now to run magicblast against this local data base, but I keep getting the error massage, I don't have an index file.This is the command I am trying to run:
$magicblast -query R2_val_2.tagged_filter.fastq.gz -db /home/BLASTdb/nr -infmt fastq -outfmt tabular -num_threads 5 -out WK56.R2.unmapped.magicblast.nr.txt
I am getting this error
BLAST Database error: No alias or index file found for nucleotide database [/home/BLASTdb/nr] in search path [/local/projects/Chad/P286::]
Can someone tell me what I am doing wrong? Do I need to somehow combine these files?
thanks
• 2,574 views
•
link
written
by
Assa Yeroslaviz
1 answer
yes it does.
my command was as such
magicblast -query tagged_filter.fastq.gz -db BLASTdb/nt/nt -infmt fastq -outfmt tabular -num_threads 5 -out unmapped.magicblast.nt.txt
while the folder with the data base is as such
$ ls -1 BLASTdb/nt
nt.00.nhd
nt.00.nhi
nt.00.nhr
nt.00.nin
nt.00.nnd
nt.00.nni
nt.00.nog
nt.00.nsq
nt.00.tar.gz.md5
nt.01.nhd
nt.01.nhi
...
nt.27.tar.gz.md5
nt.nal
nt.ndb
nt.nos
nt.not
nt.ntf
nt.nto
taxdb.btd
taxdb.bti
• 0 views
•
link
written
by
Assa Yeroslaviz
Log in to answer this question.
Hello Assa Yeroslaviz!
We believe that this post does not fit the main topic of this site.
I have downloaded the nr database, which is the non redundant protein database, but I need nt!
For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.
If you disagree please tell us why in a reply below, we'll be happy to talk about it.
Cheers!
magicblastis only meant for use against DNA databases.I know, I wanted to run it against the
ntdatabase, notnr.Heed this warning from magicblast help page. Trying to map against
ntwill not be a good idea.But I don't have a FASTA file.
I have a BLAST DB, splitted into several files
hence the 00, 01, etc.
Or do I miss something?
Documentation for magicblast only refers to using a (one) genome/transcriptome (or at most a few SRA accessions). Let us know if using the
ntblast database works. As you said above you will replacenrwithntdatabase.