code: def batchlobstr(args):
      
      p = OptionParser(batchlobstr.__doc__)
      p.add_option("--haploid", default="chrY,chrM",
                   help="Use haploid model for these chromosomes")
      p.set_cpus()
      opts, args = p.parse_args(args)
  
      if len(args) != 1:
          sys.exit(not p.print_help())
  
      bamlist, = args
      cmd = "python -m jcvi.variation.str lobstr TREDs"
      cmd += " --input_bam_path {}"
      cmd += " --haploid {}".format(opts.haploid)
      cmd += " --simulation"
      cmds = [cmd.format(x.strip()) for x in open(bamlist).readlines()]
      p = Parallel(cmds, cpus=opts.cpus)
      p.run()
docs: %prog batchlobstr bamlist
  
      Run lobSTR on a list of BAMs. The corresponding batch command for TREDPARSE:
      $ tred.py bamlist --haploid chr4 --workdir tredparse_results