To develop mathematical models describing lag times of individual bacterial cells (tau), experimental tau data were fitted to a variety of continuous distributions using BestFit. Six strains of Escherichia coli O157:H7 were used, and serial dilutions were made in Bioscreen multi-well plates to get single cells per well. Detection times (td) for individual wells were converted to tau using the maximum specific growth rate (mu) for each strain. All strains were subject to 25 trials, with up to 100 replicate wells per trial. Some strains had significantly longer t(d), and lower mu, but the tau values were not significantly different. Distributions were best fit in the order Pearson V>Pearson VI>Extreme Value>Lognormal>Lognormal2>Inverse Gaussian based on the Anderson-Darling statistic. The Lognormal distribution was selected because there was less variability in the parameter values, and parameters have specific biological meanings. Distributions could be fit to sample populations as low as six, with fittings and parameter values comparable to those obtained with larger samples (up to 89). Extreme Value, Pearson V, and Pearson VI distributions were more suitable for fitting tau values generated from a Lognormal distribution when the numbers of sample points were few, which suggested that there are similarities between the distributions. The results suggest that a Lognormal distribution can be used successfully to characterize tau.