diff --git a/compara/reconstruct.py b/compara/reconstruct.py
index <HASH>..<HASH> 100644
--- a/compara/reconstruct.py
+++ b/compara/reconstruct.py
@@ -76,14 +76,17 @@ def fuse(args):
     logging.debug("Total families: {}, Gene members: {}"\
                 .format(len(families), len(allowed)))
 
+    # TODO: Use C++ implementation of BiGraph() when available
+    # For now just serialize this to the disk
     G = BiGraph()
     for bedfile in bedfiles:
         bed = Bed(bedfile, include=allowed)
         add_bed_to_graph(G, bed, families)
 
-    for path in G.iter_paths():
-        m, oo = G.path(path)
-        print m
+    G.write(filename="graph.edges")
+    #for path in G.iter_paths():
+    #    m, oo = G.path(path)
+    #    print m
 
 
 def adjgraph(args):
[compara] Write edges in reconstruct.fuse()