We present an $L_{p}$-theory ($p\geq 2$) for time-fractional stochastic partial differential equations driven by Lévy processes of the type $$ \partial^α_{t}u=\sum_{i,j=1}^d a^{ij}u_{x^{i}x^{j}} +f+\sum_{k=1}^{\infty}\partial^β_{t}\int_{0}^{t} (\sum_{i=1}^dμ^{ik} u_{x^i} +g^k) dZ^k_{s} $$ given with nonzero intial data. Here $\partial^α_t$ and $\partial^β_t$ are the Caputo fractional derivatives, $α\in (0,2), β\in (0,α+1/p)$, and $\{Z^k_t:k=1,2,\cdots\}$ is a sequence of independent Lévy processes. The coefficients are random functions depending on $(t,x)$. We prove the uniqueness and existence results in Sobolev spaces, and obtain the maximal regularity of the solution.