[Bioperl-l] Problem using clustalw.pl

[Giuseppe Torelli]

Hi,

I'm trying to align two protein sequences in fasta format
by using clustalw.pl script in the bioperl/examples directory.
The version of bioperl is 1.0.1. This is the error:
------------- EXCEPTION  -------------
MSG: Bad input data (sequences need an id ) or less than 2 sequences in 
ARRAY(0x83b5bbc) !
STACK Bio::Tools::Run::Alignment::Clustalw::align 
/usr/lib/perl5/site_perl/5.6.1/Bio/Tools/Run/Alignment/Clustalw.pm:501
STACK main::vary_align_order bioperl-1.0.1/examples/clustalw.pl:145
STACK toplevel bioperl-1.0.1/examples/clustalw.pl:88

And this is the fasta sequence:

>A01046.PE1
MWLLLTMASLISVLGTTHGLFGKLHPGSPEVTMNISQMITYWGYPNEEYEVVTEDGYILE
VNRIPYGKKNSGNTGQRPVVFLQHGLLASATNWISNLPNNSLAFILADAGYDVWLGNSRG
NTWARRNLYYSPDSVEFWAFSFDEMAKYDLPATIDFIVKKTGQKQLHYVGHSQGTTIGFI
AFSTNPSLAKRIKTFYALAPVATVKYTKSLINKLRFVPQSLFKFIFGDKIFYPHNFFDQF
LATEVCSREMLNLLCSNALFIICGFDSKNFNTSRLDVYLSHNPAGTSVQNMFHWTQAVKS
GKFQAYDWGSPVQNRMHYDQSQPPYYNVTAMNVPIAVWNGGKDLLADPQDVGLLLPKLPN
LIYHKEIPFYNHLDFIWAMDAPQEVYNDIVSMISEDKK
>A01465.PE1
MDAMKRGLCCVLLLCGAVFVSPSQEIHARFRRGARSYQVICRDEKTQMIYQQHQSWLRPV
LRSNRVEYCWCNSGRAQCHSVPVKSCSEPRCFNGGTCQQALYFSDFVCQCPEGFAGKCCE
IDTRATCYEDQGISYRGTWSTAESGAECTNWNSSALAQKPYSGRRPDAIRLGLGNHNYCR
NPDRDSKPWCYVFKAGKYSSEFCSTPACSEGNSDCYFGNGSAYRGTHSLTESGASCLPWN
SMILIGKVYTAQNPSAQALGLGKHNYCRNPDGDAKPWCHVLKNRRLTWEYCDVPSCSTCG
LRQYSQPQFRIKGGLFADIASHPWQAAIFAKHRRSPGERFLCGGILISSCWILSAAHCFQ
ERFPPHHLTVILGRTYRVVPGEEEQKFEVEKYIVHKEFDDDTYDNDIALLQLKSDSSRCA
QESSVVRTVCLPPADLQLPDWTECELSGYGKHEALSPFYSERLKEAHVRLYPSSRCTSQH
LLNRTVTDNMLCAGDTRSGGPQANLHDACQGDSGGPLVCLNDGRMTLVGIISWGLGCGQK
DVPGVYTKVTNYLDWIRDNMRP

What is the problem ?
-- 
Giuseppe Torelli

Bioinformatic Programmer
Laboratory of Molecular Evolution
Stazione Zoologica A. Dohrn
Villa Comunale
80121 Naples - Italy
Tel.  0039 81 5833311
Fax: 0039 81 7641355