Perl reverse complement a DNA
2015-03-31 本文已影响48人
Waste_Land
sub reverse_complement {
my $dna = shift;
# reverse the DNA sequence
my $revcomp = reverse($dna);
# complement the reversed DNA sequence
$revcomp =~ tr/ACGTacgtRYKMSWrykmswBDHVbdhv/TGCAtgcaYRMKSWyrmkswVHDBvhdb/;
return $revcomp;
}