During the course of going to the covariant formulation of the Dirac equation, we have the following:
$$ i\hbar \partial_{t} \psi = [ c(\mathbf{\alpha}\cdot\mathbf{\hat{p}}) + \beta m c^{2}]\psi $$
At this point, it seems that $\psi$ is most likely a state "ket" in the Hilbert space, since everything is basis independent. In order to move to the covariant formulation of this equation, the first step is to replace $\mathbf{\hat{p}}$ by $-i\hbar\nabla$ and multiply the above by $\beta$ from the left. Then, $$ i\hbar[\gamma^{0}\partial_{0} + \gamma^{i}\partial_{i}]\psi = i \hbar \gamma^{\mu}\partial_{\mu}\psi = m c\psi $$ And, that is it.
However, as far as I am aware (for example, in Sakurai, chapter 1), when replacing $\mathbf{\hat{p}}$ by $-i\hbar\nabla$ what we actually do is this:
$$ \langle \mathbf{x} | \mathbf{\hat{p}} | \Psi \rangle = - i\hbar \mathbf{\nabla}\langle\mathbf{x}|\Psi\rangle = - i\hbar \mathbf{\nabla}\psi(\mathbf{x})$$ That is, we explicitly choose the coordinate basis and the $\psi$ here is a wavefunction and not a state vector (which is coordinate/basis independent as a vector should be) anymore. This would imply that in the Dirac equation (covariant form), we no longer have a "vector", but this is not true either since taking the conjugate of the equation gives us $\psi^\dagger$ and not $\psi^{*}$ - in fact, it is a $4 \times 1$ column vector.
So what exactly is happening here? Are we choosing the coordinate representation or is it something else? If we are choosing the representation, how are we still left with a column vector?