[Bioperl-l] Genbank/EMBL output

Alex Hudek alex@genet.sickkids.on.ca
Fri Jan 17 20:35:02 EST 2003


At first glance I thought this as well. Apparently the complement operator
is complement only rather than reverse complement as one might think. So
you do not have to reverse the ranges when moving the complment tag inside
the join. 

Alex

On Fri, 17 Jan 2003, Matthew Pocock wrote:

> Hi,
> 
> Are these the same? I would have thought that:
> 
> complement(join(1..50,60..100))
> 
> is equivalent to:
> 
> join(complement(60..100),complement(1..50))
> 
> or do the operators join & complement not work how I thought?
> 
> Matthew
> 
> 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
> > 
> 
> 
> -- 
> BioJava Consulting LTD - Support and training for BioJava
> http://www.biojava.co.uk
> 


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





More information about the Bioperl-l mailing list