We performed detailed genetic analyses of the partial hemagglutinin–neuraminidase (HN) gene in 34 human respirovirus 3 (HRV3) strains from children with acute respiratory illness during 2013–2015 in Iwate Prefecture, Japan. In addition, we performed analyses of the evolutionary timescale of the gene using the Bayesian Markov chain Monte Carlo (MCMC) method. Furthermore, we analyzed pairwise distances and performed selective pressure analyses followed by linear B-cell epitope mapping and N-glycosylation and phylodynamic analyses. A phylogenetic tree showed that the strains diversified at around 1939, and the rate of molecular evolution was 7.6 × 10−4 substitutions/site/year. Although the pairwise distances were relatively short (0.03 ± 0.018 [mean ± standard deviation, SD]), two positive selection sites (Cys544Trp and Leu555Ser) and no amino acid substitutions were found in the active/catalytic sites. Six epitopes were estimated in this study, and three mouse monoclonal antibody binding sites (amino acid positions 278, 281, and 461) overlapped with two epitopes belonging to subcluster C3 strains. Bayesian skyline plot analyses indicated that subcluster C3 strains have been increasing from 2004, whereas subcluster C1 strains have declined from 2004. Based on these results, Iwate strains were divided into two subclusters and each subcluster evolved independently. Moreover, our results suggested that some predicted linear epitopes (epitopes 3 and 5) are candidates for an HRV3 vaccine motif. To better understand the details of the molecular evolution of HRV, further studies are needed.