[Bioperl-l] Length of a sequence.
skirov
skirov at utk.edu
Fri Sep 8 14:35:20 UTC 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
sequence.
>I am doing wonderfully parsing all the features information.
>y'd think this would be simple.
>What's the best thing to do?
>
>
>use Bio::SeqIO;
># 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
>
>
>
>Nick Staffa
>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
>http://lists.open-bio.org/mailman/listinfo/bioperl-l
More information about the Bioperl-l
mailing list