Tagged with #commandlinegatk
1 documentation article | 0 announcements | 14 forum discussions


Comments (0)

A new tool has been released!

Check out the documentation at CommandLineGATK.

No posts found with the requested search criteria.
Comments (2)

Using GATK on command-line the CatVariants command fails.

Program version: GATK 3.1-1-g07a4bf8.

ERROR MESSAGE: Invalid command line: Malformed walker argument: Could not find walker with name: CatVariants

Code to invoke:

java -jar GenomeAnalysisTK-3.1-1/GenomeAnalysisTK.jar -T CatVariant -R file.fasta

Note that in the current documentation for CatVariants the example lists the name as 'org.broadinstitute.sting.tools.CatVariants' rather than just CatVariants. Trying the listed string fails with the same error.

Comments (1)

Below is the command:

java -cp $CLASSPATH/GenomeAnalysisTK.jar org.broadinstitute.sting.tools.CatVariants \
-R GATK_ref/hg19.fasta \
-V ../GATK/VQSR/parallel_batch/raw.snps_indels-1.vcf \
-V ../GATK/VQSR/parallel_batch/raw.snps_indels-2.vcf \
-V ../GATK/VQSR/parallel_batch/raw.snps_indels-3.vcf \
-out ../GATK/VQSR/parallel_batch/combined_raw.snps_indels.vcf \
-log ../GATK/VQSR/parallel_batch/log/combined.log \
-assumeSorted

After this, the combined_raw.snps_indels.vcf file only contains the header from raw.snps_indels-1.vcf, what might be wrong?

Comments (9)

I'm running the latest GATK nightly build to process human exome-seq data (has 12 samples). It seemed be faster than the older version until I run the HaplotypeCaller. The run summary shows it will take 14 days to finish. I am wondering if there's anything in my below command: How to make it faster without losing data in the output?

java -Xmx10g -Djava.io.tmpdir=/temp/GATK_temp
-jar $CLASSPATH/GenomeAnalysisTK.jar \
-T HaplotypeCaller \
-R ../GATK_ref/hg19.fasta \
-I ./compressedbam.list \
-L ../GATK_ref/hg19knownGene_UCSC_sorted.bed \
-log ../GATK/VQSR/log/HaplotypeCaller_20131018.log \
-o ../GATK/VQSR/raw.snps_indels.vcf
Comments (2)

If I put my input files as a list in the file named "input.list", how do I set the output names? or do I just need to set the output folder and the output file names will be automatically named?

Comments (4)

I started with BWA-MEM to do alignment, used Picard to process the .SAM files (converted to bam, reorder, addorreplacegroup, etc). The GATK version I'm using is version 2.5-2-gf57256b, I cannot run 2.6 because the server only has Java 6 and I cannot upgrade it to Java 7.

I got a huge stack of error message when I run this command line (RealignerTargetCrator):

java -Xmx2g -jar $CLASSPATH/GenomeAnalysisTK.jar \ -T RealignerTargetCreator \ -R /Volumes/files/Users/user1/GATK_ref/hg19.fasta \ -I sorted_Deduped_reorder_grp.bam \ -o ./GATK/forIndelRealigner.intervals>

The error messages are these (sorry, a lot): I don't know why GATK needs to connect to window server? what permission problem? I am using a Mac OS X built server (remote). Thank you

ERROR ------------------------------------------------------------------------------------------
ERROR stack trace

java.lang.InternalError: Can't connect to window server - not enough permissions. at java.lang.ClassLoader$NativeLibrary.load(Native Method) at java.lang.ClassLoader.loadLibrary0(ClassLoader.java:1827) at java.lang.ClassLoader.loadLibrary(ClassLoader.java:1724) at java.lang.Runtime.loadLibrary0(Runtime.java:823) at java.lang.System.loadLibrary(System.java:1045) at sun.security.action.LoadLibraryAction.run(LoadLibraryAction.java:50) at java.security.AccessController.doPrivileged(Native Method) at java.awt.Toolkit.loadLibraries(Toolkit.java:1605) at java.awt.Toolkit.(Toolkit.java:1627) at sun.awt.AppContext$2.run(AppContext.java:240) at sun.awt.AppContext$2.run(AppContext.java:226) at java.security.AccessController.doPrivileged(Native Method) at sun.awt.AppContext.initMainAppContext(AppContext.java:226) at sun.awt.AppContext.access$200(AppContext.java:112) at sun.awt.AppContext$3.run(AppContext.java:306) at java.security.AccessController.doPrivileged(Native Method) at sun.awt.AppContext.getAppContext(AppContext.java:287) at com.sun.jmx.trace.Trace.out(Trace.java:180) at com.sun.jmx.trace.Trace.isSelected(Trace.java:88) at com.sun.jmx.interceptor.DefaultMBeanServerInterceptor.isTraceOn(DefaultMBeanServerInterceptor.java:1830) at com.sun.jmx.interceptor.DefaultMBeanServerInterceptor.registerDynamicMBean(DefaultMBeanServerInterceptor.java:929) at com.sun.jmx.interceptor.DefaultMBeanServerInterceptor.registerObject(DefaultMBeanServerInterceptor.java:916) at com.sun.jmx.interceptor.DefaultMBeanServerInterceptor.registerMBean(DefaultMBeanServerInterceptor.java:312) at com.sun.jmx.mbeanserver.JmxMBeanServer$2.run(JmxMBeanServer.java:1195) at java.security.AccessController.doPrivileged(Native Method) at com.sun.jmx.mbeanserver.JmxMBeanServer.initialize(JmxMBeanServer.java:1193) at com.sun.jmx.mbeanserver.JmxMBeanServer.(JmxMBeanServer.java:225) at com.sun.jmx.mbeanserver.JmxMBeanServer.(JmxMBeanServer.java:170) at com.sun.jmx.mbeanserver.JmxMBeanServer.newMBeanServer(JmxMBeanServer.java:1401) at javax.management.MBeanServerBuilder.newMBeanServer(MBeanServerBuilder.java:93) at javax.management.MBeanServerFactory.newMBeanServer(MBeanServerFactory.java:311) at javax.management.MBeanServerFactory.createMBeanServer(MBeanServerFactory.java:214) at javax.management.MBeanServerFactory.createMBeanServer(MBeanServerFactory.java:175) at sun.management.ManagementFactory.createPlatformMBeanServer(ManagementFactory.java:302) at java.lang.management.ManagementFactory.getPlatformMBeanServer(ManagementFactory.java:504) at org.broadinstitute.sting.gatk.executive.MicroScheduler.(MicroScheduler.java:222) at org.broadinstitute.sting.gatk.executive.LinearMicroScheduler.(LinearMicroScheduler.java:70) at org.broadinstitute.sting.gatk.executive.MicroScheduler.create(MicroScheduler.java:169) at org.broadinstitute.sting.gatk.GenomeAnalysisEngine.createMicroscheduler(GenomeAnalysisEngine.java:443) at org.broadinstitute.sting.gatk.GenomeAnalysisEngine.execute(GenomeAnalysisEngine.java:272) at org.broadinstitute.sting.gatk.CommandLineExecutable.execute(CommandLineExecutable.java:113) at org.broadinstitute.sting.commandline.CommandLineProgram.start(CommandLineProgram.java:245) at org.broadinstitute.sting.commandline.CommandLineProgram.start(CommandLineProgram.java:152) at org.broadinstitute.sting.gatk.CommandLineGATK.main(CommandLineGATK.java:91)

ERROR ------------------------------------------------------------------------------------------
ERROR A GATK RUNTIME ERROR has occurred (version 2.5-2-gf57256b):
Comments (6)

Dear GATK Users,

Could anybody tell me how to identify the deletions from the bam file using GATK module?? Actually i used UnifiedGenotyper i am getting list like

CHROM POS ID REF ALT QUAL FILTER INFO FORMAT human

gi|262 48155 . G A 80.77 . AC=1;AF=0.500;AN=2;BaseQRankSum=0.103;DP=10;Dels=0.00;FS=0.000;HaplotypeScore =0.0000;MLEAC=1;MLEAF=0.500;MQ=28.61;MQ0=0;MQRankSum=-1.453;QD=8.08;ReadPosRankSum=-0.336 GT:AD:DP:GQ:PL 0/1:5,5:10:99:109,0,146

Thanks Sridhar

Comments (11)

Hi all: I find that among all the work flows of GATK http://www.broadinstitute.org/gatk/guide/topic?name=methods-and-workflows there are no workflows for RNA-seq analysis. I understand that GATK mainly focuses on variant calling, can anyone tell me how to use GATK for RNA-seq analysis?

thanks daniel

Comments (8)

Hello Team,

I am attempting to run GATK's PhasebyTransmission command to phase a vcf file contains a father, mother, son trio generated from complete genomics mkvcf command.

After creating the ped file and running the command I generate the error: "MESSAGE: BUG: Attempted to get likelihoods as strings and neither the vector nor the string is set!". I am not exactly sure what this means.

When I check my file and the documentation I am able to see that the 'GL' field is contained in the file, but could this not be the case? I have attached a few lines from the vcf I am using.

Any help with resolving the this issue would be of great help.

Thank you

JumaQuar

Comments (6)

Hello, I`m new to GATK and Queue. I understand that we can write a QScript in Queue to generate separate GATK jobs and run them on a cluster of several nodes. Can we implement GATK or Queue on google hadoop?

Comments (3)

I got this error message, when trying to use a file to specify at which positions to emit variants:

ERROR MESSAGE: Couldn't read file /lustre/scratch109/sanger/tc9/agv/wgs/pipeline/union4x.positions because The interval file /lustre/scratch109/sanger/tc9/agv/wgs/pipeline/union4x.positions does not have one of the supported extensions (.bed, .list, .picard, .interval_list, or .intervals). Please rename your file with the appropriate extension. Is there a GATK page describing those 5 file formats? Some of them are unknown to me; e.g. .list.

I asked my question here, but please ignore it: http://gatkforums.broadinstitute.org/discussion/2219/l-option

Thanks a lot.

Also, the error message does not mention support for vcf files, but the documentation does. Are vcf files supported?

Comments (5)

hi all! I'm trying to complete my first GATK run, I'm doing the step in the "EXECUTION STEP" following section.

please tell me if the step execution are globally correct.

------------------------------------------------------------------ERRORS----------------------------------------------------------------

the step 4.1 isn't executed without -maxCycle 1500.

when try to execute 4.2 step I got the following error:

ERROR stack trace

org.broadinstitute.sting.utils.exceptions.ReviewedStingException: Key 1036 is too large for dimension 2 (max is 1001) at org.broadinstitute.sting.utils.collections.NestedIntegerArray.put(NestedIntegerArray.java:128) at org.broadinstitute.sting.utils.recalibration.RecalibrationReport.parseAllCovariatesTable(RecalibrationReport.java:157) at org.broadinstitute.sting.utils.recalibration.RecalibrationReport.(RecalibrationReport.java:68) at org.broadinstitute.sting.utils.recalibration.BaseRecalibration.(BaseRecalibration.java:74) at org.broadinstitute.sting.gatk.GenomeAnalysisEngine.setBaseRecalibration(GenomeAnalysisEngine.java:217) at org.broadinstitute.sting.gatk.GenomeAnalysisEngine.execute(GenomeAnalysisEngine.java:253) at org.broadinstitute.sting.gatk.CommandLineExecutable.execute(CommandLineExecutable.java:113) at org.broadinstitute.sting.commandline.CommandLineProgram.start(CommandLineProgram.java:237) at org.broadinstitute.sting.commandline.CommandLineProgram.start(CommandLineProgram.java:147) at org.broadinstitute.sting.gatk.CommandLineGATK.main(CommandLineGATK.java:91)

ERROR ------------------------------------------------------------------------------------------
ERROR A GATK RUNTIME ERROR has occurred (version 2.3-9-ge5ebf34):
ERROR
ERROR Please visit the wiki to see if this is a known problem
ERROR If not, please post the error, with stack trace, to the GATK forum
ERROR Visit our website and forum for extensive documentation and answers to
ERROR commonly asked questions http://www.broadinstitute.org/gatk
ERROR
ERROR MESSAGE: Key 1036 is too large for dimension 2 (max is 1001)
ERROR ------------------------------------------------------------------------------------------

---------------------------------------------------------------EXECUTION STEP---------------------------------------------------------

2 MARKING PCR DUPLICATE

java -Xmx4g -Djava.io.tmpdir=/tmp -jar MarkDuplicates.jar INPUT=M9.bam OUTPUT=m9.marked.bam METRICS_FILE=metrics CREATE_INDEX=true VALIDATION_STRINGENCY=LENIENT

3 LOCAL REALIGNMENT AROUND INDEL

3.1

java -Xmx4g -jar GenomeAnalysisTK.jar -T RealignerTargetCreator -R ucsc.hg19.fasta -knowndbsnp_137.hg19.vcf -o m9.list -I m9.marked.bam

3.2

java -Xmx4g -Djava.io.tmpdir=/tmp -jar GenomeAnalysisTK.jar -I m9.marked.bam -R ucsc.hg19.fasta -T IndelRealigner -targetIntervals m9.list -known dbsnp_137.hg19.vcf -o m9.marked.realigned.bam

3.2.1

java -Xmx4g -jar GenomeAnalysisTK.jar -R ucsc.hg19.fasta -T ReduceReads -I m9.marked.realigned.bam -o m9.marked.realigned.reduce.bam

3.3

java -Djava.io.tmpdir=/tmp/flx-auswerter -Xmx4g -jar FixMateInformation.jar INPUT=m9.marked.realigned.reduce.bam OUTPUT=m9.marked.realigned.reduce.fixed.bam SO=coordinate VALIDATION_STRINGENCY=LENIENT CREATE_INDEX=true

4 QUALITY SCORE RECALIBRATION

4.1

java -Xmx4g -jar GenomeAnalysisTK.jar -l INFO -R ucsc.hg19.fasta -knownSites dbsnp_137.hg19.vcf -I m9.marked.realigned.reduce.fixed.bam -T BaseRecalibrator -maxCycle 1500 -cov ReadGroupCovariate -cov QualityScoreCovariate -o m9.recal_data.grp

4.2 ***********************************

java -Xmx4g -jar GenomeAnalysisTK.jar -T PrintReads -R ucsc.hg19.fasta -I m9.marked.realigned.reduce.fixed.bam -BQSR m9.recal_data.grp -o m9.marked.realigned.reduce.fixed.recal.bam

5 SNP CALLING

5.1

java -Xmx4g -jar GenomeAnalysisTK.jar -nct 4 --num_threads 4 -glm BOTH -R ucsc.hg19.fasta -T UnifiedGenotyper --sample_ploidy 5 -I m9.marked.realigned.reduce.fixed.bam -D dbsnp_137.hg19.vcf -o m9.vcf -stand_call_conf 20.0 -stand_emit_conf 20.0
-A DepthOfCoverage -A AlleleBalance

Comments (1)

HI When I run Base recabrator with the following command:

java -Xmx4g -jar /usr/bin/GenomeAnalysisTK.jar -T BaseRecalibrator -I realignedBam.bam  -R /data1/human_g1k_v37.fasta --knownSites /data1/snp132.vcf -o recalibration_report.grp

I get the following error :

INFO  07:15:53,380 HttpMethodDirector - I/O exception (javax.net.ssl.SSLException) caught when processing request: Unrecognized SSL message, plaintext connection? 
INFO  07:15:53,380 HttpMethodDirector - Retrying request 
INFO  07:15:53,386 HttpMethodDirector - I/O exception (javax.net.ssl.SSLException) caught when processing request: Unrecognized SSL message, plaintext connection? 
INFO  07:15:53,387 HttpMethodDirector - Retrying request 
INFO  07:15:53,393 HttpMethodDirector - I/O exception (javax.net.ssl.SSLException) caught when processing request: Unrecognized SSL message, plaintext connection? 
INFO  07:15:53,393 HttpMethodDirector - Retrying request 
INFO  07:15:53,398 HttpMethodDirector - I/O exception (javax.net.ssl.SSLException) caught when processing request: Unrecognized SSL message, plaintext connection? 
INFO  07:15:53,398 HttpMethodDirector - Retrying request 
INFO  07:15:53,405 HttpMethodDirector - I/O exception (javax.net.ssl.SSLException) caught when processing request: Unrecognized SSL message, plaintext connection? 
INFO  07:15:53,405 HttpMethodDirector - Retrying request 
##### ERROR ------------------------------------------------------------------------------------------
##### ERROR A USER ERROR has occurred (version 2.0-34-g07bda93): 
##### 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: Invalid command line: No tribble type was provided on the command line and the type of the file could not be determined dynamically. Please add an explicit type tag :NAME listing the correct type from among the supported types:
##### ERROR          Name        FeatureType   Documentation
##### ERROR          BCF2     VariantContext   http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_utils_codecs_bcf2_BCF2Codec.html
##### ERROR        BEAGLE      BeagleFeature   http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_utils_codecs_beagle_BeagleCodec.html
##### ERROR           BED         BEDFeature   http://www.broadinstitute.org/gatk/gatkdocs/org_broad_tribble_bed_BEDCodec.html
##### ERROR      BEDTABLE       TableFeature   http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_utils_codecs_table_BedTableCodec.html
##### ERROR EXAMPLEBINARY            Feature   http://www.broadinstitute.org/gatk/gatkdocs/org_broad_tribble_example_ExampleBinaryCodec.html
##### ERROR      GELITEXT    GeliTextFeature   http://www.broadinstitute.org/gatk/gatkdocs/org_broad_tribble_gelitext_GeliTextCodec.html
##### ERROR      OLDDBSNP    OldDbSNPFeature   http://www.broadinstitute.org/gatk/gatkdocs/org_broad_tribble_dbsnp_OldDbSNPCodec.html
##### ERROR     RAWHAPMAP   RawHapMapFeature   http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_utils_codecs_hapmap_RawHapMapCodec.html
##### ERROR        REFSEQ      RefSeqFeature   http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_utils_codecs_refseq_RefSeqCodec.html
##### ERROR     SAMPILEUP   SAMPileupFeature   http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_utils_codecs_sampileup_SAMPileupCodec.html
##### ERROR       SAMREAD     SAMReadFeature   http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_utils_codecs_samread_SAMReadCodec.html
##### ERROR         TABLE       TableFeature   http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_utils_codecs_table_TableCodec.html
##### ERROR           VCF     VariantContext   http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_utils_codecs_vcf_VCFCodec.html
##### ERROR          VCF3     VariantContext   http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_utils_codecs_vcf_VCF3Codec.html
##### ERROR ------------------------------------------------------------------------------------------
Comments (1)

Hi I´ve a strange problem with the GATK. Everytime I try to run it my Console shows the following error Message.

30.10.12 15:10:06   [0x0-0xe20e2].com.apple.JarLauncher[1114]   ##### ERROR ------------------------------------------------------------------------------------------
30.10.12 15:10:06   [0x0-0xe20e2].com.apple.JarLauncher[1114]   ##### ERROR A USER ERROR has occurred (version 2.1-13-g0f021e6): 
30.10.12 15:10:06   [0x0-0xe20e2].com.apple.JarLauncher[1114]   ##### ERROR The invalid arguments or inputs must be corrected before the GATK can proceed
30.10.12 15:10:06   [0x0-0xe20e2].com.apple.JarLauncher[1114]   ##### ERROR Please do not post this error to the GATK forum
30.10.12 15:10:06   [0x0-0xe20e2].com.apple.JarLauncher[1114]   ##### ERROR
30.10.12 15:10:06   [0x0-0xe20e2].com.apple.JarLauncher[1114]   ##### ERROR See the documentation (rerun with -h) for this tool to view allowable command-line arguments.
30.10.12 15:10:06   [0x0-0xe20e2].com.apple.JarLauncher[1114]   ##### ERROR Visit our website and forum for extensive documentation and answers to 
30.10.12 15:10:06   [0x0-0xe20e2].com.apple.JarLauncher[1114]   ##### ERROR commonly asked questions http://www.broadinstitute.org/gatk
30.10.12 15:10:06   [0x0-0xe20e2].com.apple.JarLauncher[1114]   ##### ERROR
30.10.12 15:10:06   [0x0-0xe20e2].com.apple.JarLauncher[1114]   ##### ERROR MESSAGE: Argument with name '--analysis_type' (-T) is missing.
30.10.12 15:10:06   [0x0-0xe20e2].com.apple.JarLauncher[1114]   ##### ERROR ------------------------------------------------------------------------------------------

Can you show me my mistakes please? With regards Oliver

Comments (2)

Hello,

During running of the depthOfCoverage tool, I get the error: /tmp/RsQHCt1W: No space left on device

I have tried changing the TMPDIR environment variable (and exporting) but eventually I get the same error. Is there a way to change the temporary directory that GATK uses?

I'm running GATK v2.1-8-g5efb575 on a Linux system.

Thanks, Rick