服务器建树

第一步

创建环境,接下来的软件均安装在此处
conda create -n tree conda activate tree
安装多序列比对软件
conda install -c bioconda mafft
安装序列修剪软件
conda install -c bioconda trimal
安装建树软件
conda install -c bioconda fasttree

多序列比对

mafft --auto input>output

修剪

trimal -automated1 -phylip -in 1.fasta -out 1.trimal.fasta

建树

FastTree从核苷酸或蛋白质序列的排列推断出近似最大似然系统发生树,可以直接拖到mega进行可视化
fasttree 1.trimal.fasta>1.Newick

或者用muscle比对
/public1/home/stu_zhangyingyin/project/orthofinder/muscle3.8.31_i86linux64 -in 1.fasta -out 1.afa -maxiters 1 -diags

MUSCLE还可以建树
UPGMA tree:
muscle -maketree -in seqs.afa -out seqs.phy
Neighbor-Joining tree :
muscle -maketree -in seqs.afa -out seqs.phy -cluster neighborjoining

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容