count DNA sequence reads in BAM files