The question has been resolved. The following code helps to change color of each bar in the grouped bar plot.
c = [25;55];
names = {'GR = 5.57' ; 'GR = 6.17'};
y = [604 603 ; 903 903 ; 1083 1083]';
h = bar(c,y,'FaceColor','flat');
h(1).CData(1,:) = [0 1 0]; % group 1 1st bar
h(1).CData(2,:) = [0 1 0]; % group 1 2nd bar
h(2).CData(1,:) = [0 1 0]; % group 2 1st bar
h(2).CData(2,:) = [0 1 0]; % group 2 2nd bar
h(3).CData(1,:) = [1 0 0]; % group 3 1st bar
h(3).CData(2,:) = [1 0 0]; % group 3 2nd bar