[Bioperl-l] question about Bio::DB::GFF
xianran li
xianranli78 at yahoo.com.cn
Tue Sep 4 02:11:09 UTC 2007
Hi,
I tried to load the gff3 file with load_gff.pl and extrac some information with Bio::DB::GFF. Althougth this code work properly under windows xp, the $seg got nothing when i run it under Linux.
Here is my code and the gff3 file,
####################################################################
#!/usr/local/bin/perl -w
use strict;
use Bio::SeqIO;
use Bio::DB::GFF;
my $in_gff = Bio::DB::GFF->new( -adaptor => 'dbi::mysqlopt',
-dsn => 'dbi:mysql:test',
-aggregator => ['coding'],
-user => "lixr",
-pass => "123456"
);
my $seg = $in_gff->segment'BGIOSIBCE000001.1');
print $seg->abs_start."\n";
##################################################################
##gff-version 3
##sequence-region Chr01 1 43037
Chr01 bgf mRNA 18113 20165 . + . ID=BGIOSIBCE000001.1
Chr01 bgf CDS 18113 19150 . + 0 Parent=BGIOSIBCE000001.1
Chr01 bgf CDS 19344 20165 . + 0 Parent=BGIOSIBCE000001.1
Chr01 bgf mRNA 30220 36442 . + . ID=BGIOSIBCE000002.1
Chr01 bgf CDS 30220 30387 . + 0 Parent=BGIOSIBCE000002.1
Chr01 bgf CDS 31128 31226 . + 0 Parent=BGIOSIBCE000002.1
Chr01 bgf CDS 32228 32331 . + 0 Parent=BGIOSIBCE000002.1
Chr01 bgf CDS 33907 34715 . + 1 Parent=BGIOSIBCE000002.1
Chr01 bgf CDS 34799 34921 . + 2 Parent=BGIOSIBCE000002.1
Chr01 bgf CDS 35003 35091 . + 2 Parent=BGIOSIBCE000002.1
Chr01 bgf CDS 35179 35379 . + 0 Parent=BGIOSIBCE000002.1
Chr01 bgf CDS 35981 36442 . + 0 Parent=BGIOSIBCE000002.1
Chr01 bgf mRNA 38143 39015 . - . ID=BGIOSIBCE000003.1
Chr01 bgf CDS 38143 38541 . - 0 Parent=BGIOSIBCE000003.1
Chr01 bgf CDS 38649 38813 . - 0 Parent=BGIOSIBCE000003.1
Chr01 bgf CDS 38917 39015 . - 0 Parent=BGIOSIBCE000003.1
Chr01 bgf mRNA 39545 42080 . + . ID=BGIOSIBCE000004.1
Chr01 bgf CDS 39545 40584 . + 0 Parent=BGIOSIBCE000004.1
Chr01 bgf CDS 40677 41042 . + 1 Parent=BGIOSIBCE000004.1
Chr01 bgf CDS 41130 41208 . + 1 Parent=BGIOSIBCE000004.1
Chr01 bgf CDS 41740 41920 . + 0 Parent=BGIOSIBCE000004.1
Chr01 bgf CDS 42037 42080 . + 2 Parent=BGIOSIBCE000004.1
#################################################################
I would appreaciate if any one can give me some clues/link to accomplish this.
thanks in advance ,
Xianran Li
---------------------------------
雅虎邮箱,以安全著称,是值得信赖的邮箱专家!
More information about the Bioperl-l
mailing list