搜索
查看: 978|回复: 3

GATK4运行BaseRecalibrator时报错

[复制链接]

12

主题

40

帖子

379

积分

中级会员

Rank: 3Rank: 3

积分
379
发表于 2018-5-16 14:14:17 | 显示全部楼层 |阅读模式
命令如下:
[Shell] 纯文本查看 复制代码
java  -jar ~/tools/gatk-4.0.3.0/gatk-package-4.0.3.0-local.jar BaseRecalibrator \
        -R ~/project/ref/hg38/Homo_sapiens_assembly38.fasta \
        -I N1.sorted.MarkDuplicates.bam  -L  ../aml_hg38.probe.bed  -ip 90 \
        --known-sites ~/tools/GATK/1000G_omni2.5.hg38.vcf \
        --known-sites ~/tools/GATK/1000G_phase1.snps.high_confidence.hg38.vcf \
        --known-sites ~/tools/GATK/dbsnp_146.hg38.vcf \
        --known-sites ~/tools/GATK/hapmap_3.3.hg38.vcf \
        --known-sites ~/tools/GATK/Mills_and_1000G_gold_standard.indels.hg38.vcf \
        -O $fname.realn.recal


错误信息如下:
14:41:18.808 INFO  NativeLibraryLoader - Loading libgkl_compression.so from jar:file:/public/home/lijie/tools/gatk-4.0.3.0/gatk-package-4.0.3.0-local.jar!/com/intel/
gkl/native/libgkl_compression.so
14:41:19.208 INFO  BaseRecalibrator - ------------------------------------------------------------
14:41:19.208 INFO  BaseRecalibrator - The Genome Analysis Toolkit (GATK) v4.0.3.0
14:41:19.208 INFO  BaseRecalibrator - For support and documentation go to https://software.broadinstitute.org/gatk/
14:41:19.209 INFO  BaseRecalibrator - Executing as lijie@node2 on Linux v2.6.32-504.el6.x86_64 amd64
14:41:19.209 INFO  BaseRecalibrator - Java runtime: OpenJDK 64-Bit Server VM v1.8.0_121-b15
14:41:19.210 INFO  BaseRecalibrator - Start Date/Time: May 16, 2018 2:41:18 PM CST
14:41:19.210 INFO  BaseRecalibrator - ------------------------------------------------------------
14:41:19.210 INFO  BaseRecalibrator - ------------------------------------------------------------
14:41:19.210 INFO  BaseRecalibrator - HTSJDK Version: 2.14.3
14:41:19.211 INFO  BaseRecalibrator - Picard Version: 2.17.2
14:41:19.211 INFO  BaseRecalibrator - HTSJDK Defaults.COMPRESSION_LEVEL : 2
14:41:19.211 INFO  BaseRecalibrator - HTSJDK Defaults.USE_ASYNC_IO_READ_FOR_SAMTOOLS : false
14:41:19.211 INFO  BaseRecalibrator - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_SAMTOOLS : true
14:41:19.211 INFO  BaseRecalibrator - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_TRIBBLE : false
14:41:19.211 INFO  BaseRecalibrator - Deflater: IntelDeflater
14:41:19.211 INFO  BaseRecalibrator - Inflater: IntelInflater
14:41:19.211 INFO  BaseRecalibrator - GCS max retries/reopens: 20
14:41:19.211 INFO  BaseRecalibrator - Using google-cloud-java patch 6d11bef1c81f885c26b2b56c8616b7a705171e4f from https://github.com/droazen/google-cloud-java/tree/d
r_all_nio_fixes
14:41:19.212 INFO  BaseRecalibrator - Initializing engine
14:41:20.500 INFO  FeatureManager - Using codec VCFCodec to read file file:///public/home/lijie/tools/GATK/1000G_omni2.5.hg38.vcf
14:41:20.618 INFO  BaseRecalibrator - Shutting down engine
[May 16, 2018 2:41:20 PM CST] org.broadinstitute.hellbender.tools.walkers.bqsr.BaseRecalibrator done. Elapsed time: 0.03 minutes.
Runtime.totalMemory()=1707081728

htsjdk.tribble.TribbleException: Problem detecting index type
        at htsjdk.tribble.index.IndexFactory$IndexType.getIndexType(IndexFactory.java:149)
        at htsjdk.tribble.index.IndexFactory.loadIndex(IndexFactory.java:179)
        at htsjdk.tribble.TribbleIndexedFeatureReader.loadIndex(TribbleIndexedFeatureReader.java:162)
        at htsjdk.tribble.TribbleIndexedFeatureReader.hasIndex(TribbleIndexedFeatureReader.java:227)
        at org.broadinstitute.hellbender.engine.FeatureDataSource.<init>(FeatureDataSource.java:262)
        at org.broadinstitute.hellbender.engine.FeatureManager.addToFeatureSources(FeatureManager.java:202)
        at org.broadinstitute.hellbender.engine.FeatureManager.initializeFeatureSources(FeatureManager.java:182)
        at org.broadinstitute.hellbender.engine.FeatureManager.<init>(FeatureManager.java:153)
        at org.broadinstitute.hellbender.engine.ReadWalker.initializeFeatures(ReadWalker.java:73)
        at org.broadinstitute.hellbender.engine.GATKTool.onStartup(GATKTool.java:558)
        at org.broadinstitute.hellbender.engine.ReadWalker.onStartup(ReadWalker.java:55)
        at org.broadinstitute.hellbender.cmdline.CommandLineProgram.runTool(CommandLineProgram.java:132)
        at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMainPostParseArgs(CommandLineProgram.java:179)
        at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:198)
        at org.broadinstitute.hellbender.Main.runCommandLineProgram(Main.java:160)
        at org.broadinstitute.hellbender.Main.mainEntry(Main.java:203)
        at org.broadinstitute.hellbender.Main.main(Main.java:289)
Caused by: java.io.EOFException
        at htsjdk.tribble.util.LittleEndianInputStream.readInt(LittleEndianInputStream.java:66)
        at htsjdk.tribble.index.IndexFactory$IndexType.getIndexType(IndexFactory.java:137)
        ... 16 more


请问这是哪里出问题了?



上一篇:ROSALIND上的题(18.开放阅读框架)求完整注释
下一篇:RepeatMasker 工具
回复

使用道具 举报

3

主题

11

帖子

190

积分

注册会员

Rank: 2

积分
190
发表于 2018-5-16 15:30:23 | 显示全部楼层
你用的vcf文件都构建了索引文件么?没有vcf.idx文件你就用GATK的IndexFeatureFile构建一下。
你报错的信息给的不全面,至少我看不到全部。
回复 支持 反对

使用道具 举报

12

主题

40

帖子

379

积分

中级会员

Rank: 3Rank: 3

积分
379
 楼主| 发表于 2018-5-17 11:34:39 | 显示全部楼层
zhou945876833 发表于 2018-5-16 15:30
你用的vcf文件都构建了索引文件么?没有vcf.idx文件你就用GATK的IndexFeatureFile构建一下。
你报错的信息 ...

错误的地方就是你说的这个,昨天检查了一下,发现vcf文件的索引文件大小都是0 ,原来直接用GATK的IndexFeatureFile构建就可以了吗,我试过bcftools index,不过没成功,vcf必须得是bgzip压缩的,然后我又用bgzip压缩,也没成功,想问问这个gatk的IndexFeatureFile对文件有要求吗
回复 支持 反对

使用道具 举报

3

主题

11

帖子

190

积分

注册会员

Rank: 2

积分
190
发表于 2018-5-22 18:51:31 | 显示全部楼层
elaine2017 发表于 2018-5-17 11:34
错误的地方就是你说的这个,昨天检查了一下,发现vcf文件的索引文件大小都是0 ,原来直接用GATK的IndexFe ...

是的直接用GATK的IndexFeatureFile构建index
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|手机版|小黑屋|生信技能树    

GMT+8, 2018-9-23 01:44 , Processed in 0.123329 second(s), 25 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.