function [fp, norm_fp] = freq_profile(seq, pat)
len = length(pat);
fp = zeros(1, len);
seq_uc = upper(seq);
for i = 1:len 
    hits = strfind(seq_uc, upper(string(pat{i})));
    fp(i) = length(hits);
end

% Normalize the profile
tot = sum(fp);
norm_fp = fp / tot;