Calculate the enrichment of SVs in genomic regions

bpRegionEnrichment(
  ...,
  bp_data,
  dataType = "svBreaks",
  bed_file,
  bedDir = "/Users/Nick_curie/Desktop/misc_bed/features",
  chroms = c("2L", "2R", "3L", "3R", "4", "X", "Y"),
  restrict = TRUE,
  slop = 0,
  plot = TRUE,
  genome_length = 118274340,
  intersect = FALSE,
  outDir,
  parseName = FALSE,
  minHits = 10,
  multiple = TRUE
)