Genome Analysis Toolkit

Variant Discovery in High-Throughput Sequencing Data

GATK process banner

Need Help?

Search our documentation

Community Forum

Hi, How can we help?

Developed in the Data Sciences Platform at the Broad Institute, the toolkit offers a wide variety of tools with a primary focus on variant discovery and genotyping. Its powerful processing engine and high-performance computing features make it capable of taking on projects of any size. Learn more

FilterMutectCalls didn't filtered any variants for my vcf file

Answered
0

3 comments

  • Avatar
    David Benjamin

    Xiya Song The log output simply means that no variants or reads were discarded before filtering.  The way FilterMutectCalls works is that the output vcf is the same as the input vcf, but with a list of filters or "PASS" put in the FILTERS column.  Do you see anything there?

    0
    Comment actions Permalink
  • Avatar
    Xiya Song

    Yes! Thank you. After checking filtered.vcf I found most variants are with "weak_evidence" and some are with "PASS", and also "normal_artifact", "clustered_events"...etc. So the filter did work. So if I understand correctly, I should extract variants with "PASS" by myself if I want to use these filtered variants for downstream analysis?

    0
    Comment actions Permalink
  • Avatar
    David Benjamin

    So if I understand correctly, I should extract variants with "PASS" by myself if I want to use these filtered variants for downstream analysis?

    Xiya Song That's exactly right.

    0
    Comment actions Permalink

Please sign in to leave a comment.

Powered by Zendesk