UNIX环境下合并含相同barcode标识的fastq文件问题求助

核心语法错误修正 #

你脚本里的{2 .. 96}写法有问题——bash的大括号范围扩展不允许起始/结束数字之间有空格，正确格式是{2..96}，这是导致脚本无法正常运行的直接原因。

优化后的基础脚本 #

先修复语法问题，同时调整通配符匹配更精准，避免误匹配，并且用>替代>>防止重复运行时追加重复内容：

for i in {2..96};
do
   cat *_barcode"$i"_*.fastq > barcode"$i".fastq;
done

这里把通配符改成*_barcode"$i"_*.fastq，是为了严格匹配XX_barcodeN_XX.fastq的格式，避免意外匹配到类似barcode100这类不符合规则的文件名。

增加健壮性的进阶脚本 #

如果存在某个barcode没有对应文件的情况，上面的脚本会输出cat的错误信息。可以加入文件存在检查，同时增加运行提示：

for i in {2..96};
do
   # 把匹配的文件存入数组
   target_files=(*_barcode"$i"_*.fastq)
   # 检查数组非空且第一个元素是文件（避免匹配到目录）
   if [ ${#target_files[@]} -gt 0 ] && [ -f "${target_files[0]}" ]; then
       cat "${target_files[@]}" > barcode"$i".fastq
       echo "已完成合并: barcode$i.fastq"
   else
       echo "无对应文件跳过: barcode$i"
   fi
done

额外注意事项 #

• 如果你的文件名存在变种（比如部分文件是XX_barcodeN.fastq，没有末尾的_XX），需要把通配符改回*_barcode"$i"*.fastq来适配。
• 运行脚本前，建议先用echo *_barcode2_*.fastq测试下通配符是否能正确匹配到目标文件，确保规则无误。

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