python 进化树_科学网-应用Raxml软件构建最大似然法系统进化树-王翠的博文

各种组学数据红的热火朝天，学点生物信息傍傍身已经成了当代生物学老少的必经之路。虽然信息时代的软件层出不穷，大部分软件还没能达到兼容各种格式的程度，并且运行过程中定然会有不少坑，做个小分析简直如同去西天取经，时间都花在打怪上了。

所以，为了防患于未然，我们必须要在开始跑程序的第一天就把全网搜个遍，不然半个月以后发现有个小指令忘了加，只能自己偷着哭了。今天就扒一扒raxml这个最简单的建系统进化树的软件，其中自然也少不了步步为营。

我用的是radseq序列，经过stacks前期处理后获得vcf文件。要想得到一棵自带bootstrap值的树不太容易，废话不多说，上干货。

1，raxml要求输入文件未phylip格式，所以第一步要将vcf转换成phylip文件格式。这一步可以通过python脚本vcf2phylip来完成。https://github.com/edgardomortiz/vcf2phylip

运行结束后会产生一个名为input.min4.phy的文件。

命令行：python vcf2phylip.py -i input.vcf -o OG

2.跑raxml，应用快速bootstrap分析找出得分最高的树，关键 -f a

raxmlHPC-PTHREADS-AVX -f a -s populations.snps.sorted.min4.phy -x 12345 -N 100 -m GTRCAT -n Phylo_ml -T 28 -p 5694

得到的RAxML_bipartions.[run_name]就是带有bootstrap值的最大似然树了。

3.如果不小心跑错了，用了-f b，也有补救的办法，就看这个答案吧https://www.biostars.org/p/73678/

转载本文请联系原作者获取授权，同时请注明本文来自王翠科学网博客。

分享到:

python 进化树_科学网-应用Raxml软件构建最大似然法系统进化树-王翠的博文相关推荐

cie1931 python绘制_科学网—gnuplot与CIE1931 XYZ三刺激值曲线 - 范学良的博文
1.直接生成PNG文件 gunplot script: gnuplot> set term pngcairo Terminal type set to 'pngcairo' Options ar ...
cie1931 python绘制_科学网-gnuplot与CIE1931 XYZ三刺激值曲线-范学良的博文
1.直接生成PNG文件 gunplot script: gnuplot> set term pngcairo Terminal type set to 'pngcairo' Options ar ...
python序列_科学网—Python:序列(字符串、列表、元组)和序列函数 - 刘洋洋的博文...
Python中的序列,包括字符串(String).列表(List).元组(Tuple). 序列的索引通过索引(index)访问及获得的序列的一个或多个元素,也叫切片. 正序: 0 到 N-1 倒序: ...
python 网络_科学网-python 社会网络分析工具之networkx-郗强的博文
1.networkx 2.igraph 3.SNAP 1.networkx NetworkX是一个用Python语言开发的图论与复杂网络建模工具,内置了常用的图与复杂网络分析算法,可以方便的进行复杂网 ...
bam文件读取_科学网—Pacbio Sequel两种bam文件解析 - 卢锐的博文
pacbio目前有两种主流的测序平台,RSII和Sequel,后者是前者的升级版. pacbio sequel下机是bam格式的reads文件,它和reads比对到参考基因组上生成的bam文件,内容有 ...
matlab2014的m文件画波形,科学网—用MATLAB软件绘制驻波的波形图 - 李金磊的博文...
已知驻波的运动学方程为 y=2Acos(2πx/λ)cos(ωt) 相应的MATLAB程序为 syms lambda omega; y=2*A.*cos(2*pi*x./lambda).*cos(om ...
mysql 删除数据库的所有表_科学网—mysql删除数据库中的所有表 - 高雪峰的博文...
在工作中使用的一些数据库基本操作,不是经常使用,容易忘记,放到这里做个笔记 -- 从结果集中拷贝--------------------------------------- ------------ ...
玻尔兹曼分布涨落_科学网—高分子统计物理漫谈-涨落耗散定理-2 - 苗兵的博文...
涨落耗散定理是平衡态统计物理的一个极为重要的结果,该定理将一个统计力学体系的涨落关联与其对外界刺激的响应用一个干净的等式联系了起来.有了该定理,测量响应就可以获得平衡态体系的关联性质,反之亦然.如果取 ...
linuxcnc 多轴联动_科学网—【LinuxCNC】开源机床Linux操作系统 - 刘洋的博文
LinuxCNC是昨天意外发现的一款开源实时Linux机床操作系统,实际上也适合做机器人.机器手的操作系统:重点是[实时RealTime Linux系统].官方网站: 这款系统运行于X86结构的CPU ...
火蚁机器人_科学网—红火蚁的“生命之舟” - 赵序茅的博文
红火蚁的"生命之舟" 蚂蚁在人们眼中总是一副卑微的形象,对于人类而言"碾死一只蚂蚁"无疑是再简单不过的事情了.没有人会想到小小的蚂蚁还能兴风作浪.直到20世纪3 ...

python 进化树_科学网-应用Raxml软件构建最大似然法系统进化树-王翠的博文

python 进化树_科学网-应用Raxml软件构建最大似然法系统进化树-王翠的博文相关推荐

最新文章

热门文章