美文网首页
perl 练习:DNA反向互补

perl 练习:DNA反向互补

作者: 余绕 | 来源:发表于2021-05-10 00:12 被阅读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";
}

输入内容

image.png

输出内容

image.png

相关文章

网友评论

      本文标题:perl 练习:DNA反向互补

      本文链接:https://www.haomeiwen.com/subject/vxgbdktx.html