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";
}
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";
}
本文标题:perl 练习:DNA反向互补
本文链接:https://www.haomeiwen.com/subject/vxgbdktx.html
网友评论