docs: def _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.
  
      :param feature: pybedtools.Interval object
      :param stranded: If `stranded=True`, then only counts signal on the same
          strand as `feature`.
      """
      if isinstance(feature, basestring):
          feature = helpers.tointerval(feature)
      if stranded:
          strand = feature.strand
      else:
          strand = '.'
  
      count = 0
      for al in reader[feature]:
          if stranded and al.strand != strand:
              continue
          count += 1
      return count