Let $Z_{1}$ and $Z_{2}$ be two independent normally distributed random variables with expectations $\mu_{1},\mu_{2}\in\mathbb{R}$ and variances $\sigma_{1}^2,\sigma_{2}^2\in (0,\infty)$ .
I would like to prove that the regular conditional distribution $\mathbb{P}(Z_{1}\in \cdot$ $|Z_{1}+Z_{2}=x)$, for almost every $x\in\mathbb{R}$, is equal to the normal distribution with expectation $\mu_{1}+\frac{\sigma_1^2}{\sigma_1^2+\sigma_2^2}(x-\mu_1-\mu_2)$ and variance $\frac{\sigma_2^2\sigma_1^2}{\sigma_1^2+\sigma_2^2}$.
I don't know how to start. Can anybody help, please?