Gravitational waves (GWs) emanating from unstable quasi-normal modes in neutron stars (NSs) could be accessible with the improved sensitivity of the current GW detectors or with the next-generation GW detectors and, therefore, can be employed to study the NS interior. Assuming f-mode excitation in isolated pulsars with typical energy of pulsar glitches and considering potential f-mode GW candidates for A+ (upgraded LIGO detectors operating at fifth observing run design sensitivity) and Einstein Telescope (ET), we demonstrate the inverse problem of NS asteroseismology within a Bayesian formalism to constrain the nuclear parameters and NS equation of state (EOS). We describe the NS interior within relativistic mean-field formalism. Taking the example of glitching pulsars, we find that for a single event in A+ and ET, among the nuclear parameters, the nucleon effective mass (m*) within 90% credible interval can be restricted within 10% and 5%, respectively. At the same time, the incompressibility (K) and the slope of the symmetry energy (L) are only loosely constrained. Considering multiple (10) events in A+ and ET, all the nuclear parameters are well constrained, especially m*, which can be constrained to 3% and 2% in A+ and ET, respectively. Uncertainty in the observables of a 1.4 M
⊙ NS such as radius (
R
1.4
M
⊙
), f-mode frequency (
f
1.4
M
⊙
), damping time (
τ
1.4
M
⊙
), and a few EOS properties including squared speed of sound (c
s
2) are also estimated.