码迷,mamicode.com
首页 > 其他好文 > 详细

利用phylip软件对SNP数据构建进化树

时间:2020-05-31 23:20:26      阅读:292      评论:0      收藏:0      [点我收藏+]

标签:下载   lease   命令   boot   install   make   bsp   tps   结果   

1、下载、安装 phylip软件

    官网: http://evolution.genetics.washington.edu/phylip.html 

wget http://evolution.gs.washington.edu/phylip/download/phylip-3.697.tar.gz
tar -xzvf phylip-3.697.tar.gz
cd phylip-3.697/
cd src/
make -f Makefile.unx install

技术图片

 

 

 安装成功的标志。

 可执行程序在exe路径下:

 技术图片

 

 技术图片

 

 2、下载测试数据 ,提取一小部分进行测试:

链接:https://pan.baidu.com/s/1VxgJK5asCEjukOeA1V5rlg 
提取码:e69n
head -n 5 test.ped | cat - <(tail -n 5 test.ped) > a && mv a test.ped

技术图片

 

 3、参考 下面文章进行 构建进化树  https://zhuanlan.zhihu.com/p/85978856   

   样本ID必须是10个字符,这点很重要!!!使用如下命令修改:

 

##把间隔符改为tab
plink --file test --recode tabx --sheep --out test;rm -f *.log *.nosex 

##把个体ID调整为10个字符
awk {print $2 "xxx000000000"} test.ped |cut -b 1-10|paste - test.ped |awk {OFS = "\t"}{print $2,$1,$0}| cut -f 3-5 --complement > a && mv a test.ped

##转化为vcf格式:
plink --file test --sheep --recode vcf-iid --out test;rm -f *.log *.nosex

4、下载 vcf2phylip.py脚本,用于将vcf文件转化为phylip软件需要的格式:

     下载地址:https://github.com/edgardomortiz/vcf2phylip/releases 

wget https://github.com/edgardomortiz/vcf2phylip/archive/v2.0.zip
unzip v2.0.zip
cd vcf2phylip-2.0/

  技术图片

 

 5、使用  vcf2phylip.py 将test.vcf转化为phylip软件需要的格式:

python /home/liujiaxin/software/vcf2phylip-2.0/vcf2phylip.py -i  test.vcf

 技术图片

 

 6、准备配置文件 ,参照 https://zhuanlan.zhihu.com/p/85978856

a、seqboot.par ,   test.min4.phy是上一步vcf转化格式生成的文件。 注意  test.min4.phy 后面不要加注释,否则报错!!!!

技术图片

 

 b、dnadist.par   注意:第一行不要加注释!!!!!

技术图片

 

 c、neighbor.par  注意:第一行不要加注释!!!!!

技术图片

 

 d、consense.par    注意:第一行不要加注释!!!!!

 技术图片

 

 看下 当前目录:

技术图片

 

 7、运行以下程序:

/home/liujiaxin/software/phylip-3.697/exe/seqboot < ./seqboot.par && mv ./outfile ./seqboot.out
/home/liujiaxin/software/phylip-3.697/exe/dnadist < ./dnadist.par &&  mv ./outfile ./dnadist.out
/home/liujiaxin/software/phylip-3.697/exe/neighbor < ./neighbor.par && mv  ./outfile ./nei.out && mv ./outtree ./nei.tree
/home/liujiaxin/software/phylip-3.697/exe/consense < ./consense.par && mv ./outfile ./cons.out && mv ./outtree ./constree

8、看下当前目录,constree为构建进化树的文件

技术图片

 

 9、利用 itol绘图、美化,官网:https://itol.embl.de/ 

  技术图片

 

  技术图片

 

 10、绘图 结果:

      技术图片

 

 

     

 

利用phylip软件对SNP数据构建进化树

标签:下载   lease   命令   boot   install   make   bsp   tps   结果   

原文地址:https://www.cnblogs.com/liujiaxin2018/p/13022301.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!