def index_mutations(con):
"""Builds an index of all mutations"""
cur = con.cursor()
for msaid in get_alignment_method_ids(con):
for modelid in get_phylo_modelids(con):
newick = get_anc_cladogram(con, msaid, modelid)
t = Tree()
t.read_from_string(newick, "newick")
for edge in t.preorder_edge_iter():
if edge.head_node == None or edge.tail_node == None:
continue
if edge.head_node.label == None or edge.tail_node.label == None:
continue
print msaid, modelid, edge.head_node.label, edge.tail_node.label
anc1name = "Node" + edge.head_node.label.__str__()
anc2name = "Node" + edge.tail_node.label.__str__()
index_mutations_helper(con, msaid, modelid, anc1name, anc2name)
