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.