Question: Filtering high quality reads from bowtie2
gravatar for yulia.gordon
4.3 years ago by
European Union
yulia.gordon0 wrote:

Hi everyone! I've encountered an interesting discrepancy when trying to filter the results of my paired end CHIP-SEQ. When using bowtie2 in Galaxy interface it is shown how many mapped, uniquely mapped and more than once mapped reads were found. When I filter the Bowtie2 results, leaving only reads that were mapped in a proper pair and with MAPQ >10 ( I understood MAPQ>10 is supposed to be a good enough filter to use in this case) I stay with more reads than uniquely plus non-uniquely mapped reads, which doesn’t make any sense. For example: The number of all reads is 3,864,595, The number of all mapped reads (uniquely+ more than once) is 2,156,879 The number of properly paired, MAPQ>10 reads is 2,374,371

Does anybody have any explanation for this? Thank you! Yulia

chip-seq • 1.4k views
ADD COMMENTlink modified 3.0 years ago by Biostar ♦♦ 20 • written 4.3 years ago by yulia.gordon0

That shouldn't actually be possible. Are you sure that you didn't either merge two BAM files or that you're actually filtering for MAPQ>10?

ADD REPLYlink written 4.3 years ago by Devon Ryan92k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1299 users visited in the last hour