sub revcom
{
# A subroutine to compute the reverse complement of DNA sequence
# 一个获取DNA互补序列的子程序
my($dna)=@_;
print $dna."\n";
my ($revcom)=reverse($dna);
$revcom=~tr/ACGTacgt/TGCAtgca/;
return $revcom;
}
my $plus_dna="ATATTATATATAGCGC";
my $the_result;
$the_result=revcom($plus_dna);
print $the_result;