Dr_W的个人博客分享 http://blog.sciencenet.cn/u/Dr_W

博文

linux下blast+的使用

已有 12959 次阅读 2015-4-9 21:36 |系统分类:科研笔记|关键词:学者


   按照网上的教程,经过高人指点,稍微进行修改, 终于运行成功。现总结出来,以备后阅。

安装

  我用的系统是ubuntu系统,blast+的安装可以使用命令sudo apt-get install ncbi-blast+进行安装。

运行实例:

#建立blast文件夹

mkdir blast

#进入blast文件夹

cd blast

#下载序列

wgetftp://ftp.ncbi.nih.gov/genomes/Arabidopsis_thaliana/CHR_I/NC_003070.fna

#修改文件名

mv NC_003070.fna Arabidopsis_thaliana.fna

#取前100000行建立新文件sample.01.fna

head-n 100000 Arabidopsis_thaliana.fna > sample.01.fna

#格式化数据库

makeblastdb-in Arabidopsis_thaliana.fna -out Arabidopsis_thaliana -dbtype nucl

In参数后面接将要格式化的数据库,

--dbtype 告诉程序这是核酸数据库。

#进行blastn比对

blastn -query sample.01.fna -db Arabidopsis_thaliana -out test.out -evalue0.001

query :查询序列,必须是fasta格式;

-db :格式化好的查询数据库;

-out :结果输出文件。

#查看结果

less test.out|more

成功




https://m.sciencenet.cn/blog-621869-881175.html

上一篇:真难
下一篇:如果是真的,那这个期刊可真是

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-23 21:54

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部