GATK 报错

gatk VariantFiltration : --filter-expression 显示 Input /work/my_reseq/4.snp_indel/var_qc2/BSA_no_filter.vcf   must support random access to enable queries by interval. If it's a file, please index it using the bundled tool IndexFeatureFile该怎么处理呀

请先 登录 后评论

1 个回答

omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

vcf文件没有建立索引,你需要建立一下:

gatk  --java-options "-Xmx50g" IndexFeatureFile   -I data.vcf.gz
请先 登录 后评论