In this article we present an $L_p$-theory ($p\geq 2$) for the time-fractional quasi-linear stochastic partial differential equations (SPDEs) of type $$ \partial^α_tu=L(ω,t,x)u+f(u)+\partial^β_t \sum_{k=1}^{\infty}\int^t_0 ( Λ^k(ω,t,x)u+g^k(u))dw^k_t, $$ where $α\in (0,2)$, $β<α+\frac{1}{2}$, and $\partial^α_t$ and $\partial^β_t$ denote the Caputo derivative of order $α$ and $β$ respectively. The processes $w^k_t$, $k\in \mathbb{N}=\{1,2,\cdots\}$, are independent one-dimensional Wiener processes defined on a probability space $Ω$, $L$ is a second order operator of either divergence or non-divergence type, and $Λ^k$ are linear operators of order up to two. The coefficients of the equations depend on $ω(\in Ω), t,x$ and are allowed to be discontinuous. This class of SPDEs can be used to describe random effects on transport of particles in medium with thermal memory or particles subject to sticking and trapping.