function genes = ffn2genes(filename);
lines = readfile(filename);
k = 1;
for i = 1:length(lines)
    if lines{i}(1) == '>',
        GL(k) = i;  % Gene Location, i.e. location of '>'
        k = k + 1;
    end
end

for i = 1:length(GL)
    if i ~= length(GL)
        GLs = (char({lines{GL(i)+1:GL(i+1)-1}}))';  % GLs = Gene Lines
        genes{i} = {strtrim(reshape(GLs,1,size(GLs,1)*size(GLs,2)))};
    else
        GLs = (char({lines{GL(i)+1:length(lines)}}))';  % GLs = Gene Lines
        genes{i} = {strtrim(reshape(GLs,1,size(GLs,1)*size(GLs,2)))};
    end
end
