Hello all :)
I'm trying to check/correct my ChIP-Seq reads for GC% bias, and I got the following two plots for my Input and for my H3K27ac pulldown. I have excluded repetitive regions, but not mappability or peaks:
So based on the Input plot, I obviously have a little bias which I should correct for going forwards (unless there's a chance I should have called peaks for the input too?).
However on the plot for the ChIP pull down the result is wacky because I haven't excluded my peaks from the analysis. This makes sense since you can't correct for an overall GC bias if your chip is pulling down GC-rich regions specifically - and I understand what to do going forwards, but it seems like a lot of work. Call peaks, exclude regions, normalize the BAM, call peaks again. Is it therefore acceptable to just normalize the ChIP-Seq data to the bias of the Input (since both were sequenced on the same sequencer using the same kit manufacturer)?
All the best :)