使用TBtools进行本地blast并根据blast结果提取相关序列

<meta charset="utf-8">

TBtools 是一个集合了超多生信处理的小工具,非常方便,而且是无命令行的操作,入门门槛极低。开发者是华南农业大学陈程杰博士,TBtools文章地址:https://www.biorxiv.org/content/early/2018/03/27/289660
,使用参考页面https://www.yuque.com/cjchen/hirv8i
最近要使用本地blast来寻找组内基因组的相似序列,所以准备使用tbtools来试试。

1.blast使用

1.1 数据库构建

参考网页https://zhuanlan.zhihu.com/p/418163788
首先打开软件界面,选择BLAST→BLAST GUI Wrapper→Blast Zone 会出现如下界面

image.png

然后在这个界面先构建一个本地的数据库,作为blast的数据库。
image.png

可以点击config按钮来设置数据库的存放位置。
image.png

点击图中的加号就可以添加本地的数据来构建数据库了。
image.png

点击确定之后一段时间完成数据库构建
image.png

如果要添加更多的数据库操作类似,删除数据库的话直接选中数据库然后选择 “-”就行了。
image.png

image.png

1.2 使用序列进行BLAST

直接将需要BLAST的序列粘贴到input文本框中,选择好输出文件的格式、文件名和位置点击start就可以开始了(超短序列需要选择short mode)。


image.png

查看结果


image.png

结束之后,点击右下角的 Visulize 或TextView进行查看,Visulize有几种模式查看,点击之后在左上角出现图标选项,下拉选择之后再关掉就会出现显示的界面。
如选择Aligmemt Shorer会出现如下的结果
image.png

选择Dot Plot会出现如下结果:



选择Pileup Graph会出现如下结果

而使用TextView则会显示文本信息如下:
image.png

(Frame = -1则表示是反向序列)

这样就能看到BLAST所有结果的信息。

2. 根据Blast结果提取对应序列

提取对应序列需要使用tbtools中的Fasta Extract功能。


image.png

点击Sequence→Fasta tools→Fasta Extract (Recommended)进人Fasta Extract界面


image.png

从界面上来看,分为 6 大块,其中绝大多数是可选项(即可以不做调整)。

  1. Fasta 序列文件输入文本框,用户可以直接拖拽硬盘中的 Fasta 文件并放置到文本框中,路径会自动获取;也可以点击跟随文本框的摁钮“...”,在弹出文件选择框中选取对应文件即可
  2. Initialize 摁钮,在设置 Fasta 序列文件后,可以看到 Start 摁钮仍然不可点击。需要用户点击 Initialize 摁钮,创建 Fasta 序列索引文件(如前期已有,则会软件会自动复用,节省计算时间)
  3. 输出文件设置文本框,用户同样可以拖拽放置文件或者文件夹,程序会自动获取输出文件夹,用户需要补全一个输出文件名;当然也可以直接点击跟随文本框的摁钮,在弹出的问价选择框中设置对应输出文件即可
  4. 待提取序列信息设置框,参考界面说明,接受三种类型的提取模式:
    ● 基于 ID 的完整序列记录提取,如输入 Unigene_1 ... 等完整序列 ID,每行一个,即可提取完整序列记录
    ● 基于 序列坐标信息,进行序列区间截取,如提取染色体 Chr1 上第 10000 个碱基到 20000 个碱基的一段序列,那么输入如下。如果需要提取反向互补序列,使起始坐标大于终止坐标即可。

# 注意,制表符[\t]分隔,而非空白[Space]分隔(试了下现在用空格好像也可以了)
Chr1 10000 20000
# 提取反向互补序列,则翻转碱基坐标
Chr1 20000 10000

提取坐标信息的功能,重命名区间,如我们需要提取 Peak 或者 Promoter 序列信息,并指定输出时序列名字

peak_1 Chr1 10000 10200
promoter_ATG8 Chr2 20300 22300

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 216,402评论 6 499
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,377评论 3 392
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 162,483评论 0 353
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,165评论 1 292
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,176评论 6 388
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,146评论 1 297
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,032评论 3 417
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,896评论 0 274
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,311评论 1 310
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,536评论 2 332
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,696评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,413评论 5 343
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,008评论 3 325
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,659评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,815评论 1 269
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,698评论 2 368
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,592评论 2 353

推荐阅读更多精彩内容