一键!构建最大似然树~ 简单又准确

图片

写在前面

很多人都知道“建树”,NJ树,最大似然树 - ML,贝叶斯法建树~ 建树的目的常常有两个:

  • 分类

  • 推断演化顺序

  • 其他 - 如分析基因家族扩张收缩等

很多人都会用 mega 来建树,因为他非常方便,导入序列,点击多序列比对,然后就建树(mega默认在建树前可以选择按有效位点比例剪切)。事实上,中间还缺一个关键步骤,我相信很多不接触演化的朋友完全不知道。那就是“筛选合适的氨基酸替换模型”。IQ-tree的横空出世,解决了这一问题。软件的流行,不仅仅在于其快慢,还在于其方便程度。
当然,要做最方便,那么最好是,多序列比对 和 比对结果修剪 都直接自动化起来。对于专门做演化,尤其是质体基因组的朋友,如果需要界面化工具,那么我一般推荐的是老铁张东开发的“Phylosuite”。当然,如果是 TBtools 用户,又想偷个懒,快速画个进化树看看,那么可以用今天推出的最新功能“One Step Build a ML Tree”。

一键构建进化树

一直以来,我不太想写这个功能,一者是觉得没必要,二者是懒得动。但前段时间在演示 TBtools 的使用时,我发现绝大多数人,其实跟我一样,能偷懒的地方也像偷懒。我们还是不够勤奋,无妨。有限的时间,还是要用在点上才好,于是我开放了早就写好的一键建树功能,大体如下:

  1. 输入蛋白序列(无需比对)

  2. 输出文件路径

  3. 点击开始

  4. 等待即可
    期间,该功能会调用 Muscl 进行多序列比对,使用 trimAI 修剪比对结果,最后调用 IQ-tree 自动筛选氨基酸替换模型,然后建一个ML树。完成了,就输出进化树文本,并弹窗一个进化树,可以快速分析。

使用示例

打开 TBtools,先喝一碗主界面的鸡汤,然后选择到功能

图片

总之,打开界面之后,几乎只需要做两个事情,把蛋白序列放进去,然后设置输出文件,点击 Start 即可。

图片

如果设置直接文本输出,那么结果如下

图片

建树完成会自动弹出一个绘制好的进化树

图片

写在最后

Emmm,简单的工具,其实不需要说明书。

  • 3
    点赞
  • 23
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值