function NTGL = Nontranslated(AGL,GGL,GMGL,ORFGL)
clc, close all;

% Find all the entries in ORFGL that are not in AGL:
% In other words: (ORFGL - AGL)
for i = 1:length(ORFGL)
    i
    for j = 1:length(AGL)
        if (ORFGL{i,2} == AGL{j,2})&(ORFGL{i,3} == AGL{j,3}), 
            x{i,1} = ORFGL{i,1};
            x{i,2} = ORFGL{i,2};
            x{i,3} = ORFGL{i,3};
            break;
        end
    end
end

k = 1;
for i = 1:length(x)
    if isempty(x{i,1}),
        NTORFGL{k,1} = k;  % Nontranslated ORF Gene Locations
        NTORFGL{k,2} = ORFGL{i,2};
        NTORFGL{k,3} = ORFGL{i,3};
        NTORFGL{k,4} = ORFGL{i,4};
        NTORFGL{k,5} = ORFGL{i,5};
        k = k + 1;
    end
end

clear x;
% Find all the entries in NTORFGL that are not in GGL:
% In other words: (ORFGL - AGL - GGL)
for i = 1:length(NTORFGL)
    i
    for j = 1:length(GGL)
        if (NTORFGL{i,2} == GGL{j,2})&(NTORFGL{i,3} == GGL{j,3}), 
            x{i,1} = NTORFGL{i,1};
            x{i,2} = NTORFGL{i,2};
            x{i,3} = NTORFGL{i,3};
            break;
        end
    end
end

k = 1;
for i = 1:length(x)
    if isempty(x{i,1}),
        NTORFGGL{k,1} = k;
        NTORFGGL{k,2} = NTORFGL{i,2};
        NTORFGGL{k,3} = NTORFGL{i,3};
        NTORFGGL{k,4} = NTORFGL{i,4};
        NTORFGGL{k,5} = NTORFGL{i,5};
        k = k + 1;
    end
end

clear x;
% Find all the entries in NTORFGGL that are not in GMGL:
% In other words: (ORFGL - AGL - GGL - GMGL)
for i = 1:length(NTORFGGL)
    i
    for j = 1:length(GMGL)
        if (NTORFGGL{i,2} == GMGL{j,2})&(NTORFGGL{i,3} == GMGL{j,3}), 
            x{i,1} = NTORFGGL{i,1};
            x{i,2} = NTORFGGL{i,2};
            x{i,3} = NTORFGGL{i,3};
            break;
        end
    end
end

k = 1;
for i = 1:length(x)
    if isempty(x{i,1}),
        NTGL{k,1} = k;
        NTGL{k,2} = NTORFGGL{i,2};
        NTGL{k,3} = NTORFGGL{i,3};
        NTGL{k,4} = NTORFGGL{i,4};
        NTGL{k,5} = NTORFGGL{i,5};
        k = k + 1;
    end
end