Use findpeaks with two outputs:
[pks,locs] = findpeaks(SSignal);
The ‘locs’ variable will have the index values of the peaks it returns in the ‘pks’ variable.
To plot them, use ‘locs’ as a subscript in the vectors in your second plot call:
plot(time(locs), SSignal(locs),'o')