[Bioperl-l] reverse complement - and features

Roy Chaudhuri roy at colibase.bham.ac.uk
Thu Nov 16 08:18:55 EST 2006


> 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


More information about the Bioperl-l mailing list