function genome = readfile(filename)

fid = fopen(filename);
k = 1;
while 1
    tline = fgetl(fid);
    if ~ischar(tline), break, end
    genome{k} = tline;
%     disp(tline)
    k = k+1;
end
fclose(fid);
