Code from the manuscript. The raw sequence data reported in this paper have been deposited in the Genome Sequence Archive (Genomics, Proteomics & Bioinformatics 2021) in National Genomics Data Center (Nucleic Acids Res 2022), China National Center for Bioinformation / Beijing Institute of Genomics, Chinese Academy of Sciences (GSA: CRA015021, GSA-Human: HRA006754) that are publicly accessible at https://ngdc.cncb.ac.cn/gsa. If you have any question or suggestion, please contect lilab_in_siat@163.com.