code: def fill(args):
      
      p = OptionParser(fill.__doc__)
      p.add_option("--stretch", default=3, type="int",
                   help="MAX_STRETCH to pass to FillFragments [default: %default]")
      p.set_cpus()
  
      opts, args = p.parse_args(args)
  
      if len(args) != 1:
          sys.exit(not p.print_help())
  
      fastb, = args
      assert fastb == "frag_reads_corr.fastb"
  
      pcfile = "frag_reads_corr.k28.pc.info"
      nthreads = " NUM_THREADS={0}".format(opts.cpus)
      maxstretch = " MAX_STRETCH={0}".format(opts.stretch)
      if need_update(fastb, pcfile):
          cmd = "PathReads READS_IN=frag_reads_corr"
          cmd += nthreads
          sh(cmd)
  
      filledfastb = "filled_reads.fastb"
      if need_update(pcfile, filledfastb):
          cmd = "FillFragments PAIRS_OUT=frag_reads_corr_cpd"
          cmd += " PRECORRECT_LIBSTATS=True"
          cmd += maxstretch
          cmd += nthreads
          sh(cmd)
  
      filledfasta = "filled_reads.fasta"
      if need_update(filledfastb, filledfasta):
          cmd = "Fastb2Fasta IN=filled_reads.fastb OUT=filled_reads.fasta"
          sh(cmd)
docs: %prog fill frag_reads_corr.fastb
  
      Run FillFragments on `frag_reads_corr.fastb`.