Figure 1

The UniPeak workflow. A: Sequence reads are considered as hits at their 5^′ start positions, strand-specifically. B: A global read-shift value is computed independently for each sample to align forward and reverse reads. C: The shifted reads from all samples are then used to estimate a single underlying density profile. Enriched regions are identified where this density exceeds a fixed threshold, determined as a function of sequencing depth and genome size. Shifted reads from each sample are counted within these regions, providing a read count for each sample within each genomic region.