Question

我正在尝试编辑phyloxml文件中的提示标签，因此它们只包含4位数字。例如，我目前在phyloxml文件中名称下的是左侧的ID。我想要的是右边的ID：

ACOM042150-PA ACOM
AQUA008971-PA AQUA
AGAP002137-PA AGAP
AARA006802-PA AARA

到目前为止我的代码可以使用替换方法来挑选我想要的字符：

tree = Phylo.read("GSCGT000003.xml", "phyloxml")
for i, clade in enumerate(tree.find_clades(name=True)):
    print clade.name.replace(clade.name, clade.name[0:4])

这会将所需的名称打印到终端，但我不确定如何在文件中替换，或者使用更改写入新文件。我可以使用Phylo.write写入新文件，但我似乎无法使用更改来编写文件。有什么想法吗？谢谢。

Answer 1

您必须更改每个分支的名称，然后将树保存到新文件：

tree = Phylo.read("GSCGT000003.xml", "phyloxml")

for clade in tree.find_clades(name=True):
    # This changes the clade name to only the first four chars
    clade.name = clade.name[:4]

# At this point all clade names has been changed in-place
Phylo.write(tree, "GSCGT000003_modified.xml", "phyloxml")

Phyloxml：在文件中查找和替换

1 个答案: