If you assume that the vector, which spans the null space and defines the line, i.e. its direction, already has the form ##\begin{bmatrix}0\\0\\ z\end{bmatrix}##, then you made a choice considering the basis vectors, namely, that the line is the third basis vector and the other two which all together span ##\mathbb{R}^3## are perpendicular to it. There is of course a basis, which does this. But we started with another, arbitrary orthonormal basis, for which the matrix of ##A## is given by ##\begin{bmatrix}a&b&c \\ d&e&f \end{bmatrix}## and the basis vectors have the coordinates ##\begin{bmatrix}1\\0\\0\end{bmatrix},\begin{bmatrix}0\\1\\0\end{bmatrix},\begin{bmatrix}0\\0\\1\end{bmatrix}##. In this basis, we cannot know whether the null space is identical to the third basis vector - it might or might not be the case.
$$
A \cdot \vec{x} = \begin{bmatrix}a&b&c \\ d&e&f \end{bmatrix} \cdot \begin{bmatrix}x\\y\\z\end{bmatrix} = 0 \\
\Longleftrightarrow \\
\begin{cases} ax+by+cz=0\\dx+ey+fz=0 \end{cases}
$$
The two equations each define a plane in three dimensional space. If we demand both to be true for the same choices of ##x,y,z## then this means we demand to be ##(x,y,z)## on both planes, i.e. the intersecting line - the one which spans the null space, as all points ##(x,y,z)## on the line are mapped by ##A## to the zero vector. Now row manipulation as those which lead to the row echelon form are nothing else than solving this system of linear equations by multiplying them or substracting them. Let's say we achieve the equations (and row echelon form)
$$
\begin{cases} x+gz=0\\y+hz=0 \end{cases}\\
\Longleftrightarrow \\
A' \cdot \vec{x} = \begin{bmatrix}1&0&g \\ 0&1&h \end{bmatrix} \cdot \begin{bmatrix}x\\y\\z\end{bmatrix} = 0
$$
Then ##\left\{\left. \begin{bmatrix}x \\y\\ z \end{bmatrix} \in \mathbb{R}^3 \, \right| \, A'\cdot \vec{x}= 0 \right\} = \left\{\left. \begin{bmatrix}x \\y\\ z \end{bmatrix} \in \mathbb{R}^3 \, \right| \, x=-gz\; , \;y=-hz \; \right\} = \left\{\left. \begin{bmatrix}-gz \\-hz \\ z \end{bmatrix} \in \mathbb{R}^3 \, \right| \, z \in \mathbb{R} \right\}##.
In physics the free variable ##z## is often abbreviated by ##t## as the (time) parameter of a walk along the line which here spans the null space. I'd say it is helpful to choose some numbers (small integers) for ##A## and draw the two planes given by the two equations, then see where they intersect, perform the calculations for the row echelon form and see whether the null space (intersection line) is given by the last equation as stated. E.g.
$$
A = \begin{bmatrix}-3 & 0 & -6 \\ 1 &4 &-2 \end{bmatrix}
$$
In this scenario, the null space isn't automatically spanned by a vector ##\begin{bmatrix}0&0&z\end{bmatrix}##. This is only the case, if it turns out that ##g=h=0##, but usually they aren't and the null space lies somehow diagonal according to the given basis. In the above, I didn't change the basis. (In general it could only be, that we need to shuffle the naming of ##x,y,z## in order to get the row echelon form above, which means a permutation of columns.)