使用bioawk替换无键值对的fasta表头失败，变量为何无法生效？

为什么$prefix无法生效？ #

你用单引号包裹了bioawk的处理代码，而shell会完全忽略单引号内的变量替换，导致$prefix没有被替换成shell变量的实际值，反而被bioawk当成了它自身的字段引用（类似$1、$2这种格式），但bioawk中不存在名为prefix的字段，因此抛出错误。

两种可行的解决办法 #

方法1：改用双引号包裹bioawk代码（需转义awk自身变量） #

把包裹bioawk代码的单引号换成双引号，让shell先完成${prefix}的变量替换，再将代码传递给bioawk。注意要给awk自身的变量（比如$seq）加上反斜杠转义，避免被shell误解析：

for infile in *.faa
do
        prefix=$(basename "$infile" .faa)
        bioawk -c fastx "{ print \">${prefix}_\" ++i \"\n\"\$seq }" < "$infile" > "${prefix}_hdrrn.faa"
done

方法2：用-v参数向bioawk传递变量（更推荐） #

bioawk支持-v var=value的方式从外部传入变量，这种方式不需要处理复杂的转义，也更安全（比如当prefix包含特殊字符时不会出现异常）：

for infile in *.faa
do
        prefix=$(basename "$infile" .faa)
        bioawk -v prefix="$prefix" -c fastx '{ print ">"prefix"_" ++i "\n"$seq }' < "$infile" > "${prefix}_hdrrn.faa"
done

额外补充：原脚本中用$(ls *.faa)遍历文件存在隐患，如果文件名包含空格或特殊字符会出错，直接用*.faa遍历文件更可靠。

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