perl 练习：DNA反向互补

2021-05-10 本文已影响0人余绕

open FA, "$ARGV[0]";

$/=">"; #以>进行全局分割
<FA>;
while(<FA>){
chomp;
next,if($_=~/^$/); #忽略空行
my($ID,$seq)=split /\n/, $_,2; 
#print"$ID\n";
$seq=~s/\n//g;
$seq=~tr/ATGC/TACG/; #逐一替换
$seq1=reverse$seq;
print">$ID\n$seq1\n";
}