The quark contribution to the QCD pressure, $P_q$, is evaluated up to next-to-leading order (NLO) within the renormalization group optimized perturbation theory (RGOPT) resummation approach. To evaluate the complete QCD pressure we simply add the perturbative NLO contribution from massless gluons to the resummed $P_q$. Despite this unsophisticated approximation our results for $P=P_q+P_g$ at the central scale $M\sim 2\pi T$ show a remarkable agreement with lattice predictions for $0.25\lesssim T\lesssim 1\mathrm{GeV}$. We also show that by being imbued with RG properties, the RGOPT produces a drastic reduction of the embarrassing remnant scale dependence that plagues both standard thermal perturbative QCD and hard thermal loop perturbation theory applications.