We study the zeros of sections of the form $T_k s_k$ of a large power $L^{\otimes k} \to M$ of a holomorphic positive Hermitian line bundle over a compact K\''ahler manifold $M$, where $s_k$ is a random holomorphic section of $L^{\otimes k}$ and $T_k$ is a Berezin-Toeplitz operator, in the limit $k \to +\infty$. In particular, we compute the second order approximation of the expectation of the distribution of these zeros. In a ball of radius of order $k^{-\frac{1}{2}}$ around $x \in M$, assuming that the principal symbol $f$ of $T_k$ is real-valued and vanishes transversally, we show that this expectation exhibits two drastically different behaviors depending on whether $f(x) = 0$ or $f(x) \neq 0$. These different regimes are related to a similar phenomenon about the convergence of the normalized Fubini-Study forms associated with $T_k$: they converge to the K\''ahler form in the sense of currents as $k\rightarrow + \infty$, but not as differential forms (even pointwise). This contrasts with the standard case $f=1$, in which the convergence is in the $\mathscr{C}^{\infty}$-topology. From this, we are able to recover the zero set of $f$ from the zeros of $T_k s_k$.