Variability in the numbers of bacteria remaining in saline solution and whole milk following mild heat treatment has been studied with Listeria innocua, Enterococcus faecalis, Salmonella enterica serovar Enteritidis, and Pseudomonas fluorescens. As expected, the most heat-resistant bacterium was E. faecalis, while P. fluorescens was the least heat resistant, and all bacteria showed greater thermal resistance in whole milk than in saline solution. Despite the differences in the inactivation kinetics of these bacteria in different media, the variability in the final number of bacteria was affected neither by the species nor by the heating substrate, but it did depend on the intensity of the heat treatment. The more severe the heat treatment was, the lower the average number of surviving bacteria but the greater the variability. Our results indicated that the inactivation times for the cells within a population are not identically distributed random variables and that, therefore, the population includes subpopulations of cells with different distributions for the heat resistance parameters. A linear relationship between the variability of the log of the final bacterial concentration and the logarithmic reduction in the size of the bacterial population was found.The safety and quality of many foods, such as minimally processed foods and some ready-to-eat products, are affected by a combination of multiple sublethal stresses due to processing and storage conditions (17, 32). The use of mild preservative processes instead of harsher ones affects the microbiological risk. It is therefore important to identify potential hazards and to perform accurate quantitative microbial risk assessments including probabilistic modeling approaches (17).It is unlikely that food products are contaminated by one unique bacterial strain or species. In fact, even cells from a pure culture derived from one cell are not necessarily identical (25). Nauta (16) differentiated between uncertainty, which is due to the lack of perfect knowledge about a parameter and which is reducible by further measurements, and variability, which reflects the true variability of a population and is therefore irreducible. He showed the importance of both parameters in quantitative microbial risk assessment models. The estimation of variability usually refers to the variation in the responses of the individual cells within the population, which may be homogeneous or heterogeneous. In a homogeneous population, the heat resistance parameters of the individual cells are identically distributed random variables. Heterogeneous populations include subpopulations of cells with different distributions for the heat resistance parameters. The heterogeneity of the individual cells or spores within a population is the most accepted explanation for the tails of inactivation curves (1, 30) and unexpected increases in the variability of the number of survivors after heat treatments (7,8,30).The microbial concentration in a food product at a given time is affected by the initial...