http://qs321.pair.com?node_id=197793

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; }