Entering edit mode
5.4 years ago
ch8316f5eyu
▴
10
I wonder why the log2(IP/Input) value of the distal region is not about 0. For example:
The code used to generate the bamcompare file is:
bamCompare --scaleFactorsMethod readCount --minMappingQuality 1 --samFlagExclude 256 --operation log2 --outFileFormat bigwig --binSize 20 --pseudocount 5 -b1 Input.bw -o IP.bw
I think the possible reason is: There are too many reads in some unknown region in Input sample. Is this explanation reasonable? Or there are other reasons?