Attempting to alter the color of connected components in a label matrix that are under a threshold pixel area.

Hello,

I am attempting to take this image and color the circles that are smaller than 1200 pixels cyan:

Here's what I've done so far.

%Assuming image already read in as I

[L, num] = bwlabel(I); %return label matrix and num of connected components of I
S = regionprops(I)
p = [S.Area] %create vector consisting of areas of each connected component in I

for i=1:num 
if(p(i) < 1200) %finding Areas in the vector that are less than threshold value
%not sure how to proceed from here
end
end

Any guidance would be appreciated.