Command line 爆炸+;本地配置:如何配置nt和nr数据库?

Command line 爆炸+;本地配置:如何配置nt和nr数据库?,command-line,path,bioinformatics,blast,ncbi,Command Line,Path,Bioinformatics,Blast,Ncbi,我正在mac上配置Blast+(操作系统sierra),在配置本地下载的nr和nt数据库时遇到问题。我正试图按照NCBI的指示进行操作,并且在配置和示例执行步骤上遇到了问题 他们说要更改我的.bash_配置文件,使其显示: export PATH=$PATH:$HOME/Documents/Luke/Research/Pedulla\ 17-18/blast/ncbi-blast-2.6.0+/bin 这很好,他们说为BLASTDB“类似地”配置一个路径,但要配置到我的DB所在的文件,所以我已

我正在mac上配置Blast+(操作系统sierra),在配置本地下载的nr和nt数据库时遇到问题。我正试图按照NCBI的指示进行操作,并且在配置和示例执行步骤上遇到了问题

他们说要更改我的.bash_配置文件,使其显示:

export PATH=$PATH:$HOME/Documents/Luke/Research/Pedulla\ 17-18/blast/ncbi-blast-2.6.0+/bin
这很好,他们说为BLASTDB“类似地”配置一个路径,但要配置到我的DB所在的文件,所以我已经这样做了:

export BLASTDB=$BLASTDB:$HOME/Documents/Luke/Research/Pedulla\ 17-18/blast/blastdb/nt.00
它指定了我从他们的FTP解压nt tar文件时得到的确切文件夹。使用此路径,如果我运行命令

blastn -query test_query.fa -db nt.00 -task blastn -outfmt "7 qseqid sseqid evalue bitscore" -max_target_seqs 5
然后它成功运行,我得到了结果,但我担心这些结果只会被检查整个nt.00数据库文件的nt.00部分,特别是因为如果我在Web Blast上运行test_query.fa序列,我会得到不同的结果

另外,他们的说明说,路径只需要指向包含整个数据库文件夹nt.00的文件夹,从我解压缩的tar中——而不是特定的nt.00本身——在我的例子中,它只是“blastdb/”(与“blastdb/nt.00/”相反,后者包含nt.00.nhd、nt.00.nal等)。这是有意义的,因为当我工作时,我希望通过更改命令上的-db标志,能够在nt数据库上blastn,但也能在nr数据库上blastp,等等,并且将它们全部放在这个文件夹中应该不会有问题,对吗?但是,如果我必须指定BLASTDB的路径,并在末尾添加nt.00db,那么我怎么能在同一个文件夹(BLASTDB/)中使用nr.00呢基本上,我想按照说明所说的去做,就这样:

export BLASTDB=$BLASTDB:$HOME/Documents/Luke/Research/Pedulla\ 17-18/blast/blastdb/
然后,根据我想要使用的数据库,我可以在命令的-db标志后这样说。但是当我像上面那样做路径时,它给了我这个错误:

BLAST Database error: No alias or index file found for nucleotide database [nt] in search path [/Users/LJStout::/Users/LJStout/Documents/Luke/Research/Pedulla 17-18/blast/blastdb:]
我尝试过从上面运行同一个blastn命令,并将“nt”替换为“nt.00”,并且尝试过这些命令,BLASTDB的路径以“BLASTDB/”和“BLASTDB/nt”结尾,当然还有“BLASTDB/nt.00”,这是唯一一个运行没有错误的命令

在我读到的另一篇文章中,OP担心他的执行没有检查整个nt.00文件夹,但这与我的问题不同


谢谢你的帮助

整个问题归结为nt.00&nr.00文件,即解压其各自的.tar.gz后产生的原始文件夹,在相同的父文件夹中,而它们的内容应该在相同的父文件夹中。我只是删除了它们进入的文件夹,并将内容复制到我的新文件夹中,单亲父母。我被指示误导了,这是个简单的错误。现在,我有一个文件夹,
blastdb/
,其中包含我计划使用的每个数据库的所有内容,包括nt、nr和refseq

为您提供信息,这里还有一个生物信息学stackexchange网站:谢谢,我错贴在这里了。