Question: Filtering significant somatic variants from a VCF file
gravatar for gprashant17
18 months ago by
gprashant1770 wrote:

I have used GATK's pipeline (RNA-seq alignment) and obtained a VCF file using Mutect2 and FilterMutectCalls. Filters such as PASS, Clustered_events, germline, weak-evidence, etc. were added to the variants. Before annotating the file with dbSNP, COSMIC, ANNOVAR, I would like to filter out significant somatic mutations into a separate file in order to facilitate easier analysis. Is it a good practice to exclude all the variants whose filter is not 'PASS'? Also, in case of deciding whether the variant is germline, is the presence of 'germline' filter alone sufficient or it is better to set a threshold to the GERMQ score?

Or is it better to filter out variants after annotating the files?

mutations rna-seq alignment ngs vcf • 1.6k views
ADD COMMENTlink written 18 months ago by gprashant1770

A lot of the terms you use here don't seem to be global terms. Can you please edit your question and explain what these terms mean? You can look at the ##INFO fields to understand them, except for PASS, which you will need to look at ##FILTER.

  • PASS
  • Clustered_events
  • germline
  • weak-evidence

Also, please add some details on your experiment design - if there were matched normals, panel of normals, etc.

ADD REPLYlink written 18 months ago by _r_am32k

I have aligned RNA-seq tumour samples using STAR and I am following GATK's best practices for variant calling, in order to identify somatic mutations causing cancer. Since I do not have a normal sample, I used Mutect2 with only the required arguments.

I am encountering these filters (PASS, Clustered_events, weak-evidence) for the first time and I am reading the following paper about Mutect filtering for reference:

Apparently, FilterMutectCalls labels variants which are false positives with a list of failed filters and true positives with PASS.

ADD REPLYlink written 18 months ago by gprashant1770
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2470 users visited in the last hour