Python:如何将系统发育树转换为圆形格式?

时间:2021-05-25 11:15:50

标签: python biopython

我有一个比对序列的 fasta 文件,我想为其创建系统发育树。我可以像这样创建一个普通的分支树...... Normal branched phylogenetic tree

但想将其转换为如下所示...Circular phylogenetic tree

到目前为止我使用的代码如下(out_file 是对齐的 fasta 文件的文件路径)


from Bio import AlignIO
from Bio import Phylo
from Bio.Phylo.TreeConstruction import DistanceTreeConstructor
import matplotlib
import matplotlib.pyplot as plt

align = AlignIO.read(out_file, "fasta")
# print(align)

from Bio.Phylo.TreeConstruction import DistanceCalculator
calculator = DistanceCalculator('identity')
dm = calculator.get_distance(align)


constructor = DistanceTreeConstructor()
tree = constructor.upgma(dm)


Phylo.draw_ascii(tree)

0 个答案:

没有答案
相关问题