A dispersive representation based on unitarity and analyticity is used to study the low energy 𝛾𝑁 → 𝜋𝑁 and 𝛾 * 𝑁 → 𝜋𝑁 partial wave amplitudes. Final state interactions of the 𝜋𝑁 system are critical to this analysis. The left-hand cut contribution is estimated by invoking baryon chiral perturbation theory results, while the right-hand cut contribution responsible for final state interaction effects is taken into account via an Omnès formalism with elastic phase shifts as input. It is found that a good numerical fit can be achieved with only one subtraction parameter, and the experimental data of the multipole amplitudes 𝐸 0+ , 𝑆 0+ in the energy region below the Δ(1232) are well described when the photon virtuality 𝑄 2 ≤ 0.1GeV 2 . Furthermore, we extend the partial wave amplitudes to the second Riemann sheet to extract the couplings of the subthreshold resonance 𝑁 * (890). Its couplings extracted from the multipole amplitudes 𝐸 0+ , 𝑆 0+ are comparable to those of the 𝑁 * (1535) resonance.