Perl对比两txt文件匹配子串提取指定行代码无输出求助

核心错误原因 #

• 读取文件行未去除换行符：使用<>操作符读取文本行时，每行末尾会自带换行符\n，你直接用带换行的extline作为子串去gffline中匹配，自然无法命中不含对应换行的目标子串，需要先对pegs.txt的每行执行chomp去除末尾换行。
• 错误使用exit关键字：你在内层循环匹配到第一条结果后直接调用exit，会直接终止整个Perl进程，后续打开输出文件、写入内容的代码完全不会执行，这是输出文件为空的直接原因。此处不需要exit，如果要避免同一行匹配多个关键词重复写入，可在匹配成功后用last跳出内层循环即可。
• 多余添加换行符：gffline本身已经自带读取时的换行符，不需要额外拼接\n，否则输出结果会多出空行。

修正后的完整代码 #

#!/usr/bin/perl -w
use strict;

open (FILE, "<combined_gff.txt") or die "无法打开combined_gff.txt: $!";
my @gff = <FILE>;
close FILE;

open (DATA, "<pegs.txt") or die "无法打开pegs.txt: $!";
my @ext = <DATA>;
close DATA;
chomp @ext; # 批量去除pegs所有行的末尾换行

my $str = '';

foreach my $gffline (@gff) {
    foreach my $extline (@ext) {
        if ( index($gffline, $extline) != -1) {
            $str .= $gffline;
            last; # 匹配到一个关键词就跳出内层循环，避免重复添加
        }
    }
}

open (OUT, ">", "output.txt") or die "无法创建output.txt: $!";
print OUT $str;
close (OUT);

可选优化建议 #

• 如果combined_gff.txt体积较大，不需要一次性读入内存，可以逐行读取逐行判断，降低内存占用。
• die语句中添加$!可以输出系统级错误原因，方便排查文件权限、路径不存在等问题。

内容的提问来源于stack exchange，提问作者Niravra Chakraborty