Question: bamCompare the result if 0
gravatar for HI
7 days ago by
HI0 wrote:

I am using bamCompare to normalize my Chip-seq data (antibody + input). The reference genome is bigger than 10G, and every chromosome size is bigger than 600M, I run the bamCompare for each chromosome seperately(because there is no big enough memory to do with the whole genome; SES normalization; binsize=50), when I check the results of bamCompare, I found that the log2 ratio in every chromosome site bigger than 538M is 0 (even if I set --region chr1:500000000:830000000), I didn't find any information on the official website, then what should I do with this? Or is there any other software to calculate log2 ratio? Thanks!

ADD COMMENTlink modified 7 days ago by Devon Ryan88k • written 7 days ago by HI0

What is the ultimate goal that you need these values for? There are for sure alternative approaches that will work. Did you try using the enrichment or significance scores from a peak caller such as MACS2 for your analysis?

ADD REPLYlink modified 7 days ago • written 7 days ago by ATpoint14k
gravatar for Devon Ryan
7 days ago by
Devon Ryan88k
Freiburg, Germany
Devon Ryan88k wrote:

I recently replied to an identical question on the deepTools email list. To my knowledge bigWig files don't support chromosomes larger than ~538 megabases, so there's nothing deepTools can do about that.

ADD COMMENTlink written 7 days ago by Devon Ryan88k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2204 users visited in the last hour