Error while using GATK Haplotypecaller
0
0
Entering edit mode
2.3 years ago

Hello, I want to create a single vcf file with 2 bam files using GATK's Haplotypecaller. I used the following command :

gatk HaplotypeCaller -R reference.fa -I all_bams.list --min-base-quality-score 20 -O bulk.vcf

where my all_bams.list looks like this:

bam1.bam
bam2.bam

I am runnnig this command from the same folder where the bam files are present. However I am getting the following error:

java.lang.IllegalStateException: the sample list cannot be null or empty
    at org.broadinstitute.hellbender.utils.Utils.validate(Utils.java:814)
    at org.broadinstitute.hellbender.tools.walkers.genotyper.GenotypingEngine.<init>(GenotypingEngine.java:69)
    at org.broadinstitute.hellbender.tools.walkers.genotyper.MinimalGenotypingEngine.<init>(MinimalGenotypingEngine.java:46)
    at org.broadinstitute.hellbender.tools.walkers.genotyper.MinimalGenotypingEngine.<init>(MinimalGenotypingEngine.java:35)
    at org.broadinstitute.hellbender.tools.walkers.haplotypecaller.HaplotypeCallerEngine.initializeActiveRegionEvaluationGenotyperEngine(HaplotypeCallerEngine.java:328)
    at org.broadinstitute.hellbender.tools.walkers.haplotypecaller.HaplotypeCallerEngine.initialize(HaplotypeCallerEngine.java:227)
    at org.broadinstitute.hellbender.tools.walkers.haplotypecaller.HaplotypeCallerEngine.<init>(HaplotypeCallerEngine.java:184)
    at org.broadinstitute.hellbender.tools.walkers.haplotypecaller.HaplotypeCaller.onTraversalStart(HaplotypeCaller.java:259)
    at org.broadinstitute.hellbender.engine.GATKTool.doWork(GATKTool.java:1083)
    at org.broadinstitute.hellbender.cmdline.CommandLineProgram.runTool(CommandLineProgram.java:140)
    at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMainPostParseArgs(CommandLineProgram.java:192)
    at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:211)
    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)

I don't understand what's wrong. Is something wrong with my bam file or should I change all_bams.list?

And if I add path of the bam files in all_bams.list then It throws another error:

A USER ERROR has occurred: Couldn't read file. Error was: /path/to/bam/bam1.bam with exception: Cannot read non-existent file: file:///current/folder/path/to/bam/bam1.bam

Please help me out. Thanks in advance.

Gatk Haplotypecaller Bam • 630 views
ADD COMMENT

Login before adding your answer.

Traffic: 2046 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6