于平沙如雪,构广厦万千|

鹿衔草_Yusy

园龄:2年2个月粉丝:1关注:0

使用Bioaider进行本地blast

系统环境为windows11

1. 下载blast程序

https://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST/

image.png

双击安装,记住自己的安装目录

image.png

2. 在Bioaider里配置本地Blast程序路径

image.png

image.png

选择到bin目录,add即可,此时Status会变成绿色的Installed

image.png

image.png

3. 新建数据库

首先就是新建一个数据库,这个数据库存放的是你待会儿要和他们进行比对的序列,可以理解为NCBI提供的nr库、nt库之类的

image.png

点击左下角Import database,导入数据库,然后输入你准备好的一堆序列(如果是多个多个文件夹,切记要先进行合并),在Linux下合并只需通过一条简单的命令即可实现

cat *.fas > H5_gisaid.fas

image.png

image.png

然后点击Add,等待建库完成即可

image.png

看到Success提示后我们再点开DBs,就可以看到当前多了一个H5_gisaid,此时把你的序列放在左边的Query Sequence里边,然后设置一下Outdir路径即可。至于Parameter里的参数,相信大家都知道如何设置。

image.png

然后我们就可以在对应的路径下看到txt格式的比对结果了,但是这个结果是没有表头的,表头从左到右依次是

Query id Subject id % identity alignment length mismatches gap openings q. start q. end s. start s. end e-value bit score

顺便附上表头数据的说明

1、Query id:查询序列ID标识
2、Subject id:比对上的目标序列ID标识
3、% identity:序列比对的一致性百分比
4、alignment length:符合比对的比对区域的长度
5、mismatches:比对区域的错配数
6、gap openings:比对区域的gap数目
7、q. start:比对区域在查询序列(Query id)上的起始位点
8、q. end:比对区域在查询序列(Query id)上的终止位点
9、s. start:比对区域在目标序列(Subject id)上的起始位点
10、s. end:比对区域在目标序列(Subject id)上的终止位点
11、e-value:比对结果的期望值
12、bit score:比对结果的bit score值

image.png

大功告成,离开了命令行感觉非常舒适丝滑,可能这就是GUI的意义吧。

本文作者:Yusy's Blog

本文链接:https://www.cnblogs.com/yusy/p/18725093

版权声明:本作品采用知识共享署名-非商业性使用-禁止演绎 2.5 中国大陆许可协议进行许可。

posted @   鹿衔草_Yusy  阅读(17)  评论(0编辑  收藏  举报
点击右上角即可分享
微信分享提示
评论
收藏
关注
推荐
深色
回顶
收起