function ind = codonIndex(Obs)
% codonIndex this function returns the index of the codon Obs(1:3) in the emission
% probability matrix
% Obs is a 1x3 array

    ind = 4^2*Obs(1) + 4*Obs(2) + Obs(3);

end

