In the output file of MACS2, the fold enrichment normally refers to enrichment of summit to local lambda model. I am using
--nolambda option for ATAC-seq analysis (recommended as there is no input file like Chip-seq). Is the output file now showing enrichment with respect to global lambda? If so, does anyone have a suggestion for a reasonable cutoff?