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?
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