code: def _get_regional_bed_file(data):
      
      variant_regions = bedutils.merge_overlaps(tz.get_in(["config", "algorithm", "variant_regions"], data),
                                                data)
      is_genome = data["config"]["algorithm"].get("coverage_interval", "exome").lower() in ["genome"]
      if variant_regions and utils.file_exists(variant_regions) and not is_genome:
          return variant_regions
docs: If we are running a non-genome analysis, pull the regional file for analysis.