comseq.m

function [s] = comseq(s0)
%COMSEQ - Return complement of sequences
%
% Syntax: [s] = comseq(s0)
%
% Inputs:
%    s0     - Sequences
%
% Outputs:
%    s    - New sequences
%
% See also: REVSEQ, REVCOMSEQ

% Molecular Biology and Evolution Toolbox (MBEToolbox)
% Author: James Cai
% Email: jcai@tamu.edu
% Website: http://bioinformatics.org/mbetoolbox/
% 
% $LastChangedDate: 2013-01-05 12:04:29 -0600 (Sat, 05 Jan 2013) $
% $LastChangedRevision: 327 $
% $LastChangedBy: jcai $


if (isstruct(s0))
      if ~(isvalidaln(s0,'NUCLEOTIDE'))
	error ('ERROR: Not nucleotide sequence.')
      end
      s0=s0.seq;
end

s=s0;
a=find(s==1); c=find(s==2); g=find(s==3); t=find(s==4);
s(a)=4; s(c)=3; s(g)=2; s(t)=1;