步骤一：NT/NR库全库下载

用wget 对数据库进行下载

NT/NR库：https://ftp.ncbi.nlm.nih.gov/blast/db/

NT库有76个子文件构建成NT全库

NR库有63个子文件构建成NR全库

（ncbi网站处理好的文件，不需要建库）

https://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/nt.gz

https://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/nr.gz

这两个为nr/nt的fasta文件，需要建库

写一个简单的脚本，对nr/nt库进行下载

#!/bin/bash
for i in {00..76}
dowget -c https://ftp.ncbi.nlm.nih.gov/blast/db/nt.${i}.tar.gzwget -c https://ftp.ncbi.nlm.nih.gov/blast/db/nt.${i}.tar.gz.md5md5sum -c nt.${i}.tar.gz.md5tar -zxvf nt.${i}.tar.gz -Cecho "nt.${i} has done."
done

#!/bin/bash
for i in {00..63}
dowget -c https://ftp.ncbi.nlm.nih.gov/blast/db/nr.${i}.tar.gzwget -c https://ftp.ncbi.nlm.nih.gov/blast/db/nr.${i}.tar.gz.md5md5sum -c nr.${i}.tar.gz.md5tar -zxvf nr.${i}.tar.gz -Cecho "nr.${i} has done."
done

注意：下载的时候要加-c参数，避免因网络问题下载失败

下载完后注意删除安装包

步骤二：把序列与nt/nr库进行比对

基本用法

蛋白质序列比对蛋白质数据库（blastp）

blastp -query seq.fast -out seq.blast -db dbname -outfmt 7 -evalue 1e-5 -num_des
criptions 10 -num_threads 8

核酸序列比对核算数据库（blastn）

blastn -query seq.fasta -out seq.blast -db dbname -outfmt 7 -evalue 1e-5 -num_descriptions 10 -num_threads 8

核酸序列比对蛋白质数据库（blastx）

blastx -query seq.fasta -out seq.blast -db dbname -outfmt 7 -evalue 1e-5 -num_descriptions 10 -num_threads 8

参数：
-query：输入文件路径及文件名
-out：输出文件路径及文件名
-db：格式化了的数据库路径及数据库名
-outfmt：输出文件格式，总共有12种格式，6是tabular格式对应BLAST的m8格式
-evalue：设置输出结果的e-value值
-num_descriptions：tabular格式输出结果的条数
-num_threads：线程数

创建索引数据库

makeblastdb -in mouse.faa -dbtype nucl -input_type fasta -out dna

参数：
-in: 输入数据库文件
-dbtype:数据库类型如果是蛋白质数据库则 nucl改为 prot
-out: 索引数据库名称

创建索引文件

dustmasker -in Nuc.fasta -infmt fasta -parse_seqids -outfmt maskinfo_asn1_bin -out Nuc_all.asnb

参数：
in:输入数据库文件
-infmt:输入文件的格式
-parse_seqids:解析序列的id，这个只需要提供参数，不要提供文件名
-outfmt:输出文件的格式，直接填写即可maskinfo_asnl_bin即可。
out：输出文件的文件名（索引文件名）。

创建数据库

makeblastdb -in Nuc.fasta -input_type fasta -dbtype nucl -parse_seqids -mask_data NucDatsBase.asnb -out Nuc_all

参数：
in:原本我们数据库的文件
input_type:输入文件的格式，很明显，我这里是fasta格式

本地blast与nr/nt库相关推荐

构建NCBI本地BLAST数据库 (NR NT等) | blastx/diamond使用方法 | blast构建索引 | makeblastdb...
参考链接: FTP README 如何下载 NCBI NR NT数据库? 下载blast:ftp://ftp.ncbi.nlm.nih.gov/blast/executables/blast+ 先了解 ...
NCBI下载nt/nr/swissprot库
NCBI下载nt/nr/swissprot库 1. 确定你要下载文件的位置 https://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/ 2. 进行下载方法一:使用wge ...
生物信息学习--nr/nt 数据库（总+子）构建
1. 从ncbi上下载数据下载地址:https://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/ mac端: 安装brew,运行如下命令: /bin/zsh -c &qu ...
本地blast详细用法
本地Blast的详细使用方法 blastall -p blastn -i myRNA.fasta -d humanRNA.fasta -o myresult.blastout -a 2 -F F -T ...
NCBI的NT库比对——blastn
NCBI的NT库比对--blastn NT库比对,包括对测序数据和组装后的基因组序列进行NT库比对,查看是否存在菌污染以及是否是自己的数据.这里我提供这一部分的具体操作过程. 步骤一:NT全库下载前 ...
生成mysql nt文件_linux下安装blast并创建nt数据库
BLAST 是由美国国立生物技术信息中心(NCBI)开发的一个基于序列相似性的数据库搜索程序.BLAST是"局部相似性基本查询工具"(Basic Local Alignment S ...
生信分析-本地BLAST
一. 本地blast简介本地Blast(Basic Local Alignment Search Tool),是基于本地的比对搜索工具,可以在自己建立的数据库进行blast搜索,与NCBI的在线bl ...
基因家族的鉴定-基于windows系统上的本地blast
文章首发于简书链接(https://www.jianshu.com/p/c876a5636554).发此备份. 基因家族的鉴定-基于windows系统上的本地blast 基因组的序列提取,详情请看我之 ...
本地BLAST的使用方法及基本操作步骤
文章目录 BLAST+程序下载与安装基因组数据下载制作数据库 BLAST 叮! BLAST+程序下载与安装 BLAST+安装包下载地址:ftp://ftp.ncbi.nlm.nih.gov/bla ...

本地blast与nr/nt库

步骤一：NT/NR库全库下载

步骤二：把序列与nt/nr库进行比对

本地blast与nr/nt库相关推荐

最新文章

热门文章