构建进化树可以用DNA序列,也可以用protein序列。本文就如何利用蛋白序列进行进化树的构建步骤进行详细解析:

一、利用phylip进行蛋白序列进化树分析步骤

1.       运行Clustalw,得到protein.phy文件;

2.       运行seqboot,输入protein.phy文件,输出protein_1.outfile。R选择1000,seed number: 3,然后选择Y。将得到的protein_1.outfile改名为protein_2.infile。

3.       运行protpars,输入protein_2.infile,输出protein_2.outfile和protein_2.outtree。J选择3(seed number),5(jumple number),M值选择D,然后输入1000 (sets)。最后选择Y。将得到的protein_2.outtree改为.

4.       运行consense,输入protein_3.intree,选择Y。得到protein_3.outfile和protein_3.outtree

到此为止,利用seqboot计算的bootstrap值已经完成,但是protein_3.outtree构建的tree没有branch length数据,无法直观的显示序列之间的差异。此时,应该利用tree-puzzle来计算branch length。http://java.com/en/download/index.jsp )

Tree-puzzle使用说明

利用tree-puzzle计算branch length步骤:

1. Windows下,进入tree-puzzle所在src文件夹,双击Puzzle,弹出对话框,输入文件。文件为clustalW运行的protein.phy文件。

2. 然后依次选择K、W、C,输入4,最后输入Y确定。

3. 输入consense操作得到的文件protein_3.outtree。将要生成的文件改名为protein.complz。

4. 点击enter运行。

5. 将得到的consense文件protein_3.outtree和protein.complz文件同时拷出到新文件夹。

6. 利用eATV进行bootstrap和branchlength的显示,输出PDF文件。

7. 将相应的bootstrap值,算成百分比,通过Adobe illustrator 加到puzzle计算的树上。

8. Done

eATV是一个java应用程序,需要先安装java运行环境程序(http://java.com/en/download/index.jsp )

利用DNA序列建进化树的步骤与上面大同小异,参照phylip documentation文件。(http://evolution.genetics.washington.edu/phylip/phylip.html )

这些步骤均可在windows下完成。但是如果序列较大,建议在linux系统下操作。其中protpars步骤特别耗时,经常需要过夜运行,请一定要耐心等待。祝大家顺利!

Phylip 下载:http://evolution.genetics.washington.edu/phylip/getme.html

Tree-puzzle下载:http://www.tree-puzzle.de/#download

eATV下载:http://www.softpedia.com/get/Science-CAD/ATV.shtml (rather small: 93kb)

上述方法及其参数设置得到了密苏里大学哥伦比亚分校张学成博士的悉心指导,在此表示衷心的感谢!(http://www.staceylab.missouri.edu/personnel/1232960286/ )

linux下phylip软件构建NJ树,利用phylip构建进化树详解相关推荐

  1. linux如何确定共享库路径,摘录Linux下动态共享库加载时的搜索路径详解

    对动态库的实际应用还不太熟悉的读者可能曾经遇到过类似"error while loading shared libraries"这样的错误,这是典型的因为需要的动态库不在动态链接器 ...

  2. Linux进程最大socket数,Linux下高并发socket最大连接数所受的各种限制(详解)

    1.修改用户进程可打开文件数限制 在Linux平台上,无论编写客户端程序还是服务端程序,在进行高并发TCP连接处理时,最高的并发数量都要受到系统对用户单一进程同时可打开文件数量的限制(这是因为系统为每 ...

  3. linux修改zip中文件,Linux下修改jar、zip等文件内的文件详解

    Linux下修改jar.zip等文件内的文件详解 Linux下修改jar.zip等文件内的文件详解 看到很多同事在Linux环境下修改jar包内的文本文件或zip中的文本文件时,经常是先把jar包或z ...

  4. Linux下获取详细硬件信息的工具:Dmidecode命令详解

    Dmidecode 这款软件允许你在 Linux 系统下获取有关硬件方面的信息.Dmidecode 遵循 SMBIOS/DMI 标准,其输出的信息包括 BIOS.系统.主板.处理器.内存.缓存等等.偶 ...

  5. 汇总 Linux下获取详细硬件信息的工具:Dmidecode命令详解

    Dmidecode 这款软件允许你在 Linux 系统下获取有关硬件方面的信息.Dmidecode 遵循 SMBIOS/DMI 标准,其输出的信息包括 BIOS.系统.主板.处理器.内存.缓存等等.偶 ...

  6. Linux下ps -ef 和 ps aux 的区别及格式详解

    Linux下显示系统进程的命令ps,最常用的有ps -ef 和ps aux.这两个到底有什么区别呢? 这里要说一下 BSD vs System V 这里需要说明的是:Linux不能称为"标准 ...

  7. linux下的挂载点和分区是什么关系(详解挂载点)

    百度定义: 挂载点实际上就是linux中的磁盘文件系统的入口目录,类似于windows中的用来访问不同分区的C:.D:.E:等盘符.其实winxp也支持将一个磁盘分区挂在一个文件夹下面,只是我们C:. ...

  8. Linux下ps -ef和ps aux的区别及格式详解

    Linux下显示系统进程的命令ps,最常用的有ps -ef 和ps aux.这两个到底有什么区别呢?两者没太大差别,讨论这个问题,要追溯到Unix系统中的两种风格,System V风格和BSD 风格, ...

  9. Linux下Gcc生成和使用静态库和动态库详解

    参考文章:http://blog.chinaunix.net/uid-23592843-id-223539.html 一.基本概念 1.1什么是库 在windows平台和linux平台下都大量存在着库 ...

  10. Linux下Gcc生成和使用静态库和动态库详解(转)

    一.基本概念 1.1什么是库 在windows平台和linux平台下都大量存在着库. 本质上来说库是一种可执行代码的二进制形式,可以被操作系统载入内存执行. 由于windows和linux的平台不同( ...

最新文章

  1. 计算机科学和Python编程导论(一) 计算机相关理论
  2. 解决webpack打包vue项目后,部署完成后,刷新页面页面404
  3. Linux基础-yum软件包管理
  4. 【软件安装】sublime安装方法
  5. linux加密格式化吗,linux环境下给文件加密/解密的方法
  6. 导出Excel java
  7. C++中各种智能指针的实现及弊端(三)
  8. Pandas索引操作及高级索引——reindex()方法
  9. linux 安装swoole
  10. 留言送书丨《Oracle高性能系统实战大全》包邮到家
  11. jquery给title标签赋值
  12. System Center Operations Manager 简介 [SCOM中文系列之一]
  13. 近世代数——Part2 群:基础与子群
  14. laravel文件上传与下载
  15. matlab 傅里叶变换时间窗 频率窗,短时傅里叶变换窗函数窗口宽度的选择
  16. 刘乾四川大学计算机学院刘乾,刘乾亮
  17. 前鹅厂面试官亲手整理出的高频必问软件测试基础面试题(1)!!!
  18. matlab半小提琴图,不会编程,也可以画小提琴图啦!
  19. 【晶振不起振原因分析】
  20. CDH框架搭建【小白易上手】

热门文章

  1. 百度网盘下载一直请求中问题解决
  2. Matlab:成功解决In an assignment A(I)=B,the number of elements in B and I must be the same
  3. 微信公众号、订阅号、服务号、企业号的区别
  4. [统计学理论基础] 统计方法—F检验
  5. c语言 中的求和函数,0基础理解C语言调用汇编实现求和函数
  6. mvc:annotation-driven
  7. K3 CLOUD返工生产成本方案——循环计算
  8. sed替换指定字符串为某变量的值
  9. 装系统:主分区、扩展分区、逻辑分区,引导(启动)分区、系统分区、活动分区
  10. SQLHelper蓦然回首