Entering edit mode
6.9 years ago
anu014
▴
190
Greeting of the day Biostars!
I'm dealing with variant calling in WES data nowadays. For this, I'm using GATK ,most trusted tool for variant calling. But how to know which is the true variant & which one is false positive? For this I learned there's something called 'VQSRTranche' filters & it goes like this :
FILTER=<ID=VQSRTrancheSNP99.00to99.90,Description="Truth sensitivity tranche level for SNP model at VQS Lod: -1.3154 <= x < 0.5827">
So following are the questions I've on my mind :
i. What is tranche level/tranche?
ii. Eg. VQSRTrancheSNP99.00to99.90 means the variant falls in confidence level 99.00 to 99.9 WITH false positive % being 0.1?
iii. 'VQS Lod' value for same trance changes from data to data.
iv. '--ts_filter_level 99' marks tranche level below 99 as PASS/TRUE variant?
read the doc: https://gatkforums.broadinstitute.org/gatk/discussion/39/variant-quality-score-recalibration-vqsr