correspondence of the fitgmdist output to initial mean components

Hi,

Imagine that we want to fit mixtures of Gaussians to a set of 2D points. I already have a guess of the potential mean components. Therefore, I supplied fitgmdist with my initial guess using the following code:

S.mu = guess; % guess is a nx2 matrix
S.Sigma = repmat([2 0;0 2],1,1,size(S.mu,1));
S.ComponentProportion = 1/size(S.mu,1)*ones(1,size(S.mu,1)); % equal proportion

GMModel = fitgmdist(data, size(S.mu,1), ...
                             'Options',statset('Display','off','MaxIter',1000,'TolFun',1e-6), ...
                             'Start',S);

The question:

Does GMModel.mu have the same order as S.mu? In other words, does GMModel.mu(1,:) corresponds to the initial S.mu(1,:)?