Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • GATK error

    java -jar gatk -T IndelGenotyperV2 -l INFO -bed brief.bed -R /data/results/genomes/hg18/hg18.fa -I fastq.sorted.realigned.bam -o fastq_indel.txt -verbose fastq_indel_statistics.txt

    This produces an error
    "GenomeLocParser are incorrect:The stop position 1 is less than start 141481674"

    does anyone know how to solve?


    INFO 20:34:49,281 HelpFormatter - ---------------------------------------------------------------------------
    INFO 20:34:49,283 HelpFormatter - The Genome Analysis Toolkit (GATK) v1.0.5974, Compiled 2011/06/10 13:26:59
    INFO 20:34:49,283 HelpFormatter - Copyright (c) 2010 The Broad Institute
    INFO 20:34:49,283 HelpFormatter - Please view our documentation at http://www.broadinstitute.org/gsa/wiki
    INFO 20:34:49,283 HelpFormatter - For support, please view our support site at http://getsatisfaction.com/gsa
    INFO 20:34:49,284 HelpFormatter - Program Args: -T IndelGenotyperV2 -l INFO -bed brief.bed -R /data/results/genomes/hg18/hg18.fa -I fastq.sorted.realigned.bam -o fastq_indel.txt -verbose fastq_indel_statistics.txt
    INFO 20:34:49,284 HelpFormatter - Date/Time: 2011/06/16 20:34:49
    INFO 20:34:49,284 HelpFormatter - ---------------------------------------------------------------------------
    INFO 20:34:49,284 HelpFormatter - ---------------------------------------------------------------------------
    INFO 20:34:49,295 GenomeAnalysisEngine - Strictness is SILENT
    INFO 20:36:05,521 IndelGenotyperV2Walker - No gene annotations available
    INFO 20:37:25,387 TraversalEngine - [INITIALIZATION COMPLETE; TRAVERSAL STARTING]
    INFO 20:37:25,388 TraversalEngine - Location processed.reads runtime per.1M.reads completed total.runtime remaining
    INFO 20:37:25,392 TraversalEngine - chr1:18958321 1.00e+00 79.9 s 132.1 w 0.6% 3.6 h 3.6 h
    ##### ERROR ------------------------------------------------------------------------------------------
    ##### ERROR A USER ERROR has occurred (version 1.0.5974):
    ##### ERROR The invalid arguments or inputs must be corrected before the GATK can proceed
    ##### ERROR Please do not post this error to the GATK forum
    ##### ERROR
    ##### ERROR See the documentation (rerun with -h) for this tool to view allowable command-line arguments.
    ##### ERROR Visit our wiki for extensive documentation http://www.broadinstitute.org/gsa/wiki
    ##### ERROR Visit our forum to view answers to commonly asked questions http://getsatisfaction.com/gsa
    ##### ERROR
    ##### ERROR MESSAGE: Badly formed genome loc: Parameters to GenomeLocParser are incorrect:The stop position 1 is less than start 141481674
    ##### ERROR ------------------------------------------------------------------------------------------

  • #2
    As far as I know bwa merges the reference sequences for aligning the reads. It somehow sounds as if there is a read which maps between two sequences. If this is the case you could clean the sam file by using cleanSam from picard.

    Comment


    • #3
      I'm actually having the same problem as well. CleanSam doesn't help in this case. Picard's ValidateSamFile does show multiple errors that I think are related:

      Mate unmapped flag does not match read unmapped flag of mate
      Mate not found for paired read

      Interestingly, using an older version of GATK (1.0.5777 in my case) works well, so I'm thinking it's just an issue of GATK being more stringent than it used to be.

      I'm hoping that Picard's FixMateInformation can fix it as I'd rather not be stuck with an older version of GATK, but at least we know that using the older version works.

      Comment


      • #4
        Hi,
        Sorry, it was a while ago but I think I solved this by checking out latest version of GATK from cvs (rather than download the one from their website). Also...

        # Needs the human genome to be sorted karyotypically. I did this by downloading from GATK site
        # You need to make a samtools faidx and a Picard (CreateSequenceDictionary.jar) .dict file for re-ording your bam file GENOME NEEDS A DICT AND FAIDX FILES

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Essential Discoveries and Tools in Epitranscriptomics
          by seqadmin


          The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist on Modified Bases...
          Yesterday, 07:01 AM
        • seqadmin
          Current Approaches to Protein Sequencing
          by seqadmin


          Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
          04-04-2024, 04:25 PM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 04-11-2024, 12:08 PM
        0 responses
        55 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-10-2024, 10:19 PM
        0 responses
        52 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-10-2024, 09:21 AM
        0 responses
        45 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-04-2024, 09:00 AM
        0 responses
        55 views
        0 likes
        Last Post seqadmin  
        Working...
        X