多个Sanger测序结果和多序列比对结果的可视化

写在前面

转眼五一假期就要结束。事实上,在每一个假期,我总是有一些想做的事情,比如完成某个分析任务?或者是写一个不错的工具。今次假期比较特殊,我似乎变得愿意安排出时间陪伴本应多陪伴的人。当然,这也是TBtools早期开发的目的之一。

节省可以节省的时间,用于陪伴家人。

虽然说是陪伴,主要还是在家宅着,看了几部电影,吃了罗锅,大量水果,大量肉食...。从某个角度来说,这种休假方式,确实是比较舒适。Anyway,我一样安排了时间做一些数据分析的工作。但毕竟是假期,我的习惯可能还是make something。于是,花了大半天的时间,在TBtools中,增加了两个功能。

批量查看Sanger测序峰值图(.ab1)

做分子实验的朋友,可能才会明白这个功能存在的意义。十年前(...大一...我确实有点老了),我入学并随后进去了华南农大何业华教授的课题组学习(也是我的硕士导师)。我接触的第一个分子实验,即是克隆一个基因。但凡克隆,就要测序。对于每个克隆,公司往往会返回两个文件:

  1. .ab1 文件,即测序峰值图(下机数据,信息全面)
  2. .seq或.txt 文件, 即baseCalling之后的序列(只对可信位点取峰值最好的碱基),丢失了信息,也可能出现错误Calling,但实用

往往,我们克隆并送质粒到公司做Sanger测序时,我们总会对同一个克隆片段,送几个样品与公司。一方面,Sanger测序时部分碱基也会测不准,多送几个样品,可以相互验证;另一方面,我们需要考虑SNP存在的可能...

于是,我们会同时对同一个片段,收到数个.ab1文件。


对于这些.ab1文件,存在很多工具可以直接查看,如bioedit,DNAman?但很多时候,似乎他们只支持一次展示一个样品。而我们需要的可能是同时真实四个样品,以人工查看某些位点/SNP的可靠性(我当然知道有软件可以做到,我只是看看这个功能是不是很难写,结果似乎不是)。

使用这个功能比较简单,首先是打开这个功能



将.ab1文件直接拖拽或通过摁钮选择并放入其中(这里用两个文件做示例)



点击Start,即可看到这个图

Emm...是的,我并没有将多序列比对应用进去...感觉似乎麻烦。
所以此时,应该是键盘摁住Ctrl键,鼠标拖拽SubPanel,使其对其


别忘了,鼠标滚轮可以放大缩小图片

OK,发现一个测序错误?还是?总之多了两个AA,而且应该是比较可靠
总的来说,似乎有用...尤其是当你不止有一个克隆的情况下。或者是,你需要一张图,去说服其他人,这确实有个SNP....

快速可视化多序列比对结果

很久之前,我已经将Muscle打包到TBtools中,那会只是顺便。没想到在后续我有数次还是用上了。总的来说,体验一般。主要原因是,没有可视化结果,那么其实即使比对出来,也不容易做数据观测。
于是我加了一个自动展示多序列比对结果的功能。
首先是打开这个Wrapper


随后是设置用于多序列比对的序列(可以是文件,也可以是文本区域;输出也是一样支持两种)
简单实用以上的序列作为输入

可以同时看到图片,如果你还记得这两个AA

但,事实上,这个功能更适合序列较多的情况,如

写在后面

写工具,终归是消遣时间。睡觉。

题外

课题组每年暑期有内部生信入门培训,主要是对实验室新生开展(以及湿实验为主的成员)培训。一直有收到其他课题组想要了解我们课题组生信数据分析的想法。故,在博导的提议和课题组的讨论后,我们近期计划,在本年度暑期(7~8月份之间)对外增设生信基础培训名额10枚(前面每年只是课题组内培训,而不对外)。具体请见https://mp.weixin.qq.com/s/OtmeTErd9f9rvjJPtBKjMw

欢迎访问本课题组网站

园艺植物小分子RNA与基因组研究-夏瑞课题组

课题组主页:http://xialab.scau.edu.cn/

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