Explicit bounds are given for the Kolmogorov andWasserstein distances between a mixture of normal distributions, by which we mean that the conditional distribution given some sigma-algebra is normal, and a normal distribution with properly chosen parameter values. The bounds depend only on the first two moments of the first two conditional moments given the sigma-algebra. The proof is based on Steins method. As an application, we consider the Yule-Ornstein-Uhlenbeck model, used in the field of phylogenetic comparative methods. We obtain bounds for both distances between the distribution of the average value of a phenotypic trait over n related species, and a normal distribution. The bounds imply and extend earlier limit theorems by Bartoszek and Sagitov.