*4.4. Phylogenetic Analysis*

Amino acid sequences corresponding to bacteriocin genes (class II and class III) found in this work, along with reference sequences identified by AntiSMASH 5.0 [57] and Bagel4 [58], and Uniprot databases were aligned using MAFFT [102]. Guidance2 [103] was used to filter unreliable positions and generate a mega alignment encompassing 5 alternative alignments for the sequences. The mega alignment was used to infer the evolutionary history of these proteins by using the Maximum Likelihood method, based on the VT model [104]. A discrete Gamma distribution was used to model evolutionary rate differences among sites, and the rate variation model allowed for some sites to be evolutionarily invariable [105]. Significance was assessed via aLRT [106]. All evolutionary analyses were conducted in PhyML 3.0 [107]. Tree visualization and annotation were performed on Interactive Tree Of Life (iTOL) v [108].
