*4.4. Sequencing and Phylogenetic Analysis*

Positive PCR amplicons were two-directionally sequenced at GENEWIZ (Suzhou, China). Sequences were assembled and edited using DNAstar Lasergene Editseq 7.1.0 (http://www.dnastar.com/), and reference sequences downloaded from the GenBank database were compared to determine the genotype and subtype of *Cryptosporidium* using ClustalX 2.1 (http://www.clustal.org/). The established nomenclature system was used in the naming subtype of *C. parvum* [11]. Phylogenetic analyses were conducted using neighbor-joining methods based on the Kimura-2 parameter model in MEGA 7.0 (http://www.megasoftware.net/). Seven presentative nucleotide sequences obtained in this study were submitted in the GenBank database (https://www.ncbi.nlm.nih.gov/) under the accession numbers: MH442993–MH442996, MT703861, MT703862, and MT724047.
