[Bioperl-l] Why does Bio::DB::GFF::Feature::gff3_string swap start and stop coordinates??
Mark Johnson
johnsonm at gmail.com
Wed May 16 19:11:18 UTC 2007
On 5/8/07, Chris Fields <cjfields at uiuc.edu> wrote:
> I believe all seqfeature location coordinates are designed to have
> start < stop for consistency; in cases where the strand matters (CDS,
> gene, etc.) then the strand is set to 1 or -1. When start > stop,
> the two are reversed and the strand is flipped; at least that's the
> way locations are set up in BioPerl.
>
> chris
Oh yeah? I always tend to ensure that (start < stop), regardless
of strand, when working with sequence features...the other day, I
caught Glimmer2 emitting a prediction on the plus strand with start >
stop. I was going to work up a patch for the parser, but I wonder,
should I just force everything to start < stop? Or only predictions
on the plus strand? Should all the parsers for all the ab initio
predictors ensure they emit features with coordinates like this?
More information about the Bioperl-l
mailing list