The probability of breakage in service lifetime of heat-tempered glass panes contaminated by nickel sulfide inclusions is estimated with a multiscale micromechanically motivated statistical theory, which considers the effects of the heat soak test (HST). Short and long HSTs differently affect the phase transformation of NiS of diverse chemical composition, whose increase in volume can break the glass. The main hypothesis, corroborated by experiments, is that there is a lower limit for the size of NiS stones below which no crack can be initiated from the volumetric expansion. The catastrophic propagation of nucleated fractures in the long term is modeled through a rescaled critical stress intensity factor, which accounts for the subcritical crack propagation and the slow phase transformation of NiS. A parametric analysis evidences how the failure probability is strongly affected by these parameters, depending on the holding time in the HST. Tailored experimental activity is suggested for the proper calibration of the model.