Hi, Is there any way to calculate base-specific read counts for a given genomic interval (including 1-base interval), for paired-end data aligned by Bowtie2 in BAM format? Thanks! Best, Yu Chuan