linux下phylip软件构建NJ树,SNP数据构建系统进化树

大概是长期的不锻炼使得今天的爬山运动过量了，接着悲剧就是无法入眠。也幸亏明天是周日，干脆就起床码字了。

总结下自己前面用snp构建系统进化树的方法吧。

1.构建进化树的算法

构建系统进化树的方法主要有以下几类：

基于距离矩阵的方法:NJ(邻接法)

MP(最大简约法)

ML(最大似然法)

以及贝叶斯法。

一般情况下，若有合适的模型，ML的效果较好；

近缘序列的话，一般使用MP；

远源序列，一般使用NJ或者ML。

在分析变异过滤得到SNP时，一般都会用PHYLIP构建NJ进化树。

那么具体如何操作呢？

软件安装

软件安装较为简单

wget http://evolution.gs.washington.edu/phylip/download/phylip-3.69.tar.gz ./ #下载软件

tar zxvf phylip-3.69.tar.gz #解压

cd phylip-3.69/src

make install

#以上几步即安装完软件，文件夹中的exe目录里为可执行程序

2.输入文件的格式转换

输入文件格式事例.png

其中第一行为构建进化树的样品数以及每个样品使用的snp数目。

第二行及以下为每个样品的名称及snp的具体内容。需要注意的是样品的名称必须为10个字母，如果未达到10个字母，可用tab键或者空格键代替。第11个字母后即为snp的内容，同时在这些序列中，一般每10个位点会有1个空格使其方便阅读。每个样品的用于构建snp的个数必须相同。

根据以上的规定，可以写脚本将vcf格式转化为可用于phylip的phy格式。

3.软件使用

phylip中有许多程序，大部分的程序运行方法相同，把infile作为默认的输入文件，输出结果写在outfile中。因此，在进行下一步分析前，需要重命名想要保存的文件。

seqboot: 生成随机样本，用bootstrap和jack-knife方法。需要设置选项M

dnadist：DNA距离矩阵计算器。

neighbor：NJ法的使用

consense：用多重树构建一致树。

每个程序都需要设定参数，因此还需要新建par文件。

#cat seqboot.par

all.merge.snp.phy #设定输入文件的名称，否则输入默认的名为infile的文件

r #选择bootstrap

1000 #设置bootstrap的值，即重复的replicate的数目，通常使用1000或者100，注意此处设定好后，后续两步的M值也为1000或者100

y #yes确认以上设定的参数

9 #设定随机参数，输入奇数值。

#cat dnadist.par

seqboot.out #本程序的输入文件

t #选择设定Transition/transversion的比值

2.3628 #比值大小

m #修改M值

d #修改M值

1000 #设定M值大小

2 #将软件运行情况显示出来

y #确认以上设定的参数

#cat neighbor.par

dnadist.out #本程序的输入文件

1000 #设定M值大小

9 #设定随机数，输入奇数值

y #确认以上设定的参数

# cat consense.par

nei.tree #本程序的输入文件

y #确认以上设定的参数

再运行以下命令行即可

seqboot

最后将会得到constree文件，可将该文件gaiwei*.tre文件，双击后在treeview中直接查看进化树的内容。

若要进行进一步的编辑，可使用iTOL在线的网站(http://itol.embl.de/)进行编辑，以下即为我得到的一个进化树。

进化树.png

linux下phylip软件构建NJ树,SNP数据构建系统进化树相关推荐

linux下phylip软件构建NJ树,进化树构建-NJ法lpar;megarpar;
进化树构建-NJ法 V1.2 By ZHAO Yangguo, email: sunshine.zhao@yahoo.com.cn 通过长期的分析认为对于16S rDNA序列的进化(Phyligene ...
linux中主成分分析软件,基于全基因组snp数据进行主成分分析（PCA）
现将如何基于全基因组的SNP数据进行PCA分析流程记录下来: 1)全基因组snp数据格式为 .vcf 2)利用vcftools软件进行格式转换(Linux系统下:进入 /vcftools_0.1.13 ...
Linux下的软件安装方法介绍
本文目录一.Linux下的软件安装方式二.在线安装三.使用dpkg从本地磁盘安装deb软件包本文目录一.Linux下的软件安装方式通常Linux上的软件安装主要有三种方式: 在线安装从磁 ...
高压缩比压缩软件 linux,Linux下压缩软件对比
Linux下压缩软件对比发布时间:2007-02-15 00:44:03来源:红联作者:seais 最近在做系统备传输时发现压缩时间过常,于是简单做了一下Linux下的压缩软件的测试源文件3591 ...
第十章 Linux下RPM软件的安装与卸载
第十章 Linux下RPM软件的安装与卸载第一节 RPM软件包安装 rpm命名原则如: vsftpd-3.0.2-9.el7.x86_64.rpm vsftpd 3.0.2 9 e17 x86_6 ...
linux下常用软件下载
linux下常用软件下载: wget http://admin.ooopic.com/soft/linux/APC-3.1.2.tgz wget http://admin.ooopic.com/sof ...
automake linux,Linux下automake软件编译与发布快速入门
Linux下automake软件编译与发布快速入门 2008-04-22 eNet&Ciweek 进入编辑界面,输入内容如下: AUTOMAKE_OPTIONS=foreign bin_PRO ...
linux 树状结构图,linux下tree指令的用法, 树状图列出目录, 树状图逐级列出目录...
tree命令,主要功能是创建文件列表,将所有文件以树的形式列出来 linux下的tree就比较强大了,但一般系统并不自带这个命令,需要手动下载安装, 安装 sudo apt install tree ...
Linux完全清除防火墙策略,linux下的软件防火墙iptables——规则的查看与清除、定义默认策略,...
linux下的软件防火墙iptables--规则的查看与清除.定义默认策略软件防火墙IPTABLES --规则在Linux下查看和清除,定义默认策略,防火墙意味着用户限制某些ip或用户对其主机的访问. ...
Linux下的软件安装方式+源码安装软件cmatrix代码雨
Linux下的软件安装方式 1 rpm工具安装 2 yum工具安装 3 源码编译安装软件包类型 1 二进制软件包(=rpm软件包) 无需编译,直接安装根据计算机CPU架构类型和操作系统选择合适的软 ...

linux下phylip软件构建NJ树,SNP数据构建系统进化树

linux下phylip软件构建NJ树,SNP数据构建系统进化树相关推荐

最新文章

热门文章