We investigate the fundamental oscillation modes of anisotropic dark energy stars using the Chaplygin prescription, incorporating the effects of anisotropy via the Bowers-Liang model. We study the global properties of anisotropic dark energy stars by solving stellar structure equations, and also compare the obtained mass-radius profiles with observational data from gravitational wave events and millisecond pulsars. We employ the well known Cowling approximation to compute the prominent non-radial f-modes for different values of anisotropy parameters that ranges from \(-2\)  to \(+2\) . We report that f-mode spectra of anisotropic dark energy stars are distinct from those of neutron stars and quark stars. We strongly believe that this distinct nature of the spectra would result in easy identification in the near future. Moreover, the tidal deformability characteristics of anisotropic dark energy stars have also been investigated.