美文网首页
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