Pseudomonas aeruginosabiofilms are relevant for a variety of disease settings, including pulmonary infections in people with cystic fibrosis. Biofilms are initiated by individual bacteria that undergo a phenotypic switch and produce an extracellular polymeric slime (EPS). However, the viscoelastic characteristics of biofilms at different stages of formation and the contributions of different EPS constituents have not been fully explored. For this purpose, we develop and parameterize a mathematical model to study the rheological behavior of three biofilms —P. aeruginosawild type PAO1, isogenic rugose small colony variant (RSCV), and mucoid variant biofilms against a range of experimental data. Using Bayesian inference to estimate these viscoelastic properties, we quantify the rheological characteristics of the biofilm EPS. We employ a Monte Carlo Markov Chain algorithm to estimate these properties ofP. aeruginosavariant biofilms in comparison to those of wild type. This information helps us understand the rheological behavior of biofilms at different stages of their development. The mechanical properties of wild type biofilms change significantly over time and are more sensitive to small changes in their composition than the other two mutants.