Compared with nonpolar molecules, owing to the inherent asymmetry, polar molecules exhibit rich and very complex electronic dynamics under the interaction with strong laser fields. In this work, high-order harmonic generation (HHG) of polar molecules CO is investigated by using the three-dimensional time-dependent Hartree-Fock (3D-TDHF) theory, with all electrons active. Through the high harmonic spectra and time-frequency analyses, it is found that when the laser field polarizes along the molecular axis, the ionized electrons from the two sides (C side and O side) contribute differently to the harmonic radiation. On the one hand, the harmonic intensity from the C side is greater than that from the O side, which is caused by the ionization rate. On the other hand, for the lower-order (7^th–17^th order) harmonics of plateau region, only the electrons from the C side participate in the HHG. However, for its higher part (18^th–36^th order), the electrons from both C side and O side contribute to high harmonics simultaneously. Moreover, the difference between contributions from two sides is related to the alignment angle θ between the laser polarization and the molecular axis, and it reaches a maximum value around θ = 0º and a minimum value around θ = 90º. There are two strong resonances around harmonic order H12.6 (19.5 eV) and H18 (27.9 eV) in the harmonic spectra when θ = 0º. The first resonance around H12.6 reveals that part of electrons ionized from the C side recombine to the vicinity of the further O nucleus. Near the second resonance around H18, there appears a shape resonance. Nevertheless, the shape resonances from the C and O sides are disparate. Based on the strong-field approximation theory, the ratio between photoionization cross sections from C and O sides around the shape resonance is calculated. The ratio is about 5.5 from 3D-TDHF, which is greater than the result of 3 simulated by ePloyScat, where only HOMO is considered. This discrepancy reveals that multi-electron effects enhance the asymmetry of polar molecules. This work provides an in-depth insight into the asymmetry in HHG of polar molecules, which benefits the generation of isolated attosecond pulse . It also promotes the application of high harmonic spectra in tracking the ultrafast dynamics of electrons.