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 ^thorder) harmonics of plateau region, only the electrons from the C side participate in the HHG. However, for its higher part (18 ^th–36 ^thorder), 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.