python 进化树_Python / ete3:将最密切相关的叶子定位到系统树中的特定物种

我正在使用Python包ete3。我有以下树木:

((Species1_order1,(Species2_order2,Species3_order2)),Species4_order3,Species5_order5);

我想看到与树中特定节点(此处为树Species1_order1)关系最密切的叶。在此示例中,最紧密相关的叶子是Species2_order2/Species3_order2和和Species4_order3/Species5_order5。

码:

tree=ete3.Tree('((Species1_order1, \(Species2_order2, Species3_order2)), \Species4_order3, Species5_order5);')

新的例子:

tree=ete3.Tree('((((((A,B),C),D),(E,F)),G),(H,I));')

我得到的结果是:

A B C D E F G H IA0.02.03.04.06.06.06.08.08.0B2.00.03.04.06.06.06.08.08.0C3.03.00.03.05.05.05.07.07.0D4.04.03.00.04.04.04.06.06.0E6.06.05.04.00.02.04.06.06.0F6.06.05.04.02.00.04.06.06.0G6.06.05.04.04.04.00.04.04.0H8.08.07.06.06.06.04.00.02.0I8.08.07.06.06.06.04.02.00.0

但是,例如,E和F在树中与A,B,C和D的距离相等,因此,它们似乎比D更衣。

一个好的矩阵结果应该是:

A B C D E F G H IA012

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
你可以使用一些生物信息学工具来比较物种和基因之间的形拓扑结构差异。在这里,我可以为你提供一个Python代码示例,使用Biopython库和ETE Toolkit来完成这个任务。 首先,你需要安装BiopythonETE Toolkit库。你可以使用以下命令来安装它们: ```python pip install biopython pip install ete3 ``` 然后,你可以使用以下代码来比较物种和基因形拓扑结构: ```python from Bio import Phylo from ete3 import Tree species_tree_file = "/home/dengwei/NTF_data/rotted_species_tree/species_tree.treefile" gene_tree_folder = "/home/dengwei/NTF_data/rotted_gene_tree" compare_tree_folder = "/home/dengwei/NTF_data/compare_tree" # 读取物种 species_tree = Tree(species_tree_file) # 获取基因文件列表 gene_tree_files = os.listdir(gene_tree_folder) # 遍历基因文件 inconsistent_trees = 0 for gene_tree_file in gene_tree_files: gene_tree_path = os.path.join(gene_tree_folder, gene_tree_file) # 读取基因 gene_tree = Tree(gene_tree_path) # 比较形拓扑结构 if species_tree.compare(gene_tree, unrooted=True) != 0: inconsistent_trees += 1 # 将不一致的基因保存到compare_tree_folder compare_tree_path = os.path.join(compare_tree_folder, gene_tree_file) gene_tree.write(outfile=compare_tree_path, format=5) print("共有{}棵基因物种不一致。".format(inconsistent_trees)) ``` 上述代码将遍历基因文件夹中的每个基因文件,与物种进行形拓扑结构比较。如果两者不一致,则将不一致的基因保存到compare_tree_folder中,并统计不一致的基因数量。 请确保你的文件路径正确,并且你拥有读取和写入文件的权限。希望这能帮到你!

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值