In this talk I will consider asymptotically free gauge theories with gauge group $SU(N_c)$ and $N_f$ quarks with mass $m_q <<\Lambda_{QCD}$ that undergo chiral symmetry breaking and confinement. I will described a proposal for a bootstrap method to compute the S-matrix of the pseudo-Goldstone bosons (pions) that dominate the low energy physics. For the important case of $N_c=3$, $N_f=2$, a numerical implementation of the method gives the phase shifts of the $S0$, $P1$ and $S2$ waves in good agreement with experimental results. The method incorporates gauge theory information ($N_c$, $N_f$, $m_q$, $\Lambda_{QCD}$) by using the form-factor bootstrap recently proposed by Karateev, Kuhn and Penedones together with a finite energy version of the SVZ sum rules. This requires, in addition, the values of the quark and gluon condensates. At low energy we impose constraints from chiral symmetry breaking which additionally require knowing the pion mass $m_\pi$.