Phages drive bacterial diversity -profoundly influencing diverse microbial communities, from 14 microbiomes to the drivers of global biogeochemical cycling. The vast genomic diversity of phages is 15 gradually being uncovered as >8000 phage genomes have now been sequenced. Aiming to broaden our 16 understanding of Escherichia coli (MG1655, K-12) phages, we screened 188 Danish wastewater samples (0.5 17 ml) and identified 136 phages of which 104 are unique phage species and 91 represent novel species, 18 including several novel lineages. These phages are estimated to represent roughly a third of the true diversity 19 of Escherichia phages in Danish wastewater. The novel phages are remarkably diverse and represent four 20 different families Myoviridae, Siphoviridae, Podoviridae and Microviridae. They group into 14 distinct clusters 21 and nine singletons without any substantial similarity to other phages in the dataset. Their genomes vary 22 drastically in length from merely 5 342 bp to 170 817 kb, with an impressive span of GC contents ranging 23 from 35.3% to 60.0%. Hence, even for a model host bacterium, in the go-to source for phages, substantial 24 diversity remains to be uncovered. These results expand and underlines the range of Escherichia phage 25 diversity and demonstrate how far we are from fully disclosing phage diversity and ecology. 26 27 28 108 6 R studio version 1.1.456 [48] was used for rarefaction, species diversity (q = 0, datatype: incidence_raw), 109 extrapolation hereof (estimadeD) and estimation of sample coverage. The visualisation of genome sizes and 110 GC contents was prepared in Excel version 16.31.111 3. Results and Discussion 112 3.2. Phylogenetics, taxonomy and species richness 113 The sequenced phages were analysed strictly in silico, focusing on their relatedness to known phages, 114 their taxonomy and any distinctive characteristics. Based on the confirmed morphology of closely related 115 phages, at least four different families are represented, Myoviridae (58%), Siphoviridae (25%), Podoviridae (7.4%) 116 and even the single stranded DNA (ssDNA) Microviridae (5.9%) (Table 1). Five (3.7%) of the phages are so 117 distinct that they could not with certainty be assigned to any family. A similar distribution was found by Korf 118 et al., (2019), i.e. 70% Myoviridae, 22% Siphoviridae and 8% Podoviridae, just as an analysis of the genomes of 119 Caudovirales infecting Enterobacteriaceae, by Grose & Casjens (2014), also identified more clusters belonging to 120 the Myoviridae, than the Siphoviridae and fewest of the Podoviridae [8,22]. Jurczak-Kurek et al., (2016) found more 121 siphoviruses than myovirus, but also found the Podoviridae to be the least abundant [23]. Nonetheless, these 122 distributions are just as likely to be caused by culture or isolation methods, as they are to reflect true 123 abundances. Based on DNA homology and phylogenetic analysis, the 104 unique phages identified in this 124 study group into 14 distinct clusters and nine singletons, having a <1% Gegenees score...