Hi everyone!

I have DamID data for two conditions in the form of the log2 methylation ratio of Dam LaminB1 to Dam Only, from this publication (https://genome.cshlp.org/content/27/10/1634.full).

The distribution of windows for the data are in and around 1mb, and there are about 28000 data points.

head(DamID_OIS)
chr   start     end      ratio
chr1  764621  864620 -0.7687818
chr1  864621  964620 -0.9207849
chr1  964621 1064620 -0.9493541
chr1 1064621 1164620 -0.8203121
chr1 1164621 1264620 -1.6871681
chr1 1264621 1364620 -1.6889285

From what I understand, positive values would mean that the region is attached to the lamina, based on the methylation ratio. But the range of the positive values is quite large (~15000). From this list, how can I call boundaries of Lamina Associated Domains? Should I just take higher values based on some quantile percentage, or asses the AT content of the regions and concatenate high AT regions? Or is there a tool/caller that does the job?

Thanks in advance!

Rishi