[EMBOSS] water/needle interpreting swiss-format protein sequence files as DNA

[Peter Rice]

john walshaw (JIC) wrote:
> Hello,
> 
> If I create two swiss-format protein sequence files by using seqret on
> an EMBOSS protein sequence DB, then use these two files as the input to
> water or needle, I have to use "-sprotein1 -sprotein2" - otherwise they
> are treated as DNA sequences.


Hmm ... so the ID line is being read as EMBL rather than SWISSPROT/UNIPROT.

Can you send me the files and I will take a look.

regards,

Peter