Thanks, much appreciated!
I tried with . . . -I tumor.bam -I normal.bam -normal normal_sample_name . . .
No luck, I still variants in found in tumor and Normal :(
Any idea how to get VCF outfile showing only tumor variants (filtering out Normal, not showing normal_sample_name in output vcf file)?
• 0 views
•
link
updated
by
Ram
•
written
by
m. bioinfo •
•
written
by
m. bioinfo •
HI Cocchi,
Although I haven't gone through deeply, I believe you should check out this link where the GATK group described the usage of Mutect2 for a paired t-n tissues. It's updated 14 days ago. Not sure if it serves your WGS data though. https://gatk.broadinstitute.org/hc/en-us/articles/360035894731-Somatic-short-variant-discovery-SNVs-Indels-
Field
Hi , I have used Mutect2 with tumor normal, but Why my Mutect2 output VCF shows variants found in both (Normal and Tumor) samples ?
It should only contains/show variants in Tumor sample and Normal variants should be filtered out, I am interested in variants found in Tumor (not in matched normal).
TIA
Your VCF file contains the genotype of both samples, it's your choice why,how and when to filter variants.
Thanks ! If I understand correctly, I have to filter out variants found in Normal. Would GATK FilterMutectCalls or SelectVariants be a choice ?
I read the cmd arguments but there is nothing like to filter out variants found in normal. These two can add PASS filter and keep PASS only. Any point of direction would be appreciated.
If you give both Normal-Tumor BAM file, your VCF will contain two ending columns, one for each sample. With
FilterMutectCallsyou'll apply a filtering over the previous called variants using a mixture Bayesian model, with several tunable parameters. Apply a filter means that yourFILTERfield in the VCF file will be filled with a reason (if any). With SelectVariants you could product a VCF excluding variants with specific/any reason in theFILTERfield. So typically you'd use both in the explained order. Variants found in normal file will be flagged with theGERMLINEtag, and excluded by SelectVariants.Thank you ! I tried with FilterMutectCalls but there is no Filter flag worked to remove normal sample name from VCF. In fact, no Filter flag is available for the same.
However, SelectVariant has worked for me to exclude normal_sample_name column as well as normal variants from VCF file. It resulted tumor only column showing tumor variants.
The commands worked for me is :
gatk SelectVariant -V Mutect2_Filtered.vcf.gz -R ref.fa --exclude-filtered true --exclude-sample-name normal_sample_name -O Mutect2_Filtered_PASS.vcf.gz
What I have learned that :
--exclude-filtered true # keeps PASS variants only
--exclude-sample-name normal_sample_name # Removes normal sample column and variants found in normal sample.
Do I need any corrections here ?
I still think you're misunderstanding how a VCF file is designed. As I've stated before,
FilterMutectCallshas the role of flagging variants, trying to identify sources of contamination, strand/orientation bias, Germline variants, and others. When you're calling variants in tumor sample using a matched normal sample,Mutect2will use the tumor sample as a source of variants lookup. In the same locus the search would be conducted also on normal sample, and the identified genotype will be reported in the VCF. So I'm not sure what you are referring viavariants found in normal sample. Germlines? Contamination?The flag you're using are just exclusions based on the VCF file itself, nothing that couldn't be done also via an awk command in Bash. The
--exclude-sample-nameremoves just the column containing the genotype of the passed sample name.Yes, You are right, I had typos above which I corrected. and yes, my confusion as well as you pointed out.
Actually, after running Mutect2 with Tumor and Normal, I want to keep variants found only in Tumor . For this reason, I have run:
FilterMutectCalls
SelectVariant : --exclude-filtered true (keep PASS only) and --exclude-sample-name normal_sample_name (to exclude normal variants).
Do you think, I would require to do some additional filtering ? or any corrections here please ?
Variants with no flags (PASS only) are considered somatic for the tumor sample in your setup. As it is, it's enough.
Thanks, PASS only variants are actually shown in both tumor and normal within single VCF from Mutect2 run . I want tumor somatic mutations filtering out normal PASS variants.
In other words, I get somatic PASS variants in tumor and normal from Mutec2 output. :(