code: def _finalize_merge(out_file, bam_files, config):
      
      # Ensure timestamps are up to date on output file and index
      # Works around issues on systems with inconsistent times
      for ext in ["", ".bai"]:
          if os.path.exists(out_file + ext):
              subprocess.check_call(["touch", out_file + ext])
      for b in bam_files:
          utils.save_diskspace(b, "BAM merged to %s" % out_file, config)
docs: Handle indexes and cleanups of merged BAM and input files.