Heterogeneity formation within biofilm systems