[Bioperl-l] Unusual behaviour of SeqIO::tigr
Morten Lindow
morten at binf.ku.dk
Thu Mar 4 10:57:11 EST 2004
I am trying to build a table of genomic features from the tigrxml-format
of a rice (pseudo)chromosome.
However the tigr-parser seem to behave differently from genbank/embl:
tigr.pm considers every TU/etc an individual sequence, and hence resets
its coordinate system every time it starts on a new TU.
My question is:
Is there a bioperl-way to get to the global coordinates, like when I am
parsing a genbankfile of a whole chromosome?
my $tigrin = Bio::SeqIO->new( -format => 'tigr', -file => 'chr01.xml');
while (my $seq = $tigrin->next_seq){ #I would expect this loop to run
only once
foreach my $feat ( $seq->get_SeqFeatures ) {
foreach my $location ($feat->location->each_Location){
print $location->to_FTstring, "\n";
#blablabla
}
}
}
cheers,
Morten
More information about the Bioperl-l
mailing list