前言:Blast最新版是2.10.0,同时因为项目需要,要做NT数据库,也就是核酸数据库的本地化。(NT是核酸库,NR是蛋白质库)
- Blast下载及安装
我选择用axel下载软件的源码,不知道axel的可以翻我之前的文章,解压即可使用blast
axel https://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST/ncbi-blast-2.10.0+-x64-linux.tar.gz
tar -xf ncbi-blast-2.10.0+-x64-linux.tar.gz
cd ncbi-blast-2.10.0+/bin
./blastn -help
建议添加至环境变量
-
下载NT核酸数据库
NT数据库中,有被format好的NT数据库,共分割为25个子文件,把他们下载下来,解压就可以直接用,约65G:
NT核酸库
同样使用axel进行下载,因为axel在链接断掉后,可以自动重连:
mkdir NCBI_DB
for i in `seq 0 9`
do
axel https://ftp.ncbi.nih.gov/blast/db/nt.0${i}.tar.gz
done
for i in `seq 10 24`
do
axel https://ftp.ncbi.nih.gov/blast/db/nt.${i}.tar.gz
done
下载好的文件:
NT数据库压缩包
解压:
tar -xf *.tar.gz
解压完成后把全部的文件放到/home/xxx/blast/db
,还需要在家目录下配置.ncbirc
文件,在文件中添加内容:
[BLAST]
BLASTDB=/home/xxx/blast/db
- 本地化Blast的使用
blastn -query test.fasta -db /path/to/NCBI_DB/nt -outfmt 6
网友评论