I having trouble with analyzing my ChIP-seq data. I have an IP, input with replica (ChIP1, Input1, ChIP2, Input2). I merged the alignement files (.bam hg19) with bamtools merge, then called peaks. The problem is merging summs the reads at a given position. Like:
So, the final goal is, how to merge the samples as the read numbers averaged not summed at same positions?