On Sat, 5 Jun 2004, fsms at users.sourceforge.net wrote: > >>> s = Seq('acgt', IUPAC.IUPACAmbiguousDNA) You should be using 'ACGT' in caps. HTH, -- Michael Hoffman <hoffman at ebi.ac.uk> European Bioinformatics Institute