Beefy Boxes and Bandwidth Generously Provided by pair Networks
Pathologically Eclectic Rubbish Lister
 
PerlMonks  

DNA Reverse Complement

by mousey (Scribe)
on Sep 14, 2002 at 05:02 UTC ( #197793=snippet: print w/replies, xml ) Need Help??
Description: A function to figure out the reverse complement of a DNA sequence
sub revdnacomp {
  my $dna = @_;
  my $revcomp = reverse($dna);

  $revcomp =~ tr/ACGTacgt/TGCAtgca/;

  return $revcomp;
}
Replies are listed 'Best First'.
Re: DNA Reverse Complement
by jkahn (Friar) on Sep 14, 2002 at 05:12 UTC
    I think you'll have a problem with the argument structure:
    sub revdnacomp { # my $dna = @_; # the above means $dna gets the number of # arguments in @_, since it's a scalar context! my $dna = shift; # or my $dna = shift @_; # ah, scalar context of scalar gives expected results. # my ($dna) = @_; # would work, too my $revcomp = reverse($dna); $revcomp =~ tr/ACGTacgt/TGCAtgca/; return $revcomp; }
    but why not allow the complementation over more than one element?
    sub revcompl { # operates on all elements passed in my (@dna) = @_; foreach my $segment (@dna) { my $revcomp = reverse($segment); $revcomp =~ tr/ACGTacgt/TGCAtgca/; push @done, $revcomp; } return @done; # or reverse @done; # what's best semantics? }
    Did I miss anything?
      sub revcompl { return map { (my $rev = reverse $_) =~ tr/ACGTacgt/TGCAtgca/; $rev + } @_ }

      Makeshifts last the longest.

        The above sub does not work
Log In?
Username:
Password:

What's my password?
Create A New User
Node Status?
node history
Node Type: snippet [id://197793]
help
Chatterbox?
[tye]: I don't know if you can just dismount, but I thought so.
[davido]: yeah, umount -f isn't powerful enough.
[tye]: probably something in the init subsystem that does the mounting that you could disable and reboot.
davido needs to close laptop to board flight home from yapc.
[davido]: I'll look into it in a few hours probably.
[davido]: when i do get to that point I think I'll do it in a vm snapshot just in case. :)
[oiskuu]: tye, you were right: loginuid/sessionid are part of task struct if compiled with AUDITSYSCALL. I have some doubts if you should actually depend on that feature.
[Corion]: oiskuu: Depends on what you want to do with that information
[tye]: I'm not depending on that feature. But I could in this environment. I'm using getlogin(). shrug
[Corion]: For benign logging (which user started this DB instance), it's OK

How do I use this? | Other CB clients
Other Users?
Others taking refuge in the Monastery: (5)
As of 2017-06-23 20:06 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?
    How many monitors do you use while coding?















    Results (555 votes). Check out past polls.