I noticed that whenever I use the function seqrcomplement, the line that calls this function takes about 80-90% of the total running time.
If I use
SRC = int2nt(5 - nt2int(fliplr(S)));
then this reduces the time to half, but it is still slow if I need to process millions of sequences. Is there a much faster alternative to seqrcomplement ?