To celebrate the release of GATK 4.0, we are giving away free credits for running the GATK4 Best Practices pipelines in FireCloud, our secure online analysis portal. It’s first come first serve, so sign up now to claim your free credits worth $250. Sponsored by Google Cloud. Learn more at https://software.broadinstitute.org/firecloud/documentation/freecredits

MuTect2 output error

Hi
I am facing an error in the MuTect2 output for a couple of my tumor-normal pair samples. The output shows a single merged column for tumor and normal info, that is, I am getting only 1 column for variant info instead of 2 separate ones for tumor and normal. The command I used is as follows:

java -Xmx6g -jar GenomeAnalysisTK.jar -T MuTect2 -R EXOME_REF/exome_ref.fa -I:tumor 15D_recal.bam -I:normal 15N_recal.bam --dbsnp dbsnp_138.hg19.vcf --cosmic CosmicCodingMuts.vcf --cosmic CosmicNonCodingVariants.vcf -o s15_mutect.vcf -U ALLOW_SEQ_DICT_INCOMPATIBILITY

First few lines from my vcf output (excluding the headers) is as follows:

CHROM POS ID REF ALT QUAL FILTER INFO FORMAT 15D

chr10 95188 rs201207808 G A . alt_allele_in_normal DB;ECNT=1;HCNT=2;MAX_ED=.;MIN_ED=.;NLOD=11.53;TLOD=16.48 GT:AD:AF:ALT_F1R2:ALT_F2R1:FOXOG:QSS:REF_F1R2:REF_F2R1 0/0:181,6:0.047:5:1:0.833:6051,215:97:84
chr10 300015 . G T . alt_allele_in_normal;t_lod_fstar ECNT=1;HCNT=2;MAX_ED=.;MIN_ED=.;NLOD=3.48;TLOD=4.28 GT:AD:AF:ALT_F1R2:ALT_F2R1:FOXOG:QSS:REF_F1R2:REF_F2R1 0/0:36,2:0.054:2:0:1.00:1271,67:24:12
chr10 532531 . C A . alt_allele_in_normal;t_lod_fstar ECNT=1;HCNT=8;MAX_ED=.;MIN_ED=.;NLOD=13.91;TLOD=4.08 GT:AD:AF:ALT_F1R2:ALT_F2R1:FOXOG:QSS:REF_F1R2:REF_F2R1 0/0:77,2:0.027:0:2:1.00:2452,71:29:48
chr10 1038425 . G T . alt_allele_in_normal;t_lod_fstar ECNT=1;HCNT=2;MAX_ED=.;MIN_ED=.;NLOD=7.75;TLOD=4.18 GT:AD:AF:ALT_F1R2:ALT_F2R1:FOXOG:QSS:REF_F1R2:REF_F2R1 0/0:63,3:0.038:0:3:0.00:2289,102:39:24
chr10 1089949 . G T . alt_allele_in_normal;t_lod_fstar ECNT=1;HCNT=4;MAX_ED=.;MIN_ED=.;NLOD=7.35;TLOD=4.29 GT:AD:AF:ALT_F1R2:ALT_F2R1:FOXOG:QSS:REF_F1R2:REF_F2R1 0/0:53,3:0.039:1:2:0.333:1922,104:34:19
chr10 1209599 . CTGTG C . alt_allele_in_normal;t_lod_fstar ECNT=1;HCNT=4;MAX_ED=.;MIN_ED=.;NLOD=15.60;RPA=9,7;RU=TG;STR;TLOD=4.71 GT:AD:AF:ALT_F1R2:ALT_F2R1:FOXOG:QSS:REF_F1R2:REF_F2R1 0/0:140,3:0.010:0:3:.:799,110:5:18
chr10 3182869 . C A . alt_allele_in_normal;t_lod_fstar ECNT=1;HCNT=1;MAX_ED=.;MIN_ED=.;NLOD=11.94;TLOD=4.15 GT:AD:AF:ALT_F1R2:ALT_F2R1:FOXOG:QSS:REF_F1R2:REF_F2R1 0/0:108,9:0.096:6:3:0.333:2824,115:49:59
chr10 5471029 . G T . alt_allele_in_normal;t_lod_fstar ECNT=1;HCNT=1;MAX_ED=.;MIN_ED=.;NLOD=4.58;TLOD=4.01 GT:AD:AF:ALT_F1R2:ALT_F2R1:FOXOG:QSS:REF_F1R2:REF_F2R1 0/0:38,2:0.050:1:1:0.500:1320,65:26:12
chr10 5726777 . C T . alt_allele_in_normal ECNT=1;HCNT=20;MAX_ED=.;MIN_ED=.;NLOD=15.83;TLOD=6.88 GT:AD:AF:ALT_F1R2:ALT_F2R1:FOXOG:QSS:REF_F1R2:REF_F2R1 0/0:159,13:0.094:7:6:0.462:4861,202:63:96
chr10 7614076 . CTG C . alt_allele_in_normal;str_contraction;t_lod_fstar ECNT=1;HCNT=2;MAX_ED=.;MIN_ED=.;NLOD=7.24;RPA=6,5;RU=TG;STR;TLOD=5.79 GT:AD:AF:ALT_F1R2:ALT_F2R1:FOXOG:QSS:REF_F1R2:REF_F2R1 0/0:75,4:0.049:2:2:.:2662,147:39:36
chr10 7774309 . G T . alt_allele_in_normal;t_lod_fstar ECNT=1;HCNT=2;MAX_ED=.;MIN_ED=.;NLOD=6.39;TLOD=4.17 GT:AD:AF:ALT_F1R2:ALT_F2R1:FOXOG:QSS:REF_F1R2:REF_F2R1 0/0:56,2:0.043:0:2:0.00:1823,68:30:26

Kindly suggest me what should I do to rectify the error.

Sign In or Register to comment.