code: def collect_read_counts(data, work_dir):
      
      out_file = os.path.join(work_dir, "%s-target-coverage.hdf5" % dd.get_sample_name(data))
      if not utils.file_exists(out_file):
          with file_transaction(data, out_file) as tx_out_file:
              params = ["-T", "CollectReadCounts", "-I", dd.get_align_bam(data),
                        "-L", tz.get_in(["regions", "bins", "target"], data),
                        "--interval-merging-rule", "OVERLAPPING_ONLY",
                        "-O", tx_out_file, "--format", "HDF5"]
              _run_with_memory_scaling(params, tx_out_file, data)
      return out_file
docs: Count reads in defined bins using CollectReadCounts.