The current GATK version is 3.6-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Powered by Vanilla. Made with Bootstrap.

HaplotypeCaller --fullHaplotype not outputting full haplotype

chochochocho Posts: 1Member
edited August 2012 in Ask the GATK team

I am also having a problem with HaplotypeCaller - when I set the flag to print out full haplotypes, no haplotype files are created. Here is my command:

java -jar /home/chodon/applications/GATK/GenomeAnalysisTK.jar -T HaplotypeCaller -R /home/chodon/comb_eny/abyss/final_contigs/A2/Assembly.fasta -I /home/chodon/comb_eny/abyss/final_contigs/A2/bowtie2_out/A2.sorted.bam --fullHaplotype --ignoreLaneInfo -o /home/chodon/comb_eny/abyss/final_contigs/A2/GATK_haplo2.out

Any help is much appreciated! Thanks, C

Tagged:

Best Answers

Answers

  • chochochocho Posts: 1Member

    Hi, Thanks for your quick reply!

    When I added --genotypeFullActiveRegion, full command here:

    java -jar /home/chodon/applications/GATK/GenomeAnalysisTK.jar -T HaplotypeCaller -R /home/chodon/comb_eny/abyss/final_contigs/A2/Assembly.fasta -I /home/chodon/comb_eny/abyss/final_contigs/A2/bowtie2_out/A2.sorted.bam --genotypeFullActiveRegion --fullHaplotype --ignoreLaneInfo -o /home/chodon/comb_eny/abyss/final_contigs/A2/GATK_haplo3.out

    I got an error message: Code exception - here is the stack trace

    ERROR stack trace

    java.lang.NullPointerException
    at org.broadinstitute.sting.gatk.walkers.haplotypecaller.HaplotypeCaller.map(HaplotypeCaller.java:428)
    at org.broadinstitute.sting.gatk.walkers.haplotypecaller.HaplotypeCaller.map(HaplotypeCaller.java:104)
    at org.broadinstitute.sting.gatk.traversals.TraverseActiveRegions.processActiveRegion(TraverseActiveRegions.java:246)
    at org.broadinstitute.sting.gatk.traversals.TraverseActiveRegions.callWalkerMapOnActiveRegions(TraverseActiveRegions.java:202)
    at org.broadinstitute.sting.gatk.traversals.TraverseActiveRegions.processActiveRegions(TraverseActiveRegions.java:177)
    at org.broadinstitute.sting.gatk.traversals.TraverseActiveRegions.traverse(TraverseActiveRegions.java:134)
    at org.broadinstitute.sting.gatk.traversals.TraverseActiveRegions.traverse(TraverseActiveRegions.java:27)
    at org.broadinstitute.sting.gatk.executive.LinearMicroScheduler.execute(LinearMicroScheduler.java:62)
    at org.broadinstitute.sting.gatk.GenomeAnalysisEngine.execute(GenomeAnalysisEngine.java:269)
    at org.broadinstitute.sting.gatk.CommandLineExecutable.execute(CommandLineExecutable.java:113)
    at org.broadinstitute.sting.commandline.CommandLineProgram.start(CommandLineProgram.java:236)
    at org.broadinstitute.sting.commandline.CommandLineProgram.start(CommandLineProgram.java:146)
    at org.broadinstitute.sting.gatk.CommandLineGATK.main(CommandLineGATK.java:93)

    Thanks again!
    C

  • sbrubakersbrubaker Posts: 4Member

    Hi, I am trying the same things, and I have version 2.1. I get an error that says "Unexpected bases in allele bases"

  • Geraldine_VdAuweraGeraldine_VdAuwera Posts: 10,410Administrator, Dev admin

    Hi @sbrubaker, have you tried validating your bam file? That type of error is often associated with invalid bams.

    Geraldine Van der Auwera, PhD

Sign In or Register to comment.