Shell脚本中如何将Awk输出作为twoBitToFa的输入参数？

1. 直接管道传递的正确命令（无中间文件） #

问题核心是twoBitToFa对输入格式要求严格，且需要明确指定从标准输入读取BED内容。用以下命令即可解决：

awk -v OFS="\t" '{print $1, $2, $3, $1":"$2"-"$3}' "$P/sample.bed" | twoBitToFa "$twobit" - output.fa

关键细节： #

• -v OFS="\t"：强制Awk用制表符分隔输出列，符合BED格式规范（twoBitToFa不接受空格分隔的BED内容）
• 第四列用$1":"$2"-"$3生成唯一标识（如chr1:100-200），保证每个区域唯一
• twoBitToFa命令中的-：代表从标准输入读取BED数据，这是避免错误读取全基因组的核心参数

如果需要压缩输出，可追加管道：

awk -v OFS="\t" '{print $1, $2, $3, $1":"$2"-"$3}' "$P/sample.bed" | twoBitToFa "$twobit" - stdout | gzip > output.fa.gz

2. 保存Awk输出到变量并传递的方法 #

若必须将BED内容暂存到变量，需用双引号包裹变量以保留换行符：

# 捕获Awk输出，保留换行
bed_content=$(awk -v OFS="\t" '{print $1, $2, $3, $1":"$2"-"$3}' "$P/sample.bed")
# 传递给twoBitToFa
printf "%s" "$bed_content" | twoBitToFa "$twobit" - output.fa

注意： #

• 直接用$()捕获输出时，Shell默认会把换行替换为空格，必须用双引号包裹变量才能保留原始格式
• 优先用printf "%s"而非echo -e，避免转义字符解析错误

3. 排查要点 #

• 检查原始sample.bed：确保3列是制表符分隔，无多余空格或空行
• 验证Awk输出：单独运行Awk命令，确认每行是四列制表符分隔的标准BED格式
• 确认twobit变量：指向正确的hg38.2bit文件路径

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