[Bioperl-l] Parent/parent_id attribute
Daniel Lang
Daniel.Lang at biologie.uni-freiburg.de
Mon Sep 22 18:13:57 UTC 2014
Hi,
I'm using bioperl 1.6.923-1 (Ubuntu Trusty package) and
Bio::DB::SeqFeature to store and manipulate GFF3 files.
I'm wondering why the "Parent" GFF3 attributes are stored as parent_id
values in the feature objects, but not returned as such in the gff3_string?
GFF3:
Chr01 transdecoder mRNA 5216 5627 . + .
ID=T1.Chr01.mRNA.1;Parent=T1.Chr01.gene.1;Alias=T1.asmbl_1|m.6484,T1.ORF;Name=T1.Chr01.mRNA.1
Example debugger trace after fetching stored feature:
x $f
0 Bio::DB::SeqFeature=HASH(0x3e3a798)
'attributes' => HASH(0x3e3a858)
'Alias' => ARRAY(0x3e3a8b8)
0 'T1.asmbl_1|m.6484'
1 'T1.ORF'
'load_id' => ARRAY(0x3e3aca8)
0 'T1.Chr01.mRNA.1'
'parent_id' => ARRAY(0x3e3acf0)
0 'T1.Chr01.gene.1'
'is_circular' => 0
'name' => 'T1.Chr01.mRNA.1'
'phase' => undef
'primary_id' => 2428
'ref' => 'Chr01'
'score' => undef
'source' => 'transdecoder'
'start' => 5216
'stop' => 5627
'store' => Bio::DB::SeqFeature::Store::DBI::mysql=HASH(0x39b95d0)
'class_loaded' => HASH(0x3e3a2b8)
'Bio::DB::SeqFeature' => 1
'dbh' => DBI::db=HASH(0x3dc1e40)
empty hash
'dumpdir' => '/tmp'
'is_temp' => undef
'namespace' => undef
'seqfeatureclass' => 'Bio::DB::SeqFeature'
'settings_cache' => HASH(0x3dc1d98)
'autoindex' => 1
'compress' => 0
'index_subfeatures' => 1
'serializer' => 'Storable'
'writeable' => undef
'strand' => 1
'type' => 'mRNA'
x $f->gff3_string
0
"Chr01\cItransdecoder\cImRNA\cI5216\cI5627\cI.\cI+\cI.\cIName=T1.Chr01.mRNA.1;ID=2428;Alias=T1.asmbl_1%7Cm.6484,T1.ORF"
What is the best practice to store parentage? I'm currently adding an
additional "Parent" value using add_tag_value.
Or is this a bug in the version I'm using?
Best,
Daniel
--
Dr. Daniel Lang
University of Freiburg, Plant Biotechnology
Schaenzlestr. 1, D-79104 Freiburg
fax: +49 761 203 6945
phone: +49 761 203 6989
homepage: http://www.plant-biotech.net/
http://www.cosmoss.org/
e-mail: daniel.lang at biologie.uni-freiburg.de
#################################################
My software never has bugs.
It just develops random features.
#################################################
More information about the Bioperl-l
mailing list