I have run bam files in GATK Mutect2 in "Tumor with matched normal mode" (Tumor+Normal sample) and also in "Tumor only mode" and I have also run them with and without PoN and germline resource file. Whatever I do, running vcftools afterwards for Site Quality check always returns with the value -1 (Negative 1) for each and every single position. What might I be doing wrong? Please help.
1 answer
vcftools is deprecated. How about using bcftools ?
--site-quality Generates a file containing the per-site SNP quality, as found in the QUAL column of the VCF file
bcfools query -f '%CHROM\t%POS\t%REF\t%ALT\t%QUAL\n' variant call_file.vcf
of course if variant call_file.vcf doesn't have any value for QUAL, the output will just be '.' for the QUAL column...
you don't have any QUAL set in your vcf.
grep -v "^#" variant call_file.vcf | cut -f 6
No, I don't. Running the above command returns a long list with only a dot (".") per line. I have trimmed adapter contents from fastq files using trimmomatic, then aligned with reference genome using bwa mem and then ran some additional commands to remove errors from bam files. Then I used GATK Mutect2 to call variants using tumor and matched normal bam files. The resulting vcf files have no QUAL value. I'm new at this and don't know what I'm doing wrong.
Log in to answer this question.
-
PON and germline files for Mutect2 Tumor-only mode
written by Manuel •
-
Right way to call purecn with cnvkit?
written by ww22runner •
-
GATK MUTECT2
written by Hyeong Seok •
-
Are there any repository for panel of normal for running Mutect2?
written by DareDevil
-
Somatic variant calling with/without matched normal sample (Mutect2)
written by newbio17
-
Somatic variant calling mutect2 tumour only
written by cellulebioinfobiscem •
-
Creating a panel of normal from one normal sample
written by Jordan
-
GATK4 (4.1.2.0) Mutect2 --germline-resource option
written by romagnoli.simone •
-
What is the different if I call germline SNP/Indel with mutect and HaplotypeCaller?
written by Makplus T
-
Can we omit the normal input and compare tumor to pon only in MuTect2?
written by Tania
what is the vcftools command line ?