code: def regions_coverage(bed_file, target_name, data):
      
      ready_bed = tz.get_in(["depth", target_name, "regions"], data)
      if ready_bed:
          return ready_bed
      else:
          return run_mosdepth(data, target_name, bed_file).regions
docs: Generate coverage over regions of interest using mosdepth.