我有一个如下所示的输入文件:
>Seq_1;1
AAAAAAAAAAAAAAAAAAAAA
>Seq_2;1
CCCCCCCCCCCCCCCCCCCCC
还有更多这样的成对线。我想要做的是像这样简单地打印出来:
>Seq_1;1 AAAAAAAAAAAAAAAAAAAAA
>Seq_2;1 CCCCCCCCCCCCCCCCCCCCC
但是为什么这段代码会失败:
#!/usr/bin/perl -w
while ( <> ) {
chomp;
my $line = $_;
my $rdn = "";
my $sq = "";
if ( $line =~ /^>/ ) {
$rdn = $line;
}
elsif ($line =~ /^[ATCG]/) {
$sq = $line;
}
print "$rdn $sq\n";
}
它改为打印:
>Seq_1;1
AAAAAAAAAAAAAAAAAAAAA
>Seq_2;1
CCCCCCCCCCCCCCCCCCCCC