code: def _nucmer_hits_to_assembly_coords(nucmer_hits):
          '''Input is hits made by self._parse_nucmer_coords_file.
             Returns dictionary. key = contig name. Value = list of coords that match
             to the reference gene'''
          coords = {}
          for l in nucmer_hits.values():
              for hit in l:
                  if hit.qry_name not in coords:
                      coords[hit.qry_name] = []
                  coords[hit.qry_name].append(hit.qry_coords())
  
          for scaff in coords:
              pyfastaq.intervals.merge_overlapping_in_list(coords[scaff])
  
          return coordf _nucmer_hits_to_assembly_coords(nucmer_hits):
          '''Input is hits made by self._parse_nucmer_coords_file.
             Returns dictionary. key = contig name. Value = list of coords that match
             to the reference gene'''
          coords = {}
          for l in nucmer_hits.values():
              for hit in l:
                  if hit.qry_name not in coords:
                      coords[hit.qry_name] = []
                  coords[hit.qry_name].append(hit.qry_coords())
  
          for scaff in coords:
              pyfastaq.intervals.merge_overlapping_in_list(coords[scaff])
  
          return coords
docs: Input is hits made by self._parse_nucmer_coords_file.
             Returns dictionary. key = contig name. Value = list of coords that match
             to the reference gene