Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • GATK error

    Dear All,

    I am using GATK for first time. i installed galaxy locally and running the softwares in it.
    I am using Unified Genotyper module. my Gatk version is (GenomeAnalysisTK-2.5-2-gf57256b)
    I am getting error
    ##### ERROR ------------------------------------------------------------------------------------------
    ##### ERROR A USER ERROR has occurred (version 2.5-2-gf57256b):
    ##### 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 website and forum for extensive documentation and answers to
    ##### ERROR commonly asked questions http://www.broadinstitute.org/gatk
    ##### ERROR
    ##### ERROR MESSAGE: The fasta file you specified (/tmp/tmp-gatk-NYjkbj) does not exist.

    i mentioned the path to human genome fasta file in .loc file but again i am getting same error..

    Could anyone suggest on this...

    Thanks
    Sridhar

  • #2
    What is the exact command that you used?

    Comment


    • #3
      Hi,
      Actually i am running it in galaxy interface. I am using Unified Genotyper module in GATK.
      in backend i got the command.

      python /illumina/data/galaxy/galaxy-dist/tools/gatk/gatk_wrapper.py --max_jvm_heap_fraction "1" --stdout "/illumina/data/galaxy/galaxy-dist/database/files/000/dataset_266.dat" -d "-I" "/illumina/data/galaxy/galaxy-dist/database/files/000/dataset_150.dat" "bam" "gatk_input_0" -d "" "/illumina/data/galaxy/galaxy-dist/database/files/_metadata_files/000/metadata_11.dat" "bam_index" "gatk_input_0" -p 'java -jar "/illumina/data/galaxy/galaxy-dist/tool-data/shared/jars/gatk/GenomeAnalysisTK.jar" -T "UnifiedGenotyper" --num_threads 4 --out "/illumina/data/galaxy/galaxy-dist/database/files/000/dataset_264.dat" --metrics_file "/illumina/data/galaxy/galaxy-dist/database/files/000/dataset_265.dat" -R "" --genotype_likelihoods_model "BOTH" --standard_min_confidence_threshold_for_calling "30.0" --standard_min_confidence_threshold_for_emitting "30.0"

      Comment


      • #4
        Would be easy to run the jar file than running on the Galaxy. You need to have a .dict (dictionary) file of the fasta file used for alignment. Picards tools can be used to make .dict file. You can call population level SNPs provided you have Read Groups present in the bam files.

        java -jar path to Genome analysisTK -T UnifiedGenotyper -b bam file or path

        Check the GATK forums or use -h to get details.

        java -jar path to Genome analysisTK -T UnifiedGenotyper -h

        Comment


        • #5
          Hi newbie,

          thanks for you reply..

          you mean instead of using the human genome fasta file i can use the .dict file while running the Unified Genotyper?
          Also could you tell how to have "Read Groups present in the bam files." ??

          i saw in galaxy interface the dbSNP file is needed?? to run the software dbSNP vcf file is not required??
          Last edited by sridhar28; 06-13-2013, 01:34 AM.

          Comment


          • #6
            You need both the files (fasta and .dict) to be present in the same folder/location.

            Picard tools (http://picard.sourceforge.net) will help you with adding RG groups to your bam files. Use the tool AddOrReplaceReadGroups to achieve it.

            You don't need a VCF file from dbSNP. If you give it, it only gives the SNP ids, telling that they have been previously found or unique. GATK can call SNPS without this file (dbSNP).

            Comment


            • #7
              Thanks for your suggestions. after AddorReplaceGroup from picard running GATK ..

              the output is like

              #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT human
              INFO 12:23:28,619 ProgressMeter - chr1:5545193 5.55e+06 30.0 s 5.0 s 0.2% 4.6 h 4.6 h
              INFO 12:23:58,620 ProgressMeter - chr1:10890877 1.09e+07 60.0 s 5.0 s 0.4% 4.7 h 4.7 h
              INFO 12:24:28,621 ProgressMeter - chr1:16540257 1.65e+07 90.0 s 5.0 s 0.5% 4.7 h 4.6 h
              INFO 12:24:58,622 ProgressMeter - chr1:22062765 2.21e+07 120.0 s 5.0 s 0.7% 4.7 h 4.6 h
              INFO 12:25:28,623 ProgressMeter - chr1:27533021 2.75e+07 2.5 m 5.0 s 0.9% 4.7 h 4.6 h
              INFO 12:25:58,624 ProgressMeter - chr1:33110681 3.31e+07 3.0 m 5.0 s 1.1% 4.7 h 4.6 h

              So the VCF file will be created finally after completing all process or i should give any additional option in command line to get the VCF output??

              Thanks
              Sridhar
              Last edited by sridhar28; 06-17-2013, 02:53 AM. Reason: grammer

              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