We undertake a covariance error analysis of the pion-nucleon-nucleon coupling constants from the Granada-2013 np and pp database comprising a total of 6720 scattering data below LAB energy of 350 MeV. Assuming a unique pion-nucleon coupling constant in the One Pion Exchange potential above a boundary radius r c = 3fm we obtain f 2 = 0.0763(1). The effects of charge symmetry breaking on the 3 P 0 , 3 P 1 and 3 P 2 partial waves are analyzed and we find f 2 p = 0.0761(4), f 2 0 = 0.0790(9) and f 2 c = 0.0772(6) with a strong anti-correlation between f 2 c and f 2 0 . We successfully test normality for the residuals of the fit. Potential tails in terms of different boundary radii as well as chiral Two-Pion-Exchange contributions as sources of systematic uncertainty are also investigated.