  We consider fractional parabolic equations with variable coefficients and
establish maximal $L_{q}$-regularity in Bessel potential spaces of arbitrary
nonnegative order. As an application, we show higher order regularity and
instantaneous smoothing for the fractional porous medium equation and for a
nonlocal Kirchhoff equation.
