Attention:
The frontline support team will be unavailable to answer questions on April 15th and 17th 2019. We will be back soon after. Thank you for your patience and we apologize for any inconvenience!

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.