Hello, I'mworking on a multiomic single cell RNAseq and ATAC object and I am interested in finding DARs using FindMarkers with Signac. I ran FindMarkers(object, group.by="treatment_timepoint_cluster", ident.1= "2_wt_d8", ident.2 = "2_ko_d8", min.pct= 0.05, latent.vars= "peak_region_fragments", assay= "ATAC), and I visualize the DARs in coverage plots but there are any differences . The fold changes are around 0.6 to -0.6. I don't understand why I dont see any difference. My question Is: can I find peaks with macs2 grouping by conditions and then create a new assay to run findmarkers again skipping the integration steps of this assay? Or I Need to integrated the assay peaks to the multiomic object before trying to apply findmarkers? Thanks