Para poder hablar adecuadamente de la distribución de una variable aleatoria, necesitas un espacio de medida (con una medida), y específicamente un espacio de probabilidad $(\Omega, \mathscr A, \mathbb P)$, donde por supuesto $\mathbb P (\Omega) = 1$.
La distribución de una variable aleatoria $X$ es la medida de imagen $X(\mathbb P)$, es decir,
$$ X(\mathbb P)(A) :=\mathbb P(X^{-1}(A)) \quad \text{para }A\in\mathscr A,$$
En un contexto fuera de la teoría de la medida, me encontré por primera vez con esto como la 'función de distribución acumulativa' de una variable aleatoria, es decir, la función que da la probabilidad de que $X$ sea menor o igual a $x$:
$$ F_X(x) = \mathbb P (X \leq x) = \mathbb P \left(X^{-1} (-\infty, x]\right).$$
Puedes reconocer $(-\infty, x]$ como los conjuntos que generan los conjuntos de Borel en $\mathbb R$.
Una variable aleatoria tiene una distribución normal $N(0,1)$ cuando esta medida de imagen se define como:
$$ X(\mathbb P)(A) = \int_A e^{-x^2/2} \, \mathrm \lambda(\mathrm dx),$$
donde $\lambda$ es la medida de Lebesgue unidimensional. Expresado en términos menos pesados de teoría de la medida,
$$ F_X(x) = \int_{-\infty}^x e^{t^2/2} \, \mathrm dt.$$
Si $Y = X+5$, entonces tenemos $Y(\mathbb P)(A)=X(\mathbb P)(A-5)$, lo que nos da
$$ Y(\mathbb P)(A) = \int_{A-5} e^{x^2/2} \, \lambda(\mathrm dx)\\ =\int_A e^{(x-5)^2/2} \, \lambda(\mathrm dx),$$
que es la definición de $Y$ teniendo una distribución $N(5,1)$. Aquí utilizamos el hecho de que $x \mapsto x +5$ es medible (porque es continua) para justificar el último paso.