Okay so this question was really bugging me so I spent the afternoon learning how to do the approximation.
The general Poschl-Teller Potential is $V(x) = \frac{V_o}{\cosh^2(\alpha \cdot (x-x_0))} $
I take $ V_0 = 1, \alpha =1, x_0 = 0 \ \ \ \ \rightarrow \ \ \ \ V(x) = \frac{1}{\cosh^2(x)} = \rm{sech}^2(x)$
Schutz and Will derive a general first order WKB approximation for QNMs:
$ \frac{Q_0}{\sqrt{2Q_0''}} = i(n+\frac{1}{2})$
This expands to:
$ \frac{w^2-\text{sech}^2(x)}{\sqrt{2} \sqrt{2 \text{sech}^4(x)-4 \tanh ^2(x) \text{sech}^2(x)}} = i(n+\frac{1}{2}) $
I evaluate at the peak of the potential, $x=0$, and solve for w taking the negative solution:
$ w = -\sqrt[4]{-1} \sqrt{2 n+(1-i)} $
I also reflect it across the imaginary axis to get the following spectrum, (vertical axis is the Im(w), and horizontal axis is Re(w)):
Compare this against the analytical derivation included by Berti, Cardoso, and Starinets:
$ \omega = \sqrt{1-\frac{1}{4}}-\frac{1}{2} i (2 n+1) $
Reflecting this as well yields the following spectrum, (again vertical axis is the Im(w), and horizontal axis is Re(w)):
As is known the first order WKB is not terribly acurate, the fundamental mode is close, $-1.09868 - 0.45509i$ vs $-0.866025 - 0.5i$, but the asymptotic behavior is different, and there is significant deviation in the real value for the overtones.