[Bioperl-l] reading ace files

[Susan J. Miller]

Can bioperl read ace files such as those produced by phrap?  I tried
this using the SeqIO module but next_seq and next_primary_seq both
return null on my file:

 $in = Bio::SeqIO->new('-file' => $fil, '-format' => 'ace', 'dna' =>
'dna');
 while ($seq = $in->next_seq()) {
	print "seq $seq\n";
	print "SeqID $seq->id $seq->subseq(1,10) \n";
 }


-- 
Thanks,
-susan

Susan J. Miller
Biotechnology Computing Facility
Arizona Research Laboratories
Bio West 228
University of Arizona
Tucson, AZ  85721
(520) 626-2597