We investigate the impact of the Tsallis nonextensive statistics introduced by intrinsic temperature fluctuations in $p$-Air ultrahigh energy interactions on observables of cosmic ray showers, such as the slant depth of the maximum $X_{max}$ and the muon number on the ground $N_{\mu }$. The results show that these observables are significantly affected by temperature fluctuations and agree qualitatively with the predictions of Heitler model.