multibamsummary vs multibigwigsummary

Question: multibamsummary vs multibigwigsummary

13 months ago by

skhil083 • 0 wrote:

I am in the process of comparing the similarity between two ChIP-seq conditions. According to the Deeptools documents, multibamsummary computes the read coverages for genomic regions for typically two or more BAM files whereas multibigwigsummary computes the average scores for each of the files in every genomic region.

I am having difficulty understanding the difference between the two tools and which one I should use. What is the difference between read coverage and average scores? Does the scale of the average score change depending on sequencing depth? Could someone shed some light on this please? Thank you in advance.

galaxy chip-seq • 675 views

ADD COMMENT • link •

modified 13 months ago by Devon Ryan • 1.9k • written 13 months ago by skhil083 • 0

13 months ago by

Devon Ryan • 1.9k

Germany

Devon Ryan • 1.9k wrote:

The scale of read coverage will vary with sequencing depth, so multiBamSummary is always affected by this. This is the benefit of making bigWig files and then using multiBigwigSummary, since you can use either bamCoverage or bamCompare (or bigwigCompare after bamCoverage) to first normalize sequencing depth or to input samples in some reasonable way before summarizing by region or bin. Consequently, I suggest that people use multiBigwigSummary instead of multiBamSummary.

ADD COMMENT • link written 13 months ago by Devon Ryan • 1.9k

My bigwig files were created with bedgraphtoBigwig following stretching with slopBed. Would that still be appropriate to use for multibigwigsummary? Since the scale of multibigwigsummary corresponds to average read scores calculated, would higher scores theoretically correspond to peaks with greater height or would it be width?

ADD REPLY • link written 13 months ago by skhil083 • 0

You can make the bigWig files from whatever you like. It's impossible to tell whether higher signal in a given bin corresponds to greater signal width vs. greater height. It's simply showing a greater average signal.

ADD REPLY • link written 12 months ago by Devon Ryan • 1.9k

So what is the biological translation of that greater average signal, if we are to assume the experiment was done well. Would it be greater transcription factor binding per se?

ADD REPLY • link written 12 months ago by skhil083 • 0

Yes, it's more of whatever you're IPing in that area.

ADD REPLY • link written 12 months ago by Devon Ryan • 1.9k

Similar posts • Search »