Question concerning GATK VQSR recalibration · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Question concerning GATK VQSR recalibration

Hi everyone,

I have a question concerning the GATK best practice: after joint genotyping, VQSR recalibration is performed via the resource bundle (e.g. hapmap, 1000G, etc). The resulting vcf file has some variants labelled "PASS" and some are labelled e.g. VQSRTrancheSNP99.90to100.00. As far as I understand those are considered false positives.

Is it fine to only include variants that have been labelled "PASS" in the downstream analysis?

Thanks a lot and best regards!

wgs gatk vqsr wes analysis genome

• 638 views • link

written 3.9 years ago by balaena •

0 answers

No answers yet.

Log in to answer this question.

More posts like this

GATK Mutect2
written 7 months ago by sk_24 •

I am running GATK for the first time and working with RNA-seq data from the MDA-MB-231 cell line grown under two different conditions, with three …
Can I combine variants called using T2T-CHM13 with GRCh38?
written 1.7 years ago by Rina •

I performed alignments and variant calls using GATK Best Practice guidelines for germline variant calls, against both GRCh38 and T2T-CHM13 for my paired-end WES data. …
PON and germline files for Mutect2 Tumor-only mode
written 3.0 years ago by Manuel •

Hi! I am applying Mutect2 in tumor-only.mode having hg19 as reference and samples resulting of targeted sequencing (not only coding gene parts but also some …
Questions concerning GATK
written 3.9 years ago by balaena •

Hi everyone, I wanted to ask some questions concerning medium-sized exome cohort studies with GATK and would be grateful for an answer: 1. g.vcf calling, …
GATK4 - To be known or not to be known, that is the question.
written 4.5 years ago by K.patel5 •

Dear Biostars, This question specifically is about the best practice of using GATK4 to look for snps and indels in human samples. I am fairly …
Padding out a GVCF file with 1000G exomes to get gatk VariantRecalibrator working with a small samp…
written 4.7 years ago by michael.flower.14 •

I've got sequencing data for a small 500 bp amplicon from a few samples. GATK best principles suggest running [VariantRecalibrator][1] on the GVCF files I …
Scatter / Gather for BaseRecalibrator on a single human WES dataset?
written 6.5 years ago by asg •

Hello everyone! Please excuse me if this question is a bit naïve: I'm new to bioinformatics in general and GATK in particular. I am using …
Analysis of exome samples build from different exome kits
written 7.0 years ago by Nicolas Rosewick

Question regarding germline variant calling from exome : I've a dataset composed of ~500 exome samples build using 6 differents kits (the dataset was build …
VQSR with extremely low TPs (<0.01), How low can tranche goes?
written 7.4 years ago by greyman

Hi, Sorry that i have posted this at here because for some unknown reason i am unable post a question in GATK forum. Any help …
Comparing Genome Alignment via gvcf?
written 8.3 years ago by Nick •

Hi all, I've been using a pipeline including novoalign and GATK best practices for alignment, variant calling, joint genotyping, vqsr, etc. I recently gained access …