phyml:基于最大似然法构建进化树

欢迎关注”生信修炼手册”!

phyml 是基于最大似然法原理构建系统发生树的软件,官网如下

http://www.atgc-montpellier.fr/phyml/

官网提供了在线服务,截图如下


共分成了四大部分

1. Input Data

输入文件为多序列比对的结果,支持以下两种格式

  1. phylip interleaved

  2. phylip sequential

这两种格式的文件都可以有 muscle 产生, 代码如下
phylip interleaved

muscle -phys   -in input.fa -out output.phys

phylip sequential

muscle -phyi   -in input.fa -out output.phyi

需要注意的是,muscle默认输出的phylip格式不能满足phyml的要求,需要进行调整,把序列合并成一行就可以了。其他的多序列比对软件,mafft 只支持输出fasta和clustalw格式的多序列比对结果,clustal 可以产生phylip 格式的文件。

目前存在的问题是,缺少多序列比对格式转换的脚本。如果想要全面支持mafft, clustal, muscle的结果,需要自己编写脚本修改结果文件的格式。

2. Subsitution Model

正是由于碱基或者氨基酸的替换,才导致了不同物种的演化,保证了物种的多样性。在构建进化树时,需要选择合适的替换模型,保证进化树的准确性。 采用默认值即可。

3. Tree Searching

采用迭代的方式,不断优化树的结构。 采用默认值即可。

4. Branch Support

进化树中的分支长度代表了不同物种的进化距离,这部分采用不同算法评估进化树中每个分支长度的可靠性。通常情况下,会选择bootstrap

也可以下载到本地运行,安装过程如下

wget http://www.atgc-montpellier.fr/download/binaries/phyml/PhyML-3.1.zip
unzip PhyML-3.1.zip

采用的是命令行交互式运行的方式,在命令行输入对应的程序名称,后续步骤和在线服务类似,也是分成了4个步骤。每个步骤之间通过+键进行确认,最后通过Y键运行。

默认生成的tree 文件是 Newick格式, 可以导入 figTree 或者 TreeViewer等软件中进行查看。

·end·

—如果喜欢,快分享给你的朋友们吧—

扫描关注微信号,更多精彩内容等着你!

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值