function cs = count_unique_nucleotides(seq)
%#codegen
    cs = [ sum(seq=='A'), sum(seq=='C'), sum(seq=='G'),sum(seq=='T')];
end