design: Illumina MiSeq PE250 system was used. Reads with phred scores less than 25 were trimmed out. Using bowtie2, trimmed reads were mapped backed to H. smyrnensis AAD6 genome (RefSeq: NZ_AJKS00000000.2). VarScan was used to predict the structural variants on sequenced genomes (default settings).