design: Infected gonads were pooled and DNA was extracted from the infected snail gonad using the DNeasy Blood &amp; Tissue Kit (QIAGEN, Cat#69504) according to the manufacturers protocol. ITS2 region was amplified and ITS2_F and ITS2_R primers (Hammoud, C. et al., 2022) were used for amplification. Successfully amplified PCR products were used for library preparation and library was prepared using Twist MF Library Prep Kit. Prepared libraries were quantified using Qubit DNA HS Assay (Invitrogen, Cat# Q32854). The obtained libraries were pooled and diluted to the final optimal loading concentration before cluster amplification on the Illumina flow cell. Once the cluster generation was completed, the cluster flow cell was loaded on the Illumina Miseq instrument.