Beefy Boxes and Bandwidth Generously Provided by pair Networks
We don't bite newbies here... much
 
PerlMonks  

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?
Username:
Password:

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

How do I use this? | Other CB clients
Other Users?
Others surveying the Monastery: (4)
As of 2014-11-23 06:24 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    My preferred Perl binaries come from:














    Results (128 votes), past polls