[Bioperl-l] Retrieve mRNA from Genome

Staffa, Nick (NIH/NIEHS) staffa at niehs.nih.gov
Mon Apr 9 15:43:19 UTC 2007


I have been retrieving sub-sequence from Genbank genomic records by use of
Bio::SeqIO
and ->get_SeqFeatures, ->start ->end ,
but now I'm looking for a quick way to extract CDS or mRNA from
a multi-segmented annotation, e.g.
     mRNA          
join(72458..72791,84573..84613,93279..94419,94481..94656,
                     94719..94992,95056..95350,95438..95553,95614..96056)

Is there such a method?
Please point me to appropriate documentation.


 
Nick Staffa 
Telephone: 919-316-4569  (NIEHS: 6-4569)
Scientific Computing Support Group
NIEHS Information Technology Support Services Contract
(Science Task Monitor: John D. Grovenstein (grovens1 at niehs.nih.gov)
National Institute of Environmental Health Sciences
National Institutes of Health
Research Triangle Park, North Carolina






More information about the Bioperl-l mailing list