Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • wanguan2000
    Member
    • Nov 2010
    • 24

    #1

    GATK IndelGenotyperV2 ERROR

    when I followed the pipeline:
    http://www.broadinstitute.org/gsa/wi...Whole_exome_v1

    java -Xmx4g -jar GenomeAnalysisTK.jar \
    -l INFO \
    -nt 8 \
    -R ./hg19fasta/ucsc.hg19.fasta \
    -B:dbsnp,vcf /share/user_data/wfz1/database/hg19database/dbsnp_132.hg19.vcf \
    -I ./hg19fasta/pair-end32.sorted.bam \
    -T CountCovariates \
    -cov ReadGroupCovariate \
    -cov QualityScoreCovariate \
    -cov CycleCovariate \
    -cov DinucCovariate \
    -recalFile ./hg19fasta/32.csv \
    --default_read_group HWI-ST499:5:5:8126:176904#CTTGTA \
    --default_platform illumina

    .......
    in this step:

    java -Xmx4g -jar GenomeAnalysisTK.jar \
    -T IndelGenotyperV2 \
    -l INFO \
    -R ./hg19fasta/ucsc.hg19.fasta \
    -I ./hg19fasta/cleaned31.bam \
    -bed my.brief.output.bed \
    -o my.detailed.output.vcf

    I got this ERROR:
    Badly formed genome loc: Parameters to GenomeLocParser are incorrect:The stop position 1 is less than start 150
    what is problem? I'm just following the pipeline.
    ps:my bam data is : pair-end31.sorted.bam
    Last edited by wanguan2000; 07-24-2011, 06:39 PM.
  • wanguan2000
    Member
    • Nov 2010
    • 24

    #2
    java -Xmx4g -jar /share/user_data/wfz1/wangguan/GATK/GenomeAnalysisTK-1.1-3-g1f8fc4a/GenomeAnalysisTK.jar \
    -T SomaticIndelDetector \
    -l INFO \
    -R ./hg19fasta/ucsc.hg19.fasta \
    -I ./hg19fasta/cleaned31.bam \
    -bed ./hg19fasta/my.brief.output31.bed\
    -verbose ./hg19fasta/my.detailed.output31.txt \
    -o ./hg19fasta/my.output31.vcf

    GenomeAnalysisTK-1.1-3-g1f8fc4a don't contain the IndelGenotyperV2, It has the SomaticIndelDetector.
    in my case, the pair-end data was merger in one pair-end31.sam. The SomaticIndelDetector was need -I:normal data1.bam -I:tumor data2.bam.
    what can I do?

    Comment

    • RDW
      Member
      • Oct 2008
      • 63

      #3
      'Paired sample mode' in the Somatic Indel Detector docs refers to tumour/normal pairs. If you don't have these, you can follow the instructions for 'Single Sample Mode' (the single input bam file can contain paired end data). However, in this case you'd be much better off following the current pipeline:

      http://www.broadinstitute.org/gsa/wi...th_the_GATK_v2

      This uses the Unified Genotyper to call both indels and SNPs. For help with this, see:



      The old v1 pipeline you were using is not designed to work with recent versions of the GATK.

      Comment

      Latest Articles

      Collapse

      • SEQadmin2
        Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
        by SEQadmin2



        Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

        Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
        ...
        10-14-2026, 09:55 AM
      • SEQadmin2
        New Genomics Technologies Take Aim at Long-Standing Limits
        by SEQadmin2


        Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

        We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
        ...
        09-28-2026, 10:25 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, Yesterday, 11:13 AM
      0 responses
      17 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, Yesterday, 10:48 AM
      0 responses
      18 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-15-2026, 10:52 AM
      0 responses
      36 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-14-2026, 05:39 AM
      0 responses
      42 views
      0 reactions
      Last Post SEQadmin2  
      Working...