Before the analysis, raw reads were trimmed using Trimmomatic 0.38 [29] by removing adapter sequences and ambiguous nucleotides. Reads with quality scores <20 and length below 30 bp were removed. The resulting high-quality sequences of the 9 transcriptomes were used for the detection of SNPs in...