[Bioperl-l] Genbank/EMBL output
Hilmar Lapp
hlapp@gnf.org
Mon, 23 Dec 2002 14:34:37 -0800
It's possible but would require a number of hacks to circumvent the
semantics of the object model. Those semantics define that a split
location does not carry strand information. It's the individual
locations that hold the strand instead. What you're asking for is logic
that checks whether all sub-locations are on the negative strand, and
if so promote the complement() operator to the join.
Why do you need the other form?
-hilmar
On Monday, December 23, 2002, at 02:04 PM, Alex Hudek wrote:
>
> Is it possible to have bioperl output features in Genbank format of the
> form
>
> complement(join(1..50,60..100))
>
> rather than
>
> join(complement(1..50),complement(60..100)) ?
>
> I understand the former form is more common/accepted?
>
>
> Alexander K. Hudek
> The Hospital for Sick Children
> Room 9107, The Centre for Applied Genomics
> 555 University Avenue
> Toronto, Ontario, M5G 1X8, Canada
> Lab phone: 416-813-7638
> Fax: 416-813-8319
> E-mail: alex@genet.sickkids.on.ca
>
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l@bioperl.org
> http://bioperl.org/mailman/listinfo/bioperl-l
>
--
-------------------------------------------------------------
Hilmar Lapp email: lapp at gnf.org
GNF, San Diego, Ca. 92121 phone: +1-858-812-1757
-------------------------------------------------------------