Beefy Boxes and Bandwidth Generously Provided by pair Networks
XP is just a number

Re: Another Reverse Complement Question...

by Cristoforo (Deacon)
on Jun 04, 2012 at 20:31 UTC ( #974373=note: print w/ replies, xml ) Need Help??

in reply to Another Reverse Complement Question...

I don't work with BioPerl, but they do have a reverse complement method, revcom and it can be used like this example:

#!/usr/bin/perl use strict; use warnings; use 5.014; use Bio::PrimarySeq; my $seqobj = Bio::PrimarySeq->new(-seq => 'ACTGTGGCGTCAACTG', -alphabet => 'dna', -id => 'test'); say $seqobj->seq; say $seqobj->revcom->seq; __END__ *** print out ACTGTGGCGTCAACTG CAGTTGACGCCACAGT

But this doesn't do the translation you want: A => T, C => G, G => C and T => A. So, I don't know if this is any help.

Update: After examining the reversed complement sequence, I see indeed it performs the complement as you stated but it also reverses the sequence with the first letter in the original being complemented and last in the reversed sequence, the second in the original complemented and being next to last in the reversed complement sequence, etc.

Comment on Re: Another Reverse Complement Question...
Download Code

Log In?

What's my password?
Create A New User
Node Status?
node history
Node Type: note [id://974373]
and the web crawler heard nothing...

How do I use this? | Other CB clients
Other Users?
Others scrutinizing the Monastery: (9)
As of 2015-11-25 14:21 GMT
Find Nodes?
    Voting Booth?

    What would be the most significant thing to happen if a rope (or wire) tied the Earth and the Moon together?

    Results (677 votes), past polls