We present a high statistics study of the isovector nucleon charges and form factors using seven ensembles of $2+1$-flavor Wilson-clover fermions. The axial vector and pseudoscalar form factors obtained on each of these ensembles satisfy the partially conserved axial current relation between them once the lowest energy $N\ensuremath{\pi}$ excited state is included in the spectral decomposition of the correlation functions used for extracting the ground state matrix elements. Similarly, we find evidence that the $N\ensuremath{\pi}\ensuremath{\pi}$ excited state contributes to the correlation functions with the insertion of the vector current, consistent with the vector meson dominance model. The resulting form factors are consistent with the Kelly parametrization of the experimental electric and magnetic data. Our final estimates for the isovector charges are ${g}_{A}^{u\ensuremath{-}d}=1.32(6)(5{)}_{\mathrm{sys}}$, ${g}_{S}^{u\ensuremath{-}d}=1.06(9)(6{)}_{\mathrm{sys}}$, and ${g}_{T}^{u\ensuremath{-}d}=0.97(3)(2{)}_{\mathrm{sys}}$, where the first error is the overall analysis uncertainty and the second is an additional combined systematic uncertainty. The form factors yield: (i) the axial charge radius squared, ${⟨{r}_{A}^{2}⟩}^{u\ensuremath{-}d}=0.428(53)(30{)}_{\mathrm{sys}}\text{ }\text{ }{\mathrm{fm}}^{2}$; (ii) the induced pseudoscalar charge, ${g}_{P}^{*}=7.9(7)(9{)}_{\mathrm{sys}}$; (iii) the pion-nucleon coupling, ${g}_{\ensuremath{\pi}\mathrm{NN}}=12.4(1.2)$; (iv) the electric charge radius squared, ${⟨{r}_{E}^{2}⟩}^{u\ensuremath{-}d}=0.85(12)(19{)}_{\mathrm{sys}}\text{ }\text{ }{\mathrm{fm}}^{2}$; (v) the magnetic charge radius squared, ${⟨{r}_{M}^{2}⟩}^{u\ensuremath{-}d}=\phantom{\rule{0ex}{0ex}}0.71(19)(23{)}_{\mathrm{sys}}\text{ }\text{ }{\mathrm{fm}}^{2}$; and (vi) the magnetic moment, ${\ensuremath{\mu}}^{u\ensuremath{-}d}=4.15(22)(10{)}_{\mathrm{sys}}$. All our results are consistent with phenomenological/experimental values but with larger errors. Last, we present a Pad\'e parametrization of the axial, electric, and magnetic form factors over the range $0.04<{Q}^{2}<1\text{ }\text{ }{\mathrm{GeV}}^{2}$ for phenomenological studies.