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;
}
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;
}
本文标题:Perl reverse complement a DNA
本文链接:https://www.haomeiwen.com/subject/snkoxttx.html
网友评论