  We compute the renormalized Dirac spectral density in $N_f = 2+1$ QCD at
physical quark masses, temperature $T = 230$ MeV and system size $L_s = 3.4$
fm. To that end, we perform a point-wise continuum limit of the staggered
density in lattice QCD with staggered quarks. We find, for the first time, that
a clear infrared structure (IR peak) emerges in the density of Dirac operator
describing dynamical quarks. Features of this structure are consistent with
those previously attributed to the recently-proposed IR phase of thermal QCD.
Our results (i) provide solid evidence that these IR features are stable and
physical; (ii) improve the upper bound for IR-phase transition temperature
$T_{\mathrm{IR}}$ so that the new window is $200 < T_{\mathrm{IR}} < 230\,$MeV;
(iii) support non-restoration of anomalous U$_{\mathrm A}$(1) symmetry (chiral
limit) below $T \!=\! 230\,$MeV.
