code: def parse_delta(filename):
      
      aln_length, sim_errors = 0, 0
      for line in [l.strip().split() for l in open(filename, "r").readlines()]:
          if line[0] == "NUCMER" or line[0].startswith(">"):  # Skip headers
              continue
          # We only process lines with seven columns:
          if len(line) == 7:
              aln_length += abs(int(line[1]) - int(line[0]))
              sim_errors += int(line[4])
      return aln_length, sim_errors
docs: Returns (alignment length, similarity errors) tuple from passed .delta.
  
      - filename - path to the input .delta file
  
      Extracts the aligned length and number of similarity errors for each
      aligned uniquely-matched region, and returns the cumulative total for
      each as a tuple.