[Bioperl-l] extract information with Bio::DB::GFF3 from gff3 file
Scott Cain
cain.cshl at gmail.com
Wed Jul 18 13:10:40 UTC 2007
Hi Xianran Li,
Assuming $gene is a Bio::DB::GFF::Feature object (there is no such thing
as Bio::DB::GFF3), then you can use the attributes method to get
anything in the ninth column:
my ($name) = $gene->attributes('Name');
The parenthesis are needed around $name because the attributes method
returns a list and the parens capture the first item of the list into
$name.
Scott
On Wed, 2007-07-18 at 13:55 +0800, Xianran Li wrote:
> Hi,
>
> I want to extract some infomation from the gff3 file like:
>
> 12001 . gene 854759 857385 . - . ID=12001.t00153;Name=receptor kinase ORK10, putative
>
> The gene position can be reterived as $gene->start, but how can I get the annotation infomatin (receptor kinase ORK10) ?
>
> Thanks for your help.
>
>
> Xianran Li
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
--
------------------------------------------------------------------------
Scott Cain, Ph. D. cain at cshl.edu
GMOD Coordinator (http://www.gmod.org/) 216-392-3087
Cold Spring Harbor Laboratory
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 189 bytes
Desc: This is a digitally signed message part
URL: <http://lists.open-bio.org/pipermail/bioperl-l/attachments/20070718/c66ec18b/attachment.sig>
More information about the Bioperl-l
mailing list