[Bioperl-l] reading ace files

Robson Francisco de Souza rfsouza@citri.iq.usp.br
Fri, 24 May 2002 14:23:09 -0300 (BRT)


	Hi,

	You may use Bio::Tools::Alignment::Consed, from Chad Matsalla.

	Robson

On Fri, 24 May 2002, Susan J. Miller wrote:

> 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
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l@bioperl.org
> http://bioperl.org/mailman/listinfo/bioperl-l
>