Test-drive the GATK tools and Best Practices pipelines on Terra


Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.

Mutect2 not reporting QUAL values

Hi!
I have some targeted panel data. I am using Mutect2 with the following command:

gatk Mutect2 --af-of-alleles-not-in-resource 0.0000025 --germline-resource af-only-gnomad.raw.sites.h19.vcf.gz -A BaseQualityRankSumTest -A ClippingRankSumTest -A MappingQualityRankSumTest -A QualByDepth -A RMSMappingQuality -A ReadPosRankSumTest -A FisherStrand -A StrandOddsRatio -A InbreedingCoeff -L target.bed -I input.bam -R reference.fasta -bamout bamout.bam -tumor TUMOR --disable-read-filter MateOnSameContigOrNoMappedMateReadFilter -O output.vcf.gz

I would like to perform hard filtering in these samples but I cannot obtain the QD value because QUAL value does not appear. I there anything I am missing?

Best Answer

Answers

Sign In or Register to comment.