function compseq = dnacomplement(seq)
%#codegen
    lseq = lower(seq);
    assert(all(lseq=='a' | lseq=='c' | lseq=='g' | lseq=='t'));
    compseq = seq;
    compseq(seq=='a') = 't';
    compseq(seq=='c') = 'g';
    compseq(seq=='g') = 'c';
    compseq(seq=='t') = 'a';
    compseq(seq=='A') = 'T';
    compseq(seq=='C') = 'G';
    compseq(seq=='G') = 'C';
    compseq(seq=='T') = 'A';
    
end