We present a precise lattice QCD calculation of the contribution to the neutron-proton mass splitting arising from strong isospin breaking, mn−mp|QCD=2.32±0.17 MeV. We also determine mΞ−−mΞ0|QCD=5.44±0.31 MeV. The calculation is performed at three values of the pion mass, with several values of the quark mass splitting and multiple lattice volumes, but only a single lattice spacing and an estimate of discretization errors. The calculations are performed on the anisotropic clover-Wilson ensembles generated by the Hadron Spectrum Collaboration. The omega-baryon mass is used to set the scale a−1t=6111±127 MeV, while the kaon masses are used to determine the value of the light-quark mass spitting. The nucleon mass splitting is then determined as a function of the pion mass. We observe, for the first time, conclusive evidence for non-analytic light quark mass dependence in lattice QCD calculations of the baryon spectrum. When left as a free parameter, the fits prefer a nucleon axial coupling of gA=1.24(56). To highlight the presence of this chiral logarithm in the nucleon mass splitting, we also compute the isospin splitting in the cascade-baryon system which is less sensitive to chiral dynamics. Finally, we update the best lattice QCD determination of the CP-odd pion-nucleon coupling that would arise from a non-zero QCD theta-term, g¯0/(2‾√fπ)=(14.7±1.8±1.4)⋅10−3θ¯. The original lattice QCD correlation functions, analysis results and extrapolated quantities are packaged in HDF5 files made publicly available including a simple Python script to access the numerical results, construct effective mass plots along with our analysis results, and perform the extrapolations of various quantities determined in this work.