Normalise read counts between two bam files
Entering edit mode
6.4 years ago
nanana ▴ 120

I'm trying to calculate the read depth ratio between two bam files in a specific genomic region.

One of these samples was sequenced deeper than the other - what's the recommend way to normalise read counts between the two?

Deeptools seems like a great standalone option for normalising between two bam files, can anyone outline the method they use to create a scaling factor?

For example, I am interested in calculating the read depth ratio in a particular window (window1 below) normalised for sequencing depth differences between two bam files, bam1 and bam2

Total number of mapped reads in bam1: 42973364
Reads in window1 in bam1: 9359
Total number of mapped reads in bam2: 101526426
Reads in window1 in bam2: 28347

How can I calculate normalised read depth ratio in this example?

genome next-gen sequencing • 3.3k views
Entering edit mode
6.4 years ago
ATpoint 84k

Have a look at deeptools bamCompare. It normalized two BAM files to each other, offering a number of normalization methods.

Entering edit mode

I just updated my question. I know that deep tools handles this, my question is, in the hypothetical example I give above, what is the appropriate way to calculate this by hand


Login before adding your answer.

Traffic: 3711 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6