We theoretically study the dependence of below-threshold harmonic generation (BTHG) of atoms on the carrier-envelope phase (CEP) driven by few-cycle mid-infrared laser pulses. The BTHG spectra can be accurately and efficiently calculated by solving the three-dimensional time-dependent Schrödinger equation using the time-dependent generalized pseudospectral method. We present the BTHG spectra as a function of the laser-field CEP. CEP-dependent enhancement or suppression occurred at low laser field intensities owing to the changes in the resonant effects associated with multiple quantum trajectories. However, the BTHG of atoms driven by high laser intensities is insensitive to the CEP. The synchrosqueezing time-frequency transform of the BTHG and extended semiclassical analysis are performed to elucidate the underlying physical mechanism.