Test-drive the GATK tools and Best Practices pipelines on Terra


Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.

CombineGVCFs : ERROR MESSAGE: Line 116400953: there aren't enough columns for line

Dear the GATK team,

There are a problem when running CombineGVCF:

The ERROR MESSAGE is : Line 116400953: there aren't enough columns for line.
But I check the 116400944-116400954 line of input files, it seems not having something wrong.

chr3 188853626 . T . . END=188853636 GT:DP:GQ:MIN_DP:PL 0/0:13:33:11:0,33,374
chr3 188853637 . C . . END=188853640 GT:DP:GQ:MIN_DP:PL 0/0:11:30:11:0,30,450
chr3 188853641 . T . . END=188853649 GT:DP:GQ:MIN_DP:PL 0/0:11:27:10:0,27,405
chr3 188853650 . C . . END=188853654 GT:DP:GQ:MIN_DP:PL 0/0:10:24:10:0,24,360
chr3 188853655 . G . . END=188853664 GT:DP:GQ:MIN_DP:PL 0/0:9:21:8:0,21,315
chr3 188853665 . G . . END=188853666 GT:DP:GQ:MIN_DP:PL 0/0:7:15:7:0,15,225
chr3 188853667 . A . . END=188853680 GT:DP:GQ:MIN_DP:PL 0/0:7:12:4:0,12,140
chr3 188853681 . T . . END=188853688 GT:DP:GQ:MIN_DP:PL 0/0:5:15:5:0,15,174
chr3 188853689 . C . . END=188853691 GT:DP:GQ:MIN_DP:PL 0/0:5:9:5:0,9,135
chr3 188853692 . A . . END=188853699 GT:DP:GQ:MIN_DP:PL 0/0:6:12:4:0,12,135

The commend line is :

/ifshk5/PC_PA_EU/USER/zhangbaifeng/software/java/jre1.8.0_45/bin/java -Xmx4g -Djava.io.tmpdir=/ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/java_tmp -jar /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/software/GATK/GenomeAnalysisTK.jar -et NO_ET -K /ifshk5/PC_PA_EU/USER/zhangbaifeng/software/gatk/zhangbaifeng_genomics.cn.key -T CombineGVCFs -R /ifshk1/BC_CANCER/01bin/DNA/software/pipeline/CSAP_v5.2.4/Database/human_19/hg19_fasta_GATK/hg19.fasta \
-V /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/haplotypecaller/vcf/M54_001_FA.recal.raw_variants.g.vcf \
-V /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/haplotypecaller/vcf/M54_001_MA.recal.raw_variants.g.vcf \
-V /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/haplotypecaller/vcf/M54_002_MA.recal.raw_variants.g.vcf \
-V /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/haplotypecaller/vcf/M54_003_MA.recal.raw_variants.g.vcf \
-V /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/haplotypecaller/vcf/M54_004_FA.recal.raw_variants.g.vcf \
-V /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/haplotypecaller/vcf/M54_006_FA.recal.raw_variants.g.vcf \
-V /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/haplotypecaller/vcf/M54_006_MA.recal.raw_variants.g.vcf \
-V /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/haplotypecaller/vcf/M54_007_MA.recal.raw_variants.g.vcf \
-V /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/haplotypecaller/vcf/M54_008_FA.recal.raw_variants.g.vcf \
-V /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/haplotypecaller/vcf/M54_008_MA.recal.raw_variants.g.vcf \
-V /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/haplotypecaller/vcf/M54_009_MA.recal.raw_variants.g.vcf \
-V /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/haplotypecaller/vcf/M54_010_MA.recal.raw_variants.g.vcf \
-V /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/haplotypecaller/vcf/M54_012_FA.recal.raw_variants.g.vcf \
-V /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/haplotypecaller/vcf/M54_014_FA.recal.raw_variants.g.vcf \
-V /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/haplotypecaller/vcf/M54_015_FA.recal.raw_variants.g.vcf \
-V /ifshk7/BC_RES/TECH/PMO/zhangbaifeng/330.snp.analysis/330.sample.GVCF/haplotypecaller/vcf/M54_016_FA.recal.raw_variants.g.vcf \
-o M54.g.vcf.gz

Thanks.

Tagged:

Answers

Sign In or Register to comment.