批量处理salmon生成的quant.sf，提取lnc条目至对应目录lnc.sf

你的脚本问题在于每次循环都将结果输出到当前工作目录的lnc.sf，导致后续循环的结果覆盖之前的内容，最终只保留最后一次的输出，且没有在每个quant.sf所在目录生成对应文件。

以下是两种可行的解决方案：

方法一：循环中定位目标目录 #

find . -name "quant.sf" | while read quant_file
do
  # 获取当前quant.sf文件所在的目录路径
  target_dir=$(dirname "$quant_file")
  # 直接用grep读取文件并输出到对应目录的lnc.sf
  grep 'lnc' "$quant_file" > "$target_dir/lnc.sf"
done

• 用dirname "$quant_file"精准定位每个quant.sf的所在目录
• 去掉冗余的cat命令，直接让grep读取文件，提升执行效率
• 变量用双引号包裹，避免文件名含空格或特殊字符时出现错误

方法二：用find的-exec选项简化脚本 #

find . -name "quant.sf" -exec sh -c 'grep "lnc" "$0" > "$(dirname "$0")/lnc.sf"' {} \;

• 利用find的-exec参数直接处理每个找到的文件
• {}代表当前找到的quant.sf文件，sh -c执行后续命令逻辑，实现和方法一完全一致的效果

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