Whole-genome sequencing was performed using the Illumina HiSeq 2500 platform. Genomic DNA (gDNA) was used for paired-end libraries (2 × 150 bp) with NEBNext® Fast DNA Fragmentation and Library Preparation Kit (New England Biolabs Inc., MA, USA). The final library quality was assessed using the 2100 Bioanalyzer (Agilent Technologies, CA, USA) and agarose gel electrophoresis. We performed the quality control of raw sequences obtained from Illumina HiSeq sequencing by using FastQC v0.12.1, considering as good sequences that ones with phred score equal or higher than 30. These sequences underwent de-novo assembly using SPAdes v3.15.5. We employed QUAST v4.4 to measure the assembled genomes' quality metrics and perform statistical analysis. Pair-end reads were aligned to the assemblies using the Burrows-Wheeler Aligner (BWA) tool. The average sequencing coverage was then calculated from the resulting BAM file using SAMtools. Taxonomic classification was first performed with GTDB-tk. When GTDB-tk was insufficient for detailed species taxonomy, FastANI v1.33 was implemented as a supplementary tool. When species-level classification was not achieved using GTDB-tk, we downloaded all available assemblies for each identified genus from NCBI. Species were considered the same if their Average Nucleotide Identity (ANI) values exceeded 95%. Finally, the annotation of functional genes in each genome was carried out using RAST-tk.