[Bioperl-l] reading ace files

Susan J. Miller sjmiller@u.arizona.edu
Fri, 24 May 2002 09:39:15 -0700


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