This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Normalizing BAM Files

Hi!

I am back again with one more naive question.

What according to you is the best way to normalize multiple BAM files to one another?

My goal is to simply normalize a set of BAM files, then plot the avg read coverage.

I've attached an example of a read metagene plot showing the avg read coverage with BAM files that are not normalized to one another. To reiterate, my goal is to normalize these BAM files, then make the same type of plot.

Read Metagene Plot

bam
• 3,312 views • link
written by kstangline • kstangline

Can you please define "normalize"? Do you mean to tweak the total number of aligned reads so all BAMs have comparable numbers of aligned reads?

• 0 views • link
written by Ram Ram

So the plot above is average read coverage, but I'm needing to make a new plot (the same type of figure as the one above) showing averaged normalized read coverage, which is normalized to the depth of the library.

• 0 views • link
written by kstangline • kstangline

1 answer

You have a few options:

  1. 1x normalize the BAM files with something like bamCoverage
  2. Use multiBamSummary with the --scalingFactors to compute appropriate scaling factors in a more robust way and then use those to scale the per-BAM signals.
• 0 views • link
written by Devon Ryan Devon Ryan

Log in to answer this question.

More posts like this