Let $M_2(\Gamma_0(N))$ be the space of weight $2$ modular forms for $\Gamma_0(N)$ and let $f_1, \dots, f_r$ be a basis of normalized eigenforms for $M_2(\Gamma_0(N))$. Given a rational prime $p$, I'll say that the eigenforms $f_1, \dots, f_r$ are linearly independent mod $p$ if the mod $p$ reductions of the Fourier expansions of $f_1, \dots, f_r$ are $\mathbf{F}_{p^k}$-linearly independent in $\mathbf{F}_{p^k} [[ q ]]$.
My question is: what are sufficient conditions on $N$ and $p$ that ensure that $f_1, \dots, f_r$ are linearly independent mod $p$? For example, if $N$ is prime and $p$ divides the numerator of $(N-1)/12$, then $f_1, \dots, f_r$ are not linearly independent mod $p$, because there is an Eisenstein series which is congruent to a cuspform in $M_2(\Gamma_0(N))$. Are there sufficient conditions on $N$ and $p$ that ensure linear independence of the forms $f_1, \dots, f_r$?