different number of reads mapped to plus strand and minus strand

[gfmgfm]

Hello all,

I have data of capture (sure select) followed by Illumina sequencing.
I see that I have more reads on the minus strand compared to the plus strand for example, for a certain position:
33 reads mapped to the plus strand and 96 reads mapped to the minus strand
88 reads mapped to the plus strand and 115 reads mapped to the minus strand

Does anyone has explanations for that?

[seb567]

Hello all,

I have data of capture (sure select) followed by Illumina sequencing.
I see that I have more reads on the minus strand compared to the plus strand for example, for a certain position:
33 reads mapped to the plus strand and 96 reads mapped to the minus strand
88 reads mapped to the plus strand and 115 reads mapped to the minus strand

Does anyone has explanations for that?

DNA fragmentation is random, and so is cluster generation.

Did you look at the entire set of mapped reads or just a few ?

[gfmgfm]

Hi,

Thanks for the reply.

I looked at 2 parameters:
1. The number of reads mapped (on targets) - on plus strand and on minus strand:
For some samples it is the same, for many samples it is very different (~x2).
2. For each target base - I looked at a histogram of the number of reads mapped to that base on the plus strand divided by the total number of reads mapped to that base. It seems that there is some bias - and that there are bases that are covered only or mainly by one of the strands.