在对输入文件建立索引时报错: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
报错: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文件做索引嘛;
但我尝试用