I am using samtools mpileup to generate a pileup for paired-end RNA sequencing data. I am curious about how samtools handles pair mates whose read mappings overlap. With regard to the simple example below, for positions 7-11, are both pair mates enumerated?
position 1 6 11 16
reference ATGCATGCATGCATGC
pair mate 1 ATGCATGCATG
pair mate 2' GCATGCATGC (reverse complemented)
I am parsing the pileup output to quantify RNA editing at some positions and do not want to count an 'RNA molecule' twice just because the read pair mates 'overlapped.'
Thanks.