Abstract Aims This study was conducted to investigate the presence of Shiga toxin-producing O157 and non-O157 E. coli in raw water buffalo milk, as well as to determine the virulence gene profiles, phylogroups, sequence types, and serotypes of the isolated strains. Methods and results A total of 200 hand-milked raw water buffalo milk samples were collected from 200 different water buffaloes over a period of three months from 20 different farms. Isolation of STEC was performed using CHROMagar STEC. Presence of stx1, stx2, and eaeA genes were investigated by mPCR. Phylogroups and sequence types of E. coli strains were determined by Clermont phylotyping and MLST. Serotyping was performed using PCR or WGS. According to the results, two milk samples obtained from two different farms were found as STEC-positive. All Stx-positive E. coli isolates belonged to phylogenetic group A and were assigned to ST10. WGS results indicated that serotype of two isolates was O21:H25 and average nucleotide identity was detected at 99.99%. Thirteen additional registered E. coli O21:H25 assembled WGS data were obtained from EnteroBase and a phylogenetic tree was constructed. Conclusions With this study, the presence of stx2 harboring E. coli O21:H25 in milk was identified for the first time. Although the identified serotype is considered a non-pathogen seropathotype, we conclude it could play an important role in the environmental circulation of Stx-phages and consequently contribute to the emergence of new STEC-related outbreaks.