code: def zip_sequences(G, allseqs, color="white"):
      
      for s in zip(*allseqs):
          groups = defaultdict(list)
          for x in s:
              part = x.split('_', 1)[1]
              groups[part].append(x)
          for part, g in groups.items():
              with G.subgraph(name="cluster_" + part) as c:
                  for x in g:
                      c.node(x)
                  c.attr(style="invis")
docs: Fuse certain nodes together, if they contain same data except for the
      sequence name.