群体遗传进化进化树分析方法:

群体遗传进化进化树分析方法:

cd $workdir  #回到工作目录
mkdir 01.phylo_tree
cd 01.phylo_tree
#文件格式转换
run_pipeline.pl  -Xmx5G -importGuess  $workdir/00.filter/clean.vcf.gz  \
    -ExportPlugin -saveAs supergene.phy -format Phylip_Inter
#最大似然法构建进化树
#方法1:fasttree 构建进化树
fasttree -nt -gtr  supergene.phy   >  fasttree.nwk
#方法2:iqtree 构建进化树  设置boots值,最大似然法
iqtree2 -s supergene.phy -st DNA -T 2  -mem 8G \
    -m  GTR  -redo \
    -B 1000 -bnni \
    --prefix iqtree 
    
#方法3:raxml 构建进化树 最大似然法
#raxml-ng  -msa supergene.phy --model GTR  --prefix raxml_tree \
    #    --threads 2 --seed 123   1>raxml.log 2>raxml.err
## with bootstrap
#raxml-ng -all  -msa supergene.phy --model GTR --bs-trees 1000 \
    #    --prefix raxml_tree_bootstrap --threads 2 --seed 123   1>raxml_bs.log 2>raxml_bs.err
#方法4:phylip NJ法构建进化树  Phylip 格式对样品ID字符要求不超过10个,不然会截断
cp supergene.phy infile
echo -e "Y\n" |dnadist
cp infile.dist infile
echo -e "Y\n"  |neighbor
mv outtree outtree.nwk

#进化树美化:https://www.omicsclass.com/article/671



视频课程推荐:
  • 发表于 2021-09-01 13:25
  • 阅读 ( 2327 )
  • 分类:遗传进化

0 条评论

请先 登录 后评论
omicsgene
omicsgene

生物信息

652 篇文章

作家榜 »

  1. omicsgene 652 文章
  2. 安生水 321 文章
  3. Daitoue 167 文章
  4. 生物女学霸 120 文章
  5. 红橙子 78 文章
  6. CORNERSTONE 72 文章
  7. rzx 63 文章
  8. xun 61 文章