what is the difference between mixture of two normal distributions and sum of two independent variables The following denotes a mixture of a standard normal with a
normal with the same mean but 100 times the variance:
$0.95 \mathrm{~N}(0,1)+ 0.05 \mathrm{~N}(0,100)$
Let Y = 0.95 X + 0.05 Z with X,Z are independent normal distribution with X ~ N(0,1) Z~N(0,100)
I remember Y = X + Z if X ~ N(u1,s1) and Z ~ N(u2,s2), then Y would follow N(u1 + u2,s1+s2)
These two conclusions seems different and makes me confused.
Hope someone can explain

If $X_1 \sim \mathsf{Norm}(\mu_1,\sigma_1)$ and, independently,
$X_2 \sim \mathsf{Norm}(\mu_1,\sigma_2)$ then
$$X_1+X_2 \sim \mathsf{Norm}\left(\mu_1+\mu_2,\sqrt{\sigma_1^2+\sigma_2^2}\right).$$
In other words, it is the variances that add,
not the standard deviations.
Also, there are some things about mixture distributions you may not understand. Although
I try to illustrate a 50-50 mixture of two
normal distributions (with different standard
deviations) below, you might want to look at the
Wikipedia page on 'mixture distributions' and
the 'Related' links in the margin of this page
that mention 'mixture' in their titles.
[The Wikipedia page gives examples of mixture
distributions with different means, but the formulas for PDFs, means, and variances are relevant.]
Because it is easier to visualize a 50:50 mixture
graphically than a 95:5 mixture, I use a 50:50 mixture below. [In the R functions for normal
distributions, the parameters are mean and SD.]
set.seed(2022); n = 10^5
p = rbinom(n, 1, .5)     # for 50:50 mixture
x1 = rnorm(n, 0, 1)
x2 = rnorm(n, 0, 10)
y.mix = p*x1 + (1-p)*x2
y.sum = x1 + x2

For samples of size $n = 100,000,$ sample means
and standard deviations should match population means and SDs to a couple of significant digits.
summary(y.mix); var(y.mix)
     Min.   1st Qu.    Median      Mean   3rd Qu.      Max. 
-40.11645  -1.52578   0.00029   0.01102   1.53021  48.13074 
[1] 50.40555   # sample SD

summary(y.sum); var(y.sum)
     Min.   1st Qu.    Median      Mean   3rd Qu.      Max. 
-42.67832  -6.75347   0.00270  -0.00142   6.81791  48.53083 
[1] 100.6077

In the figure below the (blue) histograms nearly
approximate the (black) density curves of the
distributions for the mixture and the sum.
par(mfrow=c(2,1))
hist(y.mix, prob=T, br=100, col="skyblue2")
 curve(.5*dnorm(x,0,1)+.5*dnorm(x,0,10), add=T, lwd=3)
hist(y.sum, prob=T, br=100, col="skyblue2")
 curve(dnorm(x, 0, sqrt(1+100)), add=T, lwd=3)
par(mfrow=c(1,1))


If you have two independent normal random variables, then their sum is normal, as shown at the beginning of this Answer. However, the mixture of two independent
normal random variables need not be even nearly normal.
