雖然簡單陨囊,但絕對首發(fā)……
該給自己的系統(tǒng)發(fā)育構(gòu)建流程升級了爷怀,新版本軟件不僅速度快,而且對用戶友好匀伏!
#mafft比對
linsi allnew.fasta > allnew.fasta.mafft
#Gblocks修整
Gblocks allnew.fasta.mafft -t=d -e=.gb
#修改序列名
sed -i '/^>/s/ /_/g' allnew.fasta.mafft.gb
#模型檢驗(yàn)
wget https://github.com/ddarriba/modeltest/files/3790700/modeltest-ng-0.1.6-static-linux64.tar.gz
unzip modeltest-ng-0.1.6-static-linux64.tar.gz
modeltest-ng-static -i allnew.fasta.mafft.gb -d nt
#構(gòu)建ML樹
wget https://github.com/amkozlov/raxml-ng/releases/download/0.9.0/raxml-ng_v0.9.0_linux_x86_64.zip
unzip raxml-ng_v0.9.0_linux_x86_64.zip
mv data/ INSTALL.txt LICENSE.txt raxml-ng raxml-ng_v0.9.0
/home/zhangjj/software/raxml-ng_v0.9.0/raxml-ng --all --msa allnew.fasta.mafft.gb --model TVM+G4 --force
#allnew.fasta.mafft.gb.raxml.support即為最終結(jié)果
modeltest-ng 的速度是真的飛快洒忧!
modeltest-ng 和 raxml-ng簡直是絕配!
Reference
Diego D , David P , Kozlov A M , et al. ModelTest-NG: A New and Scalable Tool for the Selection of DNA and Protein Evolutionary Models[J]. Molecular Biology and Evolution.
RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference