mutect2 error in strand bias count?

I have some examples showing something wrong in calculating ALT_F1R2 and ALT_F2R1 in MuTect2. In this example, the alignment in the locally re-assemblied BAM (by MuTect2) shows that all 12 reads are aligned one directional (reverse) but ALT_F1R2 and ALT_F2R1 counts in the vcf file are 9 and 3, respectively. This SNP is actually rejected by MuTect 1 due to 'strand_artifact' which looks right but passed in MuTect2. Is this a bug in the MuTect2?

ECNT=1;HCNT=8;MAX_ED=.;MIN_ED=.;NLOD=9.93;TLOD=27.01 GT:AD:AF:ALT_F1R2:ALT_F2R1:FOXOG:QSS:REF_F1R2:REF_F2R1
0/1:83,12:0.141:9:3:0.250:2695,364:55:28
0/0:35,0:0.00:0:0:.:1120,0:29:6

Best Answer

Answers

Sign In or Register to comment.