docs: def _nucmer_hits_to_percent_identity(nucmer_hits):
          '''Input is hits made by self._parse_nucmer_coords_file.
             Returns dictionary. key = contig name. Value = percent identity of hits to that contig'''
          percent_identities = {}
          max_lengths = {}
  
          for contig in nucmer_hits:
              max_length = -1
              percent_identity = 0
              for hit in nucmer_hits[contig]:
                  if hit.hit_length_qry > max_length:
                      max_length = hit.hit_length_qry
                      percent_identity = hit.percent_identity
              percent_identities[contig] = percent_identity
  
          return percent_identities