I understand that to transform any normal distribution to the standard normal distribution we use the Z-transformation: $$Z = \frac{x - \mu}{\sigma}$$
What I am looking for is a basic and intuitive explanation of how does this transformation work. i.e how does finding the Z-score of a value of a normal distribution means finding its equivalent in the standard normal distribution?
$$f_X(x)=\frac{1}{\sigma\sqrt{2\pi}}e^{-\frac{1}{2\sigma^2}(x-\mu)^2}$$
with $x \in \mathbb{R}$, $\mu \in \mathbb{R}$, $\sigma >0$
Let's take the new random variable
$$Z=\frac{X-\mu}{\sigma}$$
Let's calculate the CDF of Z by definition:
$$F_Z(z)=\mathbb{P}[Z \leq z]=\mathbb{P}[\frac{X-\mu}{\sigma} \leq z]=\mathbb{P}[X \leq \sigma z+\mu]=F_X(\sigma z+\mu)=\int_{-\infty}^{\sigma z+\mu}f_X(x)dx$$
to get $f_Z(z)$ simply derive the CDF
$$f_Z(z)=\frac{d}{dz}F_Z(z)=f_X(\sigma z+\mu)\cdot \sigma=\frac{1}{\sqrt{2\pi}}e^{-\frac{z^2}{2}}$$
$$g(x|\mu,\sigma)=\frac{1}{\sigma}\psi(\frac{x-\mu}{\sigma})$$
where $\psi(.)$ is the Standard Normal Density