Control-1.R
The libraries were sequenced on an llumina Novaseq 6000 platform and 150 bp paired-end reads were generated. About  raw reads for each sample were generated. Raw reads of fastq format were firstly processed using fastp1 and the low quality reads were removed to obtain the clean reads.