perl 使用正则表达式检查有效的 DNA 序列

Question

我想编写一个以 FASTA 文件作为参数并打印出序列（不带标题）的子例程。 子程序应检查序列是否包含除 DNA 碱基（A、T、G、C）以外的任何其他字母。

这是我的代码：

scalar_sequence ("sequence.fa");

sub scalar_sequence {
    my $file = $_[0];
    my $sequence;
    open (READ, $file) || die "Cannot open $file: $!.\n";
    while (<READ>){
        if (/^>/){
            next;
        } 
        if (/^[ATCG]/){
            $sequence .= $_;
        } else {
            die "invalid sequence\n";
        }
    }
    print $sequence, "\n";
}

当我运行此代码时，我得到“无效序列”作为输出。 当我将“else”省略时，即使序列包含另一个字母，它也会打印出序列。

有什么问题？

提前致谢！

Answer 1

问题出在/^[ATCG]/这行应该是/^[ATCG]+$/

你的代码应该是

chomp;  
next if (/^>/); # skip for header
next if(/^\s*$/);  #skip for empty line
if (/^[ATCG]+$/){
        $sequence .= $_;
    } else {
        die "invalid sequence\n";
    }

您只考虑以 A 或 T 或 G 或 C 开头的行的开头。您应该扩展匹配项。