建立本地的Blast数据库
Blast(basic local alignment search tool) 局部序列比对基本检索工具,是NCBI开发的一款基于序列相似性的数据库搜索程序。主要的BLAST程序有以下几种:
BLAST的在线版https://blast.ncbi.nlm.nih.gov/Blast.cgi可以方便的进行单一的同源性序列搜索,但是不方便进行大批量的数据操作,也不能建立自己的数据库。接下来小编就要教大家如何建立本地的BLAST数据库。
基于BLAST优秀的算法,BLAST程序可以轻松的在普通的个人电脑中运行。
BLAST程序的官方下载地址:ftp://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST
在网页中我们可以看到很多文件,其中*.md5是效验文件,我们可以忽略它,现在一般是x64的操作系统,我们可以根据自己的电脑系统,下载x64的win或Linux或macosx版本。
小编下载的是windows的一个压缩包版本,解压后得到这样一个文件夹:
进去后可以看到目录结构并不复杂,readme文件对BLAST作了一个简要的介绍:
bin文件夹通常是主程序的文件夹,我们进去后看到很多*.exe的可执行程序,但是双击后并不能运行。
这种情况一般是需要用命令提示符来运行。
通过win+R 输入cmd命令,再通过cd 命令进入当前的bin 文件夹,或者直接在bin文件夹内按住shift和鼠标右键,选择在此文件夹打开命令提示符
我们通过makeblastad -help命令,查看创建本地数据库的帮助文件
可以发现需要输入文件,在NCBI数据库,我们可以方便地下载我们需要的序列文件,下载后复制到bin文件夹内,方便操作。
小编在这里下载了一个红曲基因组的fasta文件
并复制到bin文件夹下
之后执行命令 makeblastdb in sequence.fasta -dbtype mucl -out sesameESTdb
这里取名为sesameESTdb是受某教程误导,大家可以按习惯起名字。
-in 后面是空格加输入数据库文件的名称,-dbtype后面是空格加数据库类型(核苷酸是nucl,蛋白质是prot),-out后面是空格加输出数据库名称
接下来就可以拿目标序列来比对了,由于手头上没有相似的序列,小编直接从sequence.fasta中复制了一段,另存为test.fasta进行比对。
输入 blastn -query test.fasta -db sesameESTdb -out 001.txt
可以看到在bin文件夹下出现了一个001.txt文件,打开后可以查看结果:
至于结果的含义大家可以参考BLSAT官方文档
ref:
Zheng Zhang, Scott Schwartz, Lukas Wagner, and Webb Miller (2000), "A greedy algorithm for aligning DNA sequences", J Comput Biol 2000; 7(1-2):203-14.
- 【译】Spring 官方教程:使用 Restdocs 创建 API 文档
- c#中设置快捷键
- 程序员你为什么这么累【续】:编码习惯之工具类规范
- IntelliJ IDEA 复杂的重构技巧
- 打造属于自己的支持版本迭代的Asp.Net Web Api Route
- 分布式消息队列 RocketMQ 源码分析 —— Message 拉取与消费(上)
- 分布式消息队列 RocketMQ 源码分析 —— Message 拉取与消费(下)
- Spring Boot中使用RabbitMQ
- Spring Cloud构建微服务架构:消息驱动的微服务(入门)【Dalston版】
- 哪类人适合当产品经理?
- 产品经理·杂谈
- Python机器学习中的特征选择
- Android学习第八弹之改变状态栏的颜色使其与APP风格一体化
- 手把手教你 MongoDB 的安装与详细使用(二)
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- 【Java8新特性】01 函数式接口和Lambda表达式你真的会了吗
- 只出现一次的元素
- 【排序】堆排序
- redis灵魂拷问:怎样搭建一个哨兵主从集群
- 是时候给你的鼠标指针更换样式换下风格体验了!如何卸载与安装鼠标指针?
- Linux设置虚拟内存
- SSM 单体框架 - 教育平台后台管理系统:接口文档
- java学习应用篇|使用环境变量做一些工具
- 【Java8新特性】06 新的日期和时间:LocalDate LocalTime LocalDateTime
- UNSAFE和Java 内存布局
- 【玩转Redis面试第1讲】Redis数据结构和常用命令速记
- SSM 单体框架 - 教育平台后台管理系统:课程模块
- 【玩转Redis面试第2讲】面试官再问Redis事务把这篇文章扔给他
- SSM 单体框架 - 教育平台后台管理系统:Maven 进阶使用
- SSM 单体框架 - 教育平台后台管理系统:广告和用户模块开发