The effect of the existence of several prototropic tautomers of cytosine on its UV/vis spectra and the excited state decay dynamics is studied by spectral and nonadiabatic molecular dynamics (NAMD) simulations in connection with the mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) method. Simulated UV/vis spectra provide a strong indication that the H3N keto–amino cytosine tautomer (the least anticipated species) may be present under experimental conditions. The NAMD simulations yield a wide range of excited state decay constants for various tautomers of cytosine, ranging from ∼1.3 ps for the biologically relevant H1N keto–amino tautomer to ∼0.1 ps for the keto–imino tautomer. The slowness of the H1N decay dynamics follows from the presence of a barrier on the excited state energy surface separating the Franck–Condon structure from the major decay funnel, the conical intersection seam. It is suggested that the experimentally observed photodecay dynamics may result from a combination of the decay processes of various tautomers (H3N in particular) present simultaneously under the experimental conditions.