我有一个比对序列的 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)