In typical statistical mechanical systems the grand canonical partition function at finite volume is proportional to a polynomial of the fugacity $e^{\mu/T}$. The zero of this Lee-Yang polynomial closest to the origin determines the radius of convergence of the Taylor expansion of the pressure around $\mu=0$. The computationally cheapest formulation of lattice QCD, rooted staggered fermions, with the usual definition of the rooted determinant, does not admit such a Lee-Yang polynomial. We argue that the radius of convergence is then bounded by the spectral gap of the reduced matrix of the unrooted staggered operator. This is a cutoff effect that potentially affects all estimates of the radius of convergence with the standard staggered rooting. We suggest a new definition of the rooted staggered determinant at finite chemical potential that allows for a definition of a Lee-Yang polynomial, and, therefore of the numerical study of Lee-Yang zeros. We also describe an algorithm to determine the Lee-Yang zeros and apply it to configurations generated with the 2-stout improved staggered action at $N_t = 4$. We perform a finite-volume scaling study of the leading Lee-Yang zeros and estimate the radius of convergence of the Taylor expansion extrapolated to an infinite volume. We show that the limiting singularity is not on the real line, thus giving a lower bound on the location of any possible phase transitions at this lattice spacing. In the vicinity of the crossover temperature at zero chemical potential, the radius of convergence turns out to be $\mu_B/T \approx 2$ and roughly temperature independent. Our simulations are performed at strange quark chemical potential $\mu_s=0$, but the method can be straightforwardly extended to strangeness chemical potential $\mu_S=0$ or strangeness neutrality.