code: def _population_load_script(work_bams, names, chrom, pairmode, items):
      
      bed_file = _get_regional_bed_file(items[0])
      if bed_file:
          return _population_prep_targeted.format(bam_file_str=",".join(work_bams), names_str=",".join(names),
                                                  chrom=chrom, num_cores=0, pairmode=pairmode, bed_file=bed_file)
      else:
          return _population_prep.format(bam_file_str=",".join(work_bams), names_str=",".join(names),
                                         chrom=chrom, num_cores=0, pairmode=pairmode)
docs: Prepare BAMs for assessing CNVs in a population.