Genomic Technologies shared resource on two lanes of Illumina HiSeq 2000 or HiSeq 4000 obtaining 2 × 125 or 2 × 151 single- or paired-end reads. Since only eight samples were sequenced in the single-end mode and all showed different epigenomics profile based on multidimensional scaling analysis, we decided to exclude them. After additional quality controls (see below), 269 ATAC-seq libraries were retained21.