In general relativity when deriving the geodesic equation
$$\ddot{x}^\mu + \Gamma^\mu_{\alpha\beta}\dot{x}^\alpha\dot{x}^\beta = 0\tag{1}$$
from the action
$$S = \int d\tau \sqrt{|g_{\mu\nu} \dot{x}^\mu\dot{x}^\nu|}\tag{2}$$
(using (-+++) convention) one assumes that it is possible to choose an affine parametrization
$$-g_{\mu\nu} \dot{x}^\mu\dot{x}^\nu = 1.\tag{3}$$
This parametrization implies that the trajectory is timelike. For all other curves (e.g. spacelike curves) an affine parametrization is not possible.
Is it possible to derive that any solution of the geodesic equation (using an arbitrary parametrization) is timelike (probably not because there also exist spacelike geodesics)?
Or am I restricting to timelike curves when deriving the Euler-Lagrange equation out of the action? This is also strange because then I cannot allow for arbitrary $\delta x^\mu$.
Or am I just assuming that the curve is timelike after I derived the geodesic equation, then choosing the affine parametrization and then simplifying the geodesic equation to the form it is usually stated? But then why can I insert this into the action before deriving the Euler-Lagrange equation?
Or is there some other solution?