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

上一篇 下一篇

猜你喜欢

热点阅读