If you need a point forecast of $X_{t+1}$ (which you indicate in the body of your question), you should build a conditional mean model for $X_t$.
If you need a range forecast of $X_{t+1}$ (which you indicate in the title of your question), you would need a predicted density or predicted quantiles of $X_{t+1}$.
A GARCH model involves a conditional mean and a conditional variance specification, e.g.
$$
\begin{aligned}
X_t &= \sigma_t \varepsilon_t \\
\sigma_t^2 &= \omega + \alpha_1 X_{t-1}^2 + \beta_1 \sigma_{t-1}^2
\end{aligned}
$$
with a density specification for $\varepsilon_t$, e.g. $\varepsilon_t \sim i.i.d.(0,1,\theta)$ for some density function $d$ with zero mean, unit variance and possibly extra parameters $\theta$.
Given the estimated model, you can obtain the predicted density of $X_{t+1}$, $d(0,\hat\sigma^2_{t+1},\hat\theta)$, where hats denote fitted values and
$$ \hat\sigma_{t+1}^2 = \hat\omega + \hat\alpha_1 X_{t}^2 + \hat\beta_1 \hat\sigma_{t}^2. $$
You can obtain an $x\%$ confidence interval for $X_{t+1}$ using the quantiles of the predicted density of $X_{t+1}$. This will be the likely range of $X_{t+1}$, where "likely" is quantified as $x\%$ confidence.