metaseq.array_helpers._local_count

metaseq.array_helpers._local_count(reader, feature, stranded=False)

The count of genomic signal (typcially BED features) found within an interval.

Usually this only makes sense for BED or BAM (not bigWig) files.

Parameters:

• feature – pybedtools.Interval object
• stranded – If stranded=True, then only counts signal on the same strand as feature.