Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • GATK unified SNP caller error

    Hi,

    Thanks for nice piece of software. I have been trying out the UnifiedGenotyper but am having some issues with it. It crashes if I specify a region larger than a few Mbs, like so:

    $ java -Xms32g -Xmx64g -jar /projects/Crisp/bin/GenomeAnalysisTK-1.0.3056/GenomeAnalysisTK.jar -l INFO -R hg18.fasta -T UnifiedGenotyper -I SRR017931.bwa_mapped.sorted.bam -varout bwa_all_test.geli -vf GELI -confidence 50.0 --platform SOLEXA -L chr20:1-5000000
    INFO 16:53:06,596 HelpFormatter - ---------------------------------------------------------------------------
    INFO 16:53:06,598 HelpFormatter - The Genome Analysis Toolkit (GATK) v1.0.3056, Compiled 2010/03/22 12:02:00
    INFO 16:53:06,598 HelpFormatter - Copyright (c) 2009 The Broad Institute
    INFO 16:53:06,599 HelpFormatter - Please view our documentation at http://www.broadinstitute.org/gsa/wiki
    INFO 16:53:06,599 HelpFormatter - For support, please view our support site at http://getsatisfaction.com/gsa
    INFO 16:53:06,599 HelpFormatter - Program Args: -l INFO -R hg18.fasta -T UnifiedGenotyper -I SRR017931.bwa_mapped.sorted.bam -varout bwa_all_test.geli -vf GELI -confidence 50.0 --platform SOLEXA -L chr20:1-5000000
    INFO 16:53:06,599 HelpFormatter - Date/Time: 2010/03/26 16:53:06
    INFO 16:53:06,600 HelpFormatter - ---------------------------------------------------------------------------
    INFO 16:53:06,607 GenomeAnalysisEngine - Strictness is SILENT
    INFO 16:53:06,674 ReferenceOrderedData - * Adding rod class GFF
    INFO 16:53:06,675 ReferenceOrderedData - * Adding rod class dbSNP
    INFO 16:53:06,676 ReferenceOrderedData - * Adding rod class HapMapAlleleFrequencies
    INFO 16:53:06,676 ReferenceOrderedData - * Adding rod class SAMPileup
    INFO 16:53:06,677 ReferenceOrderedData - * Adding rod class GELI
    INFO 16:53:06,677 ReferenceOrderedData - * Adding rod class RefSeq
    INFO 16:53:06,678 ReferenceOrderedData - * Adding rod class Table
    INFO 16:53:06,679 ReferenceOrderedData - * Adding rod class PooledEM
    INFO 16:53:06,680 ReferenceOrderedData - * Adding rod class CleanedOutSNP
    INFO 16:53:06,680 ReferenceOrderedData - * Adding rod class Sequenom
    INFO 16:53:06,681 ReferenceOrderedData - * Adding rod class SangerSNP
    INFO 16:53:06,682 ReferenceOrderedData - * Adding rod class SimpleIndel
    INFO 16:53:06,682 ReferenceOrderedData - * Adding rod class PointIndel
    INFO 16:53:06,683 ReferenceOrderedData - * Adding rod class HapMapGenotype
    INFO 16:53:06,683 ReferenceOrderedData - * Adding rod class Intervals
    INFO 16:53:06,684 ReferenceOrderedData - * Adding rod class Variants
    INFO 16:53:06,685 ReferenceOrderedData - * Adding rod class GLF
    INFO 16:53:06,686 ReferenceOrderedData - * Adding rod class VCF
    INFO 16:53:06,686 ReferenceOrderedData - * Adding rod class PicardDbSNP
    INFO 16:53:06,687 ReferenceOrderedData - * Adding rod class HapmapVCF
    INFO 16:53:06,688 ReferenceOrderedData - * Adding rod class Beagle
    INFO 16:53:06,689 ReferenceOrderedData - * Adding rod class Plink
    INFO 16:53:06,690 ReferenceOrderedData - Processing ROD bindings: 1 -> interval : Intervals : chr20:1-5000000
    INFO 16:53:06,690 ReferenceOrderedData - Created binding from interval to chr20:1-5000000 of type class org.broadinstitute.sting.gatk.refdata.IntervalRod
    INFO 16:53:06,695 MicroScheduler - Creating linear microscheduler
    INFO 16:53:06,702 GenomeAnalysisEngine - Processing 5000000 bases in intervals
    INFO 16:53:07,459 TraversalEngine - [PROGRESS] Traversed to chr20:1, processing 1 sites in 0.76 secs (760000.00 secs per 1M sites)
    #
    # An unexpected error has been detected by Java Runtime Environment:
    #
    # SIGSEGV (0xb) at pc=0x00002b8dbc15c821, pid=11809, tid=1085090112
    #
    # Java VM: OpenJDK 64-Bit Server VM (1.6.0-b09 mixed mode linux-amd64)
    # Problematic frame:
    # V [libjvm.so+0x20a821]
    #
    # An error report file with more information is saved as:
    # /projects/Crisp/bwa_vs_bfast/hs_err_pid11809.log
    #
    # If you would like to submit a bug report, please visit:
    # http://icedtea.classpath.org/bugzilla
    # The crash happened outside the Java Virtual Machine in native code.
    # See problematic frame for where to report the bug.
    #
    Aborted

    Any ideas on this? It occurs if I use -L and specify a region > a few Mbs or if I specify no size at all. I works fine if U use -L chr1:1-500,000 for example. Any ideas on this?

    Thanks
    Daniel

  • #2
    I had a similar problem. The GATK people advised me to update the version of java, and this fixed the issue. pfs

    Comment


    • #3
      Hi,

      I updated the java version to jre1.6.0_18 and it did the trick.

      Thanks a lot for the info.
      Daniel

      Comment


      • #4

        Comment


        • #5
          Hi,

          the image i posted doesn't appear. so i have to write my question!
          my code is as follows:
          [renml@microcystis ~]$ java -jar GenomeAnalysisTK.jar -R exampleFASTA.fasta -I exampleBAM.bam -T CountLoci
          INFO 14:03:05,468 HelpFormatter - ---------------------------------------------------------------------------
          INFO 14:03:05,470 HelpFormatter - The Genome Analysis Toolkit (GATK) v1.0.4705, Compiled 2010/11/18 14:06:00
          INFO 14:03:05,470 HelpFormatter - Copyright (c) 2010 The Broad Institute
          INFO 14:03:05,470 HelpFormatter - Please view our documentation at http://www.broadinstitute.org/gsa/wiki
          INFO 14:03:05,471 HelpFormatter - For support, please view our support site at http://getsatisfaction.com/gsa
          INFO 14:03:05,471 HelpFormatter - Program Args: -R exampleFASTA.fasta -I exampleBAM.bam -T CountLoci
          INFO 14:03:05,471 HelpFormatter - Date/Time: 2010/11/24 14:03:05
          INFO 14:03:05,471 HelpFormatter - ---------------------------------------------------------------------------
          INFO 14:03:05,471 HelpFormatter - ---------------------------------------------------------------------------
          INFO 14:03:05,472 AbstractGenomeAnalysisEngine - Strictness is SILENT
          INFO 14:03:05,588 TraversalEngine - [INITIALIZATION COMPLETE; TRAVERSAL STARTING]
          2026
          INFO 14:03:05,718 TraversalEngine - [PROGRESS] Traversed 2,026 sites in 0.16 secs (79.47 secs per 1M sites)
          INFO 14:03:05,719 TraversalEngine - Total runtime 0.16 secs, 0.00 min, 0.00 hours
          INFO 14:03:05,720 TraversalEngine - 0 reads were filtered out during traversal out of 33 total (0.00%)
          INFO 14:03:05,722 GATKRunReport - Aggregating data for run report



          java version "1.6.0_21"
          Java(TM) SE Runtime Environment (build 1.6.0_21-b06)
          Java HotSpot(TM) 64-Bit Server VM (build 17.0-b16, mixed mode


          could anyone help me?

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Strategies for Sequencing Challenging Samples
            by seqadmin


            Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
            03-22-2024, 06:39 AM
          • seqadmin
            Techniques and Challenges in Conservation Genomics
            by seqadmin



            The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

            Avian Conservation
            Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
            03-08-2024, 10:41 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, Yesterday, 06:37 PM
          0 responses
          10 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, Yesterday, 06:07 PM
          0 responses
          9 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-22-2024, 10:03 AM
          0 responses
          51 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-21-2024, 07:32 AM
          0 responses
          67 views
          0 likes
          Last Post seqadmin  
          Working...
          X