Start with a matrix whose columns are the vectors you have.
Then reduce this matrix to row-echelon form.
A basis for the columnspace of the original matrix is given by the columns in the original matrix that correspond to the pivots in the row-echelon form.
What you are doing does not really make sense because elementary row operations do not preserve the column space, but you are looking for a basis of the column space. In the end, what you are finding is a basis for the nullspace of the matrix, but what you are looking for is a basis for the columnspace.
Now, the nullspace and the rowspace are orthogonal to each other (in fact, orthogonal complements), so you could use your two vectors to find a basis for the column space, but it's too much extra work.
So: starting with the matrix of vectors you have,
$\begin{align*} \left(\begin{array}{rrrr} 2 & 7 & 1 & 3\\ -4 & -2 & 2 & -2\\ -1 & 7 & 3 & 2\\ -2 & 2 & 2 & 0 \end{array}\right) &\to \left(\begin{array}{rrrr} -1 & 7 & 3 & 2\\ -4 & -2 & 2 & -2\\ 2 & 7 & 1 & 3\\ -2 & 2 & 2 & 0 \end{array}\right) \to \left(\begin{array}{rrrr} -1 & 7 & 3 & 2\\ 0 & -30 & -10 & -10\\ 0 & 21 & 7 & 7\\ 0 & -12 & -4 & -4 \end{array}\right)\\ &\to \left(\begin{array}{rrrr} -1 & 7 & 3 & 2\\ 0 & 3 & 1 & 1\\ 0 & 3 & 1 & 1\\ 0 & 3 & 1 & 1 \end{array}\right) \to \left(\begin{array}{rrrr} \color{red}{-1} & 7 & 3 & 2\\ 0 & \color{red}3 & 1 & 1\\ 0 & 0 & 0 & 0\\ 0 & 0 & 0 & 0 \end{array}\right) \end{align*}$ (Same thing you had)
Now, since the pivots of the row echelon form are in columns 1 and 2 (marked in red), that means that columns 1 and 2 of the original matrix are independent and span the column space. That is, you can take $\left(\begin{array}{r} 2 \\-4\\-1\\-2\end{array}\right)\quad\text{and}\quad \left(\begin{array}{r} 7\\-2 \\ 7 \\ 2 \end{array}\right).$ It should be easy to see that the two vectors are linearly independent. Also, $\left(\begin{array}{r} 1\\2\\3\\2\end{array}\right) = -\frac{2}{3}\left(\begin{array}{r} 2\\-4\\-1\\-2\end{array}\right) + \frac{1}{3}\left(\begin{array}{r}7\\-2\\7\\2\end{array}\right)$ and $\left(\begin{array}{r} 3\\-2\\2\\0\end{array}\right) = \frac{1}{3}\left(\begin{array}{r} 2\\-4\\-1\\-2\end{array}\right) + \frac{1}{3}\left(\begin{array}{r}7\\-2\\7\\2\end{array}\right)$ so the first two columns span the same space as all four.
The reason this works is that any linear combination between the columns is respected by elementary row operations (equivalently: elementary row operations do not change the nullspace). Since the columns that contain the pivots are necessarily a basis for the columnspace of the row-echelon form, the corresponding columns in the original form a basis. In fact, you can find the expressions I found above to show the other two columns are linear combinations by looking in the reduced-row echelon form instead of the original matrix; this is often easier.
This method extracts a subset of your original set that is a basis for the span of the set; that is, it gives you a way to pare your set down to a basis. If all you need is to find a basis, you can also begin by putting the vectors as the rows of a matrix and then performing elementary row operations until you get a matrix in row-echelon form. The (transpose of the) rows that contain the pivots in the row-echelon form give you a basis for the span of the original vectors. The reason this works is that elementary row operations don't change the row space, and the rows that contain the pivots are clearly a basis for the row space of the row-echelon form of the matrix. But you would not look for the nullspace in any case.