[Bioperl-l] reverse complement - and features

[Roy Chaudhuri]

> do you know of an easy way to reverse complement the sequence of an
> embl entry and the features associated with that sequence.
> I'm looking at genomic syntenic regions, and in some organisms the
> gene is in the plus strand while in others is in the minus.

Try the revcom_with_features method from Bio::SeqUtils:

$revcom=Bio::SeqUtils->revcom_with_features($seq);

(Apparently the documentation isn't correct for this. I have submitted a 
patch to Bugzilla.)

Roy.
--
Dr. Roy Chaudhuri
Bioinformatics Research Fellow
Division of Immunity and Infection
University of Birmingham, U.K.

http://xbase.bham.ac.uk