We consider finite element error approximations of the steady incompressible Navier-Stokes equations defined on a randomly perturbed domain, the perturbation being small. Introducing a random mapping, these equations are transformed into PDEs on a fixed reference domain with random coefficients. Under suitable assumptions on the random mapping and the input data, in particular the so-called small data assumption, we prove the well-posedness of the problem. We assume then that the mapping depends affinely on $L$ independent random variables and adopt a perturbation approach expanding the solution with respect to a small parameter $\varepsilon$ that controls the amount of randomness in the problem. We perform an a posteriori error analysis for the first order approximation error, namely the error between the exact (random) solution and the finite element approximation of the first term in the expansion with respect to ". Numerical results are given to illustrate the theoretical results and the effectiveness of the error indicator.