We introduce and analyze a prototype model for chemotactic effects in biofilm formation. The model is a system of quasilinear parabolic equations into which two thresholds are built in. One occurs at zero cell density level, the second one is related to the maximal density which the cells cannot exceed. Accordingly, both diffusion and taxis terms have degenerate or singular parts. This model extends a previously introduced degenerate biofilm model by combining it with a chemotaxis equation. We give conditions for existence and uniqueness of weak solutions and illustrate the model behavior in numerical simulations.