For the previously published HSF123 and H9-Laurent24 data sets, the GEO sequence read data were mapped using the MethylC-Seq pipeline (H9-Laurent) and BS Seeker (HSF1)36 (settings: -e ≤ 55, -m 3), and post-processing and methylcytosine identification was performed using MethylC-Seq pipeline as described earlier.