for file_name in os.listdir(target_directory):
file_name = os.path.join(target_directory, file_name)
with open(file_name, "r") as source:
bg = GRIMMReader.get_breakpoint_graph(source)
bg_graphs[file_name.split(".")[0]] = bg
bg = BreakpointGraph()
for br_gr in bg_graphs.values():
bg.update(br_gr, merge_edges=True)
target_multicolor = Multicolor("Anguilla_japonica")
print("Breakpoint graph stats:")
print(
"\t", "non-infinity nodes count:",
len(
list(node for node in bg.nodes()
if not BGVertex.is_infinity_vertex(node))))
normal_edges, infinity_edges = [], []
for edge in bg.edges():
if edge.is_infinity_edge:
infinity_edges.append(edge)
else:
normal_edges.append(edge)
print("\t", "non-infinity edges count:", len(normal_edges))
print("\t", "infinity edges count:", len(infinity_edges))
print("\t", "connected component count:",
len(list(bg.connected_components_subgraphs())))
print()
print()
