perl应用:DNA互补序列的获取



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;


你可能感兴趣的:(perl应用:DNA互补序列的获取)