用GATK call SNP时报错

报错:A USER ERROR has occurred: Traversal by intervals was requested but some input files are not indexed.

Please index all input files:

samtools index SRR3274663.new.sorted.bam

Set the system property GATK_STACKTRACE_ON_USER_EXCEPTION (--java-options '-DGATK_STACKTRACE_ON_USER_EXCEPTION=true') to print the stack trace.

这是我的命令:gatk --java-options "-Xmx20g -Djava.io.tmpdir=./tmp" HaplotypeCaller -R genome.fasta -l SRR3274663.new.sorted.bam -ERC GVCF -O SRR3274663.g.vcf 1>663.log 2>&1

在这里我已经对基因组文件建了索引,这里是还要对输入的bam文件做索引嘛;

但我尝试用

samtools index SRR3274663.new.sorted.bam此命令对输入文件做索引也不成功
请先 登录 后评论

1 个回答

甜辣辣

在对输入文件建立索引时报错:Read 'SRR3274663.23205715' with ref_name='genome', ref_length=584211917, flags=113, pos=536889210 cannot be indexed

samtools index: failed to create index for "SRR3274663.new.sorted.bam": Numerical result out of range

请先 登录 后评论