design: The pooled pulldown DNA (~10ng determined by Agilent TapeStation4200) samples were submitted to Welgene Co. (Taiwan) for NGS conducted on Illumina NovaSeq (151 cycles paired-end reads). The NGS library was constructed with an SSEL XT Low Input Reagent Kit with IDT for Illumina  TruSeq DNA UD Indexes following the manufacturers protocol by Welgene Co. core facility. The sequencing experiments were conducted once for Human genome. Raw FASTQ reads were quality and adaptor trimmed by fastp 0.20.0. Trimmed reads were mapped to the human reference genome GRCh38 by bowtie2 2.3.4.3. The control reads provided background information for deducing the folds of peak enrichment by MACS 2.2.4 with minimum FDR (q-value) cutoff value set to 0.05. The distribution of genomic elements (Promoter, 5'-UTR, 3'-UTR, Exon, Downstream, Intron, Distal Intergenic) and the percentage of GC (GC%) in the reference genome and enriched peaks (>1.36 folds, q<0.05) were analyzed in R 3.5.1.