GATK BaseRecalibrator--known-sites参数缺失报错求助

错误根源分析 #

你遇到的"Argument 'known-sites' is required"错误，本质是**--known-sites参数后没有提供有效的文件路径**，两个命令的问题分别是：

• 第一个命令中${Known_Sites}变量未正确赋值（为空），导致参数无对应值；
• 第二个直接调用jar的命令里，不仅--known-sites没跟路径，连-R（参考基因组）参数也遗漏了对应路径，语法完全错误。

具体解决步骤 #

1. 检查并修复变量赋值
   先确认命令中用到的变量是否正确定义，比如用echo查看${Known_Sites}的值：

   echo ${Known_Sites}
   

   如果输出为空，先通过export定义变量（替换为你的实际文件路径）：

   export Known_Sites="/path/to/your/known_variants.vcf"
   export reference="/path/to/reference_genome.fasta"  # 注意你之前拼写的referance是错误的，应该是reference
   export aligned_reads="/path/to/aligned_reads_dir"
   export data1="/path/to/output_dir"
   

2. 修正gatk命令（推荐使用官方封装的gatk命令）
   确保所有必填参数都有对应值，正确命令示例：

   gatk BaseRecalibrator \
     -I ${aligned_reads}/file_copy_sorted.bam \
     -R ${reference} \
     --known-sites ${Known_Sites} \
     -O ${data1}/sorted_Read1.table
   

3. 直接调用jar包的正确写法
   若必须直接调用jar，每个参数后都要跟实际路径，不能只写参数名：

   java -Dsamjdk.use_async_io_read_samtools=false \
     -Dsamjdk.use_async_io_write_samtools=true \
     -Dsamjdk.use_async_io_write_tribble=false \
     -Dsamjdk.compression_level=2 \
     -jar /home/hi/Documents/gatk-4.4.0.0/gatk-package-4.4.0.0-local.jar BaseRecalibrator \
     -I /path/to/file_copy_sorted.bam \
     -R /path/to/reference_genome.fasta \
     --known-sites /path/to/known_variants.vcf \
     -O /path/to/sorted_Read1.table
   

   注意：--known-sites需要传入已知变异位点的VCF文件（比如dbSNP数据集），不能留空。

4. 验证参数格式
   若对参数有疑问，可查看工具的详细帮助文档：

   gatk BaseRecalibrator --help
   

   或直接调用jar时查看帮助：

   java -jar /home/hi/Documents/gatk-4.4.0.0/gatk-package-4.4.0.0-local.jar BaseRecalibrator --help
   

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