Hello,

I am using scATACseq from the catlas. I would like to select some open chromatin regions (OCRs) in given cell types. Catlas provide these OCRs as follow:

#chrom  start end id score
chr1    180570  181069  cCRE_2  65.884831996234
chr1    633780  634279  cCRE_7  90640.346231263
chr1    778387  778886  cCRE_8  369.014168253344
chr1    827320  827819  cCRE_14 50.392152324195
chr1    904491  904990  cCRE_24 29.9988771375576
chr1    940078  940577  cCRE_34 13.7404112137246
chr1    1019289 1019788 cCRE_46 56.5497501656823
chr1    1019854 1020353 cCRE_47 16.0218937141982

From what I grasped, these scores were obtained with Macs2 and would correspond to -log10(pvalue). But the values seem quite high to me... Am I missing something there ? And if I am not, is there a general recommendation on a threshold to apply to these scores ?

Thanks in advance,

Raphaël