示例#1
0
    #name_path = sys.argv[2]

    #read nodes
    (nodes, taxes_ids) = read_nodes()

    #get nodes by input rank
    filtered_nodes_dic = filter_by_rank(nodes, taxes_ids, rank)

    filtered_nodes_dic = nodes
    #group id
    #nodes_dic[group_id].rank

    #for each node find its lineage
    find_lineages(filtered_nodes_dic, taxes_ids)

    lineages = []
    for id, node in nodes.items():
        lineages.append(node.lineage)
        #print (str(node))

    tree = csv_to_weightless_newick(lineages)

    t = Tree(tree + ";", format=8)
    print(t)

    for node in t.traverse("postorder"):

        print(node.name)

    print(t.get_node_by_name("12"))