It is true that the light coming from a laser is not precisely a coherent state. One can measure the photon statistics to see that it only approximates the Poisson statistics. However, the OP is not concerned with the precise modeling of the light coming form a laser. Instead, the question is how a coherent state can come about as a result of the stimulated emission that occurs in a laser. To address this aspect, I present here a simplistic view of the process. It neglects the possibility that the excited atom remains excited and not radiate.
When the laser is switch on, it starts with the creation of a population inversion. Then one of the excited atoms decays spontaneously. The photon that is produced by the spontaneous decay stimulates other atoms to decay producing more photons. However, one can also see the subsequent stimulated decays together with the initial spontaneous decay as multiple spontaneous decays together with a normalization process.
Each spontaneous decay will effectively produce a superposition
$$ |\psi_{spon}\rangle = |\text{vac}\rangle\beta + |1\rangle\zeta , $$
where $|\text{vac}\rangle$ is the vacuum state, $|1\rangle$ is a single photon state with the correct cavity mode, and $\beta$ and $\zeta$ are coefficients. The reason for the vacuum state is not because the atom did not radiated, which I exclude here, but because some of the photons are produced in modes that would not survive in the long run. These I remove and replace by the vacuum state.
For $n$ such spontaneous decays, one gets
$$ |\psi_{n-spon}\rangle = \mathcal{N}\left(|\text{vac}\rangle\beta + |1\rangle\zeta\right)^n = \mathcal{N}\sum_{p=0}^n \frac{n!}{p! (n-p)!}
\left(|1\rangle\right)^p \frac{\zeta^p}{\beta^p}, $$
where $\mathcal{N}$ is a normalization constant. In the last expression I absorbed $\beta^n$ in $\mathcal{N}$ and removed the tensor products with the vacuum states. Note that a tensor product of $p$ single photon states becomes a $p$-photon Fock state
$$\left(|1\rangle\right)^p= \sqrt{p!} |p\rangle . $$
In the limit $n\rightarrow\infty$, we have
$$ \frac{n!}{(n-p)!}\rightarrow 1 . $$
So, the state becomes
$$ |\psi\rangle = \mathcal{N}\sum_{p=0}^{\infty} \frac{\alpha^p}{\sqrt{p!}} |p\rangle , $$
where we defined
$$ \frac{\zeta}{\beta} \equiv \alpha . $$
What remains is to compute the normalization constant of the state, which would then lead to the well-known expression for the coherent state.