- Divide into equal-size (region_size/avg_bin_size) portions
- Emit the (chrom, start, end) coords of each portion
target_chroms = group_coords(ngfrills.parse_regions(target_bed, True))
if access_bed:
// Chromosome accessible sequence regions are given -- use them
access_chroms = group_coords(ngfrills.parse_regions(access_bed, True))
After Change
- Divide into equal-size (region_size/avg_bin_size) portions
- Emit the (chrom, start, end) coords of each portion
target_chroms = group_coords(RA.read(target_bed).coords())
if access_bed:
// Chromosome accessible sequence regions are given -- use them
access_chroms = group_coords(RA.read(access_bed).coords())