count reads with a combination of specific SNPs in bam

Question

[genya35]

Hello,

I have a small bam file and would like to count reads that have specific bases in two different positions(A and G), 25 bp apart. Which tool can I use?

Thanks