Quantcast
Channel: Post Feed
Viewing all articles
Browse latest Browse all 3231

Samtools Mpileup And Overlapping Paired-End Reads

$
0
0

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.


Viewing all articles
Browse latest Browse all 3231

Trending Articles