LocalAlignment.m

% Global Alignmnet
% Fixed scores
clc;
close all;
clear;

T="CTGATGAT";
S="AGCTGC";

gap=-2;
mismatch=-1;
match=1;

T = convertStringsToChars(T);
S = convertStringsToChars(S);
m = length(T);
n = length(S);

M = zeros(n+1,m+1);
M(1,:)=0*(0:m);
M(:,1)=0*(0:n);

p=max(m,n)-1;
B=zeros(p,2);
B(1,:)=[1,1];
fprintf('(%d)    ',M(1,:));
fprintf('\n')
for i=2:n+1
    s='';
    fprintf('(%d)',M(i,1));
    for j=2:m+1
        up=M(i-1,j)+gap;
        left=M(i,j-1)+gap;
        if T(j-1)==S(i-1)
            diags=M(i-1,j-1)+match;
        else
            diags=M(i-1,j-1)+mismatch;
        end
        comp=[up,left,diags,0];
        M(i,j)=max(comp);
        max_m=M(i,j);
        if max_m==up && max_m==left
            fprintf(' <|(%d) ',M(i,j));
            
        elseif max_m==up && max_m==diags
            fprintf(' \\|(%d) ',M(i,j));
            
        elseif max_m==diags && max_m==left
            fprintf(' <\\(%d) ',M(i,j));
            
        elseif max_m==up
            i_max=i-1;
            j_max=j;
            fprintf(' | (%d) ',M(i,j));
            s(j-1)='|';
        elseif max_m==left
            i_max=i;
            j_max=j-1;
            fprintf(' < (%d) ',M(i,j));
            s(j-1)=' ';
        elseif max_m==diags
            i_max=i-1;
            j_max=j-1;
            fprintf(' \\ (%d) ',M(i,j));
            s(j-1)='.';
        else
            fprintf('   (%d) ',M(i,j));
        end
        
    end
    fprintf('\n');
end
disp("M=")
disp(M);
[~, ALIGNMENT, ~] = nwalign(T,S);
disp("Alignment:")
disp(ALIGNMENT);