[Bioperl-l] *major* error in genbank parser or am i just insane?

[Lincoln Stein]

Hilmar Sez:
>  feature) are /not/ the same. You can't collapse this all into one 

I've been trying to stay out of this, but Hilmar's last comment is driving me 
nuts.  How do these two statements come to be semantically different?

	complement(join(2691..4571,4918..5163))
	join(complement(4918..5163),complement(2691..4571))

They both produce the same DNA sequence in the same orientation.  If GenBank 
chooses to represent numeric 0 as +0 sometimes and as -0 others, are we 
morally obligated to maintain the distinction in the data model and round 
trip it?

Lincoln

-- 
========================================================================
Lincoln D. Stein                           Cold Spring Harbor Laboratory
lstein@cshl.org			                  Cold Spring Harbor, NY
========================================================================