This script merges the loop locations of different different resolutions.
Loops need to have format as follows:
chr start end chr start end
A merge happens if x and y position of a loop overlap with x and y position of another loop; all loops are considered as an overlap within +/- the bin size of the lowest resolution. I.e. for a loop with coordinates x and y, the overlap to all other loops is search for (x - lowest resolution) and (y + lowest resolution). If two or more locations should be merged, the one with the lowest resolution is taken as the merged loop.
$ hicMergeLoops -i gm12878_10kb.bedgraph gm12878_5kb.bedgraph gm12878_25kb.bedgraph -o merged_result.bedgraph -r 25000
Please recall: We work with binned data i.e. the lowest resolution is therefore the one where we merge the most bases into one bin. In the above example the lowest resultion is 25kb, the highest resolution is 5kb.
usage: hicMergeLoops --inputFiles INPUTFILES [INPUTFILES ...] --outFileName OUTFILENAME --lowestResolution LOWESTRESOLUTION [--help] [--version]
- --inputFiles, -i
The loop files from hicDetectLoops. To use files from other sources, please follow ‘chr start end chr start end’ format and remove any header.
- --outFileName, -o
The name of the merged loop file.
- --lowestResolution, -r
The lowest resolution of all loop files, i.e. 5kb, 10kb and 25kb, please use 25000.
show program’s version number and exit