This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Why Mutect2 calling only germline SNV from normal-tumor datasets ?

Hi I called variants using mutect2 from normal and tumor samples. I'm surprised it showed all germ-line variants cant find single somatic variant . Whether I missed any parameters ? how to sort out this problem.

Thanks *Input *

Cancer type - breast
Normal bam size - 67 GB
Tumor bam size - 77 GB
Ref- hg38.fa

CMD-

java -Xmx32G -jar gatk-package-4.0.2.1-local.jar Mutect2 -I /data/cancer_samples/sample1/SAM/markduplicatesfixed_Realigned.normal.bam -normal NA12878 -O sample1_mutect2_normaltumor.vcf -R /data/Ref/hg38.fa -I /data/cancer_samples/sample1/SAM//markduplicatesfixed_Realigned.tumor.bam -tumor NA12878

Output-
sample1_mutect2_normaltumor.vcf - 25 mb

**VCF output **
CHROM POS ID REF ALT QUAL FILTER INFO FORMAT NA12878

chr1 16206 . T A . . DP=99;ECNT=2;NLOD=5.63;N_ART_LOD=10.07;POP_AF=1.000e-03;P_GERMLINE=-2.329e+00;TLOD=10.07 GT:AD:AF:F1R2:F2R1:MBQ:MFRL:MMQ:MPOS:SA_MAP_AF:SA_POST_PROB 0/0:88,8:0.100:52,6:36,2:40:321,344:22:25
chr1 16263 . G C . . DP=102;ECNT=2;NLOD=15.02;N_ART_LOD=3.98;POP_AF=1.000e-03;P_GERMLINE=-1.374e+01;TLOD=3.98 GT:AD:AF:F1R2:F2R1:MBQ:MFRL:MMQ:MPOS:SA_MAP_AF:SA_POST_PROB 0/0:94,6:0.071:49,3:45,3:39:318,317:20:44
chr1 55249 . C CTATGG . . DP=102;ECNT=1;NLOD=5.68;N_ART_LOD=12.14;POP_AF=1.000e-03;P_GERMLINE=-2.380e+00;TLOD=12.14 GT:AD:AF:F1R2:F2R1:MBQ:MFRL:MMQ:MPOS:SA_MAP_AF:SA_POST_PROB 0/0:86,5:0.102:40,2:46,3:36:321,306:40:14
chr1 77568 . C A . . DP=48;ECNT=1;NLOD=4.58;N_ART_LOD=3.45;POP_AF=1.000e-03;P_GERMLINE=-3.830e+00;TLOD=3.45 GT:AD:AF:F1R2:F2R1:MBQ:MFRL:MMQ:MPOS:SA_MAP_AF:SA_POST_PROB 0/0:40,2:0.105:22,1:18,1:36:325,332:43:3
chr1 84133 . A T . . DP=208;ECNT=1;NLOD=11.73;N_ART_LOD=25.43;POP_AF=1.000e-03;P_GERMLINE=-8.429e+00;TLOD=25.43 GT:AD:AF:F1R2:F2R1:MBQ:MFRL:MMQ:MPOS:SA_MAP_AF:SA_POST_PROB 0/0:181,12:0.094:85,3:96,9:39:320,352:37:17
chr1 91486 . G C . . DP=48;ECNT=2;NLOD=5.44;N_ART_LOD=4.26;POP_AF=1.000e-03;P_GERMLINE=-3.892e+00;TLOD=4.26 GT:AD:AF:F1R2:F2R1:MBQ:MFRL:MMQ:MPOS:PGT:PID:SA_MAP_AF:SA_POST_PROB 0/0:46,2:0.042:17,0:29,2:39:320,321:27:15:0|1:91486_G_C
chr1 91490 . A G . . DP=51;ECNT=2;NLOD=6.05;N_ART_LOD=4.20;POP_AF=1.000e-03;P_GERMLINE=-4.550e+00;TLOD=4.20 GT:AD:AF:F1R2:F2R1:MBQ:MFRL:MMQ:MPOS:PGT:PID:SA_MAP_AF:SA_POST_PROB 0/0:48,2:0.049:19,0:29,2:39:319,321:27:11:0|1:91486_G_C
chr1 98753 . T C . . DP=173;ECNT=1;NLOD=37.48;N_ART_LOD=3.19;POP_AF=1.000e-03;P_GERMLINE=-3.699e+01;TLOD=3.19 GT:AD:AF:F1R2:F2R1:MBQ:MFRL:MMQ:MPOS:SA_MAP_AF:SA_POST_PROB 0/0:159,3:0.050:78,0:81,3:39:326,327:38:11
chr1 98921 . AG A . . DP=117;ECNT=3;NLOD=22.24;N_ART_LOD=4.49;POP_AF=1.000e-03;P_GERMLINE=-2.046e+01;RPA=5,4;RU=G;STR;TLOD=4.49 GT:AD:AF:F1R2:F2R1:MBQ:MFRL:MMQ:MPOS:PGT:PID:SA_MAP_AF:SA_POST_PROB 0/0:112,4:0.037:53,2:59,2:34:320,326:34:24:0|1:98921_AG_A
chr1 98945 . C T . . DP=152;ECNT=3;NLOD=27.49;N_ART_LOD=5.85;POP_AF=1.000e-03;P_GERMLINE=-2.457e+01;TLOD=5.85 GT:AD:AF:F1R2:F2R1:MBQ:MFRL:MMQ:MPOS:PGT:PID:SA_MAP_AF:SA_POST_PROB 0/0:138,4:0.060:68,2:70,2:33:321,324:32:7:0|1:98921_AG_A
chr1 98999 . TTTTATTTA T . . DP=208;ECNT=3;NLOD=13.59;N_ART_LOD=3.50;POP_AF=1.000e-03;
genome sequencing alignment next-gen
• 1,802 views • link
updated by GenoMax GenoMax • written by pinn pinn

Please provide feedback on your previous posts. You have multiple people waiting for your input.

If an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted. You can accept more than one if they work.

I will close this and future question until you take time to revisit older threads. Please let us know if anything is unclear.

• 1 views • link
written by WouterDeCoster WouterDeCoster

Please put some more effort in formatting questions, using markup for increased readability. You can do this by selecting the text and clicking the 101010 button. When you compose or edit a post that button is in your toolbar, see image below:

101010 Button

• 1 views • link
written by WouterDeCoster WouterDeCoster

Hello pinninti1991reddy!

We believe that this post does not fit the main topic of this site.

Feedback required on other threads

For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.

If you disagree please tell us why in a reply below, we'll be happy to talk about it.

Cheers!

• 1 views • link
written by WouterDeCoster WouterDeCoster

0 answers

No answers yet.

Log in to answer this question.

More posts like this