Bioperl: start, end, strand

Steve Chervitz sac@neomorphic.com (Steve A. Chervitz)
Sat, 17 Jun 2000 18:15:22 -0700 (PDT)


Ewan Birney writes:
 > On Tue, 13 Jun 2000, Bradley Marshall wrote:
 > 
 > > 
 > > In game, the minus strand is handled implicitly if the
 > > start is larger than the end.  Should I maintain this
 > > convention in bioperl or explicitly switch the coord's
 > > and set the strand to -1?
 > 
 > 
 > No... 
 > 
 > Bioperl has a strong convention of
 > 
 > 	start is *strictly* less than end
 > 
 > and
 > 
 > 	strand is one of 1,-1,0 (0 being not strand specific, eg, low
 > complexity repeat or a variation).
 > 
 > 
 > Although start > end => reverse strand is sometimes useful you both
 > 
 > 	(a) have difficult in expressing strand 0 type features
 > 
 > 	(b) end up with just as many times that you need to reverse
 > the coordinates as in the case where start is strictly less than end.
 > 

Some other reasons why we can't rely on start vs. end comparisons to
determine the strand:

1) Single-residue features, (start = end).

2) Circular sequences, (the last base is 5' of the first).

Steve

 > 
 > > 
 > > Brad
 > > 
 > > __________________________________________________
 > > Do You Yahoo!?
 > > Yahoo! Photos -- now, 100 FREE prints!
 > > http://photos.yahoo.com
 > > =========== Bioperl Project Mailing List Message Footer =======
 > > Project URL: http://bio.perl.org/
 > > For info about how to (un)subscribe, where messages are archived, etc:
 > > http://www.techfak.uni-bielefeld.de/bcd/Perl/Bio/vsns-bcd-perl.html
 > > ====================================================================
 > > 
 > 
 > -----------------------------------------------------------------
 > Ewan Birney. Mobile: +44 (0)7970 151230, Work: +44 1223 494420
 > <birney@ebi.ac.uk>. 
 > -----------------------------------------------------------------
 > 
 > =========== Bioperl Project Mailing List Message Footer =======
 > Project URL: http://bio.perl.org/
 > For info about how to (un)subscribe, where messages are archived, etc:
 > http://www.techfak.uni-bielefeld.de/bcd/Perl/Bio/vsns-bcd-perl.html
 > ====================================================================
 > 
=========== Bioperl Project Mailing List Message Footer =======
Project URL: http://bio.perl.org/
For info about how to (un)subscribe, where messages are archived, etc:
http://www.techfak.uni-bielefeld.de/bcd/Perl/Bio/vsns-bcd-perl.html
====================================================================