Human coronaviruses (HCoV) are RNA viruses that cause respiratory tract infections with viral replication of limited duration. The host and viral population heterogeneity could influence clinical phenotypes. Employing long RT-PCR with Illumina sequencing, we quantified the gene mutation load at 0.5% mutation frequency for the 4529 bp-domain spanning the Spike gene (4086 bp) of HCoV-OC43 in four upper respiratory clinical specimens obtained during acute illness. There were a total of 121 mutations for all four HCoV samples with the average number of mutations at 30.3 ± 10.2, which is significantly higher than that expected from the Illumina sequencing error rate. There were two mutation peaks, one at the 5' end and the other near position 1 550 in the S1 subunit. Two coronavirus samples were genotype B and two were genotype D, clustering with HCoV-OC43 strain AY391777 in neighbor-joining tree phylogenetic analysis. Nonsynonymous mutations were 76.1 ± 14% of mutation load. Although lower than other RNA viruses such as hepatitis C virus, HCoV-OC43 did exhibit quasi-species. The rate of nonsynonymous mutations was higher in the HCoV-OC43 isolates than in hepatitis C (HCV) virus genotype 1a isolates analyzed for comparison in this study. These characteristics of HCoV-OC43 may affect viral replication dynamics, receptor binding, antigenicity, evolution, transmission, and clinical illness.