使用Nextflow运行流程时出现fail to read the header from "-"错误求助

核心原因分析 #

错误里的两个提示存在关联：

• 容器平台不匹配（linux/amd64 vs linux/arm64/v8）：arm64架构主机上运行amd64容器时，若未开启合适的架构模拟，会导致bwa等工具运行异常，无法输出有效的SAM格式数据。
• [main_samview] fail to read the header from "-"：samtools从标准输入读取不到有效SAM头部，说明上游的bwa mem没有生成正常输出。

手动执行.command.run正常，大概率是因为手动启动容器时Docker自动启用了qemu架构模拟，而Nextflow默认配置未明确指定平台，导致模拟机制未生效。

具体解决步骤 #

1. 强制指定容器运行平台 #

在align进程中添加platform参数，让Docker强制以amd64架构运行容器，借助qemu完成架构兼容：

process align {
   container = 'quay.io/biocontainers/mulled-v2-fe8faa35dbf6dc65a0f7f5d4ea12e31a79f73e40:8110a70be2bfe7f75a2ea7f2a89cda4cc7732095-0'
   platform = 'linux/amd64'  # 添加该参数指定运行平台

   input:
   tuple val(sampleId), val(path),path(read1), path(read2)
   tuple path(file1),path(file2)
  
   exec:
   String strBam = sampleId + '.bam'
  
   script:
   """
   bwa mem -v 0 $params.align_ref/hg19.fa $file1 $file2 | samtools view -Sb -u > $strBam
   """
}

也可在全局配置文件nextflow.config中设置默认平台，避免每个进程重复配置：

docker {
    enabled = true
    platform = 'linux/amd64'
}

2. 验证工具输出有效性 #

若指定平台后问题仍存在，可临时拆分管道排查具体环节：

bwa mem -v 2 $params.align_ref/hg19.fa $file1 $file2 > temp.sam
echo "bwa输出行数: $(wc -l temp.sam)"
samtools view -Sb -u temp.sam > $strBam
rm temp.sam

通过该脚本可以查看bwa是否生成了有效SAM文件：如果temp.sam为空或无头部信息，需进一步确认$params.align_ref/hg19.fa在容器内是否可访问，以及file1/file2是否为有效测序文件。

3. 替换为arm64兼容镜像 #

如果qemu模拟带来性能问题，可寻找支持arm64架构的bwa+samtools镜像，比如使用biocontainers官方的arm64镜像，或自行构建适配arm64的镜像，从根源解决平台兼容问题。

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