[Bioperl-l] Length of a sequence.
skirov at utk.edu
Fri Sep 8 10:35:20 EDT 2006
Surprisingly the method is length:
perldoc -f length Bio::Seq
>===== Original Message From "Staffa, Nick (NIH/NIEHS) [C]"
<staffa at niehs.nih.gov> =====
>Wondering why I'm having trouble finding out how to get the total length of a
>I am doing wonderfully parsing all the features information.
>y'd think this would be simple.
>What's the best thing to do?
># Get the sequence
>$seqio_object = Bio::SeqIO->new(-file => "$filename" );
>$seq_object = $seqio_object->next_seq;
>>> What's the magic code to put HERE????<<<<<<<<<<<<<<<<<<<
>print "Sequence $fileroot has length $len\n";
># Get the features as an array of objects
># This seems to get them all. Makes an array of objects.
>my @features = $seq_object->get_SeqFeatures(); # just top level
>Telephone: 919-316-4569 (NIEHS: 6-4569)
>Scientific Computing Support Group
>NIEHS Information Technology Support Services Contract
>(Science Task Monitor: Jack L. Field( field1 at niehs.nih.gov )
>National Institute of Environmental Health Sciences
>National Institutes of Health
>Research Triangle Park, North Carolina
>Bioperl-l mailing list
>Bioperl-l at lists.open-bio.org
More information about the Bioperl-l