  We study the non-radial $f$-mode oscillations of both isotropic and
anisotropic dark energy stars by using the modified Chaplygin prescription of
dark energy to model the stellar matter. The anisotropic pressure in the system
is modeled with Bowers-Liang prescription. By solving the stellar structure
equations in presence of anisotropy, we study the global properties of the dark
energy star and compare the mass-radius profiles with data from GW events and
milli-second pulsars. We proceed to determine the prominent non-radial $l=2$
$f$-mode frequencies of the anisotropic dark energy star by employing the
Cowling approximation and analyse and quantify the spectra by varying the
anisotropic parameter. We report that $f$-mode spectra of dark energy star have
distinctly different behaviour compared to neutron star and quark star, and
this may possibly help in its future identification. Further, the tidal
deformability factors of the anisotropic dark energy stars have also been
analyzed.
