Normalise read counts between two bam files
1
3
Entering edit mode
6.2 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.2k views
ADD COMMENT
2
Entering edit mode
6.2 years ago
ATpoint 82k

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

ADD COMMENT
0
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

ADD REPLY

Login before adding your answer.

Traffic: 2572 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6