code: def merge(args):
      
      p = OptionParser(merge.__doc__)
      p.set_outfile(outfile="merged_results.delta")
      opts, args = p.parse_args(args)
  
      if len(args) < 3:
          sys.exit(not p.print_help())
  
      ref, query = args[:2]
      deltafiles = args[2:]
      outfile = opts.outfile
  
      ref = get_abs_path(ref)
      query = get_abs_path(query)
      fw = must_open(outfile, "w")
      print(" ".join((ref, query)), file=fw)
      print("NUCMER", file=fw)
      fw.close()
  
      for d in deltafiles:
          cmd = "awk 'NR > 2 {{print $0}}' {0}".format(d)
          sh(cmd, outfile=outfile, append=True)
docs: %prog merge ref.fasta query.fasta *.delta
  
      Merge delta files into a single delta.