Ultralight axions and other bosons are dark matter candidates present in many high energy physics theories beyond the Standard Model. In particular, the string axiverse postulates the existence of up to \U0001d4aa(100) light scalar bosons constituting the dark sector. We test the validity of the effective field theory of large-scale structure approach to mixed ultralight axion dark matter by making our own mock galaxy catalogs and find an anisotropic ultralight axion signature in the galaxy quadrupole. We use the Fourier-space galaxy clustering statistics from the Baryon Oscillation Spectroscopic Survey (BOSS) and obtain upper bounds for the axion relic density Ω_a h2<0.004 for axions of mass 10-31 eV≤ m_a≤ 10-26 eV at 95% confidence. We also improve existing constraints by a factor of over 4.5 and 2.1 for axion masses of 10-25 eV and 10-32 eV compared to CMB alone. To improve the convergence speed of our analysis, we develop an augmented interpolation scheme allowing a fast computation of the axion contribution to the linear matter power spectrum. This method leads to a 70% reduction of the computational cost for the full Monte Carlo Markov chains analysis.