这个代码没有测试过,不一定能运行,数据量不大用mega也行的;
使用的脚本是来自于课程中的NJ_tree_phylip.py文件
整个命令为:
python ~/NJ_tree_phylip.py -i ~/gene_family_demo/03.gene_tree_analysis/pep_aligned_trimed.fa -t prot -o ~/gene_family_demo/03.gene_tree_analysis -n NJ
报错为:os.rename('/'+args.outdir.strip()+'/outfile','/'+args.outdir.strip()+'/infile')
FileNotFoundError: [Errno 2] No such file or directory: '//work/gene_family_demo/03.gene_tree_analysis/outfile' -> '//work/gene_family_demo/03.gene_tree_analysis/infile'