function [Delta_G_PAM] = energy_PAM(DNA_target,PAMs)
[m,n] = size(DNA_target);
Delta_G_PAM = PAMs{2}(strmatch(DNA_target(n-2:n),PAMs{1}));
end

