The resulting data were compared with the reference genome assembly GCF_022869665.1 (
According to previous studies [31 (link),32 (link)], we annotated the unmatched genes with the NR database using Rockhopper. Unannotated genes with 50–500 nt of length and stable secondary structures were listed as candidate sRNAs [46 (link)]. In addition, we used the sRNAMap database (version 2009) and the Rfam database (version 13) to match known sRNAs.