[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
>