GATK_UnifiedGenotyper_Unable to merge temporary Tribble output file

douymdouym Posts: 18Member
edited January 2013 in Ask the GATK team

Hi all,

I've been analyzing some illumina whole exome sequencing data these days. Yesterday I used GATK(version 2.0) UnifiedGenotyper to call snps and indels with the following commands: -T UnifiedGenotyper -R GRCh37/human_g1k_v37.fasta -I GATK_recal_result.bam -glm BOTH --dbsnp reference/dbsnp_134.b37.vcf -stand_call_conf 50 -stand_emit_conf 10 -o raw2.vcf -dcov 200 --num_threads 10

After running theses commands, I got a vcf file which is very small(when I checked the vcf file, I found these called snps and indels are all from Chromosome1) The error message is as follows:

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

org.broadinstitute.sting.utils.exceptions.ReviewedStingException: Unable to merge temporary Tribble output file. at org.broadinstitute.sting.gatk.executive.HierarchicalMicroScheduler.mergeExistingOutput( at org.broadinstitute.sting.gatk.executive.HierarchicalMicroScheduler.execute( at org.broadinstitute.sting.gatk.GenomeAnalysisEngine.execute( at org.broadinstitute.sting.gatk.CommandLineExecutable.execute( at org.broadinstitute.sting.commandline.CommandLineProgram.start( at org.broadinstitute.sting.commandline.CommandLineProgram.start( at org.broadinstitute.sting.gatk.CommandLineGATK.main( Caused by: org.broad.tribble.TribbleException$MalformedFeatureFile: Unable to parse header with error: /rd/tmp/ (Too many open files), for input source: /rd/tmp/ at org.broad.tribble.TribbleIndexedFeatureReader.readHeader( at org.broad.tribble.TribbleIndexedFeatureReader.( at org.broad.tribble.AbstractFeatureReader.getFeatureReader( at at at org.broadinstitute.sting.gatk.executive.OutputMergeTask.merge( at org.broadinstitute.sting.gatk.executive.HierarchicalMicroScheduler.mergeExistingOutput( ... 6 more Caused by: /rd/tmp/ (Too many open files) at Method) at at org.broad.tribble.util.ParsingUtils.openInputStream( at org.broad.tribble.TribbleIndexedFeatureReader.readHeader( ... 12 more

ERROR ------------------------------------------------------------------------------------------
ERROR A GATK RUNTIME ERROR has occurred (version 2.0-39-gd091f72):
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
ERROR MESSAGE: Unable to merge temporary Tribble output file.
ERROR ------------------------------------------------------------------------------------------

Would you please help me solve it ? Thanks a lot

Post edited by Geraldine_VdAuwera on

Best Answer


Sign In or Register to comment.