Hi all
I've tried snp calling using GATK unifiedgenotyper but got the error message - it seems that there's a problem in fasta sequence dictionary. I am currently following <Exome analysis> in How-to Wiki section (http://seqanswers.com/wiki/How-to/ex..._recalibration) However, it seems that some of procedures does not use the latest program. So, I changed few things for those.
Please give any advice, it will be really appreciate.
Here is command
and error message
I've tried snp calling using GATK unifiedgenotyper but got the error message - it seems that there's a problem in fasta sequence dictionary. I am currently following <Exome analysis> in How-to Wiki section (http://seqanswers.com/wiki/How-to/ex..._recalibration) However, it seems that some of procedures does not use the latest program. So, I changed few things for those.
Please give any advice, it will be really appreciate.
Here is command
$ java -Xmx4g -jar /usr/local/bin/gatk/GenomeAnalysisTK-1.4-30-gf2ef8d1/GenomeAnalysisTK.jar -glm BOTH -R hg19.fa -T UnifiedGenotyper -I input.marked.realigned.fixed.recal.bam -D dbsnp132.txt -o snps.vcf -metrics snps.metrics -stand_call_conf 50.0 -stand_emit_conf 10.0 -dcov 1000 -A DepthOfCoverage -A AlleleBalance -L target_intervals.bed
INFO 20:44:34,610 RodBindingArgumentTypeDescriptor - Dynamically determined type of target_intervals.bed to be BED
INFO 20:44:34,647 HelpFormatter - ---------------------------------------------------------------------------------
INFO 20:44:34,648 HelpFormatter - The Genome Analysis Toolkit (GATK) v1.4-30-gf2ef8d1, Compiled 2012/02/17 20:18:04
INFO 20:44:34,648 HelpFormatter - Copyright (c) 2010 The Broad Institute
INFO 20:44:34,648 HelpFormatter - Please view our documentation at http://www.broadinstitute.org/gsa/wiki
INFO 20:44:34,648 HelpFormatter - For support, please view our support site at http://getsatisfaction.com/gsa
INFO 20:44:34,649 HelpFormatter - Program Args: -glm BOTH -R hg19.fa -T UnifiedGenotyper -I input.marked.realigned.fixed.recal.bam -D dbsnp132.txt -o snps.vcf -metrics snps.metrics -stand_call_conf 50.0 -stand_emit_conf 10.0 -dcov 1000 -A DepthOfCoverage -A AlleleBalance -L target_intervals.bed
INFO 20:44:34,649 HelpFormatter - Date/Time: 2012/02/24 20:44:34
INFO 20:44:34,649 HelpFormatter - ---------------------------------------------------------------------------------
INFO 20:44:34,649 HelpFormatter - ---------------------------------------------------------------------------------
INFO 20:44:34,669 RodBindingArgumentTypeDescriptor - Dynamically determined type of dbsnp132.txt to be VCF
INFO 20:44:34,682 GenomeAnalysisEngine - Strictness is SILENT
INFO 20:44:34,737 SAMDataSource$SAMReaders - Initializing SAMRecords in serial
INFO 20:44:34,752 SAMDataSource$SAMReaders - Done initializing BAM readers: total time 0.01
INFO 20:44:34,762 RMDTrackBuilder - Loading Tribble index from disk for file dbsnp132.txt
INFO 20:44:38,171 GATKRunReport - Uploaded run statistics report to AWS S3
##### ERROR ------------------------------------------------------------------------------------------
##### ERROR A USER ERROR has occurred (version 1.4-30-gf2ef8d1):
##### 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: Contig chr6_apd_hap1 given as location, but this contig isn't present in the Fasta sequence dictionary
##### ERROR ------------------------------------------------------------------------------------------
INFO 20:44:34,647 HelpFormatter - ---------------------------------------------------------------------------------
INFO 20:44:34,648 HelpFormatter - The Genome Analysis Toolkit (GATK) v1.4-30-gf2ef8d1, Compiled 2012/02/17 20:18:04
INFO 20:44:34,648 HelpFormatter - Copyright (c) 2010 The Broad Institute
INFO 20:44:34,648 HelpFormatter - Please view our documentation at http://www.broadinstitute.org/gsa/wiki
INFO 20:44:34,648 HelpFormatter - For support, please view our support site at http://getsatisfaction.com/gsa
INFO 20:44:34,649 HelpFormatter - Program Args: -glm BOTH -R hg19.fa -T UnifiedGenotyper -I input.marked.realigned.fixed.recal.bam -D dbsnp132.txt -o snps.vcf -metrics snps.metrics -stand_call_conf 50.0 -stand_emit_conf 10.0 -dcov 1000 -A DepthOfCoverage -A AlleleBalance -L target_intervals.bed
INFO 20:44:34,649 HelpFormatter - Date/Time: 2012/02/24 20:44:34
INFO 20:44:34,649 HelpFormatter - ---------------------------------------------------------------------------------
INFO 20:44:34,649 HelpFormatter - ---------------------------------------------------------------------------------
INFO 20:44:34,669 RodBindingArgumentTypeDescriptor - Dynamically determined type of dbsnp132.txt to be VCF
INFO 20:44:34,682 GenomeAnalysisEngine - Strictness is SILENT
INFO 20:44:34,737 SAMDataSource$SAMReaders - Initializing SAMRecords in serial
INFO 20:44:34,752 SAMDataSource$SAMReaders - Done initializing BAM readers: total time 0.01
INFO 20:44:34,762 RMDTrackBuilder - Loading Tribble index from disk for file dbsnp132.txt
INFO 20:44:38,171 GATKRunReport - Uploaded run statistics report to AWS S3
##### ERROR ------------------------------------------------------------------------------------------
##### ERROR A USER ERROR has occurred (version 1.4-30-gf2ef8d1):
##### 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: Contig chr6_apd_hap1 given as location, but this contig isn't present in the Fasta sequence dictionary
##### ERROR ------------------------------------------------------------------------------------------
Comment