Question: How Do I Use Variantfiltration To Filter By Vqslod And/Or The Tranches?
1
gravatar for Jeremy Leipzig
8.3 years ago by
Philadelphia, PA
Jeremy Leipzig19k wrote:

After running VariantRecalibrator and ApplyRecalibration in GATK's VQSR module, I have a VCF with VQSLOD tags, as well as some labels such as "TruthSensitivityTranche99.00to99.90"

How do I use the VariantFiltration to obtain an elite VCF file? What is a reasonable minimum VQSLOD (or should I stick with tranches?)

vcf gatk • 3.6k views
ADD COMMENTlink modified 6 weeks ago by geneart$$40 • written 8.3 years ago by Jeremy Leipzig19k
1
gravatar for Frédéric Bigey
8.3 years ago by
Montpellier, France
Frédéric Bigey290 wrote:

Hi Jeremy,

In our study we used SelectVariant to extract variant position with VQSLOD > 4.0.

$pathToGATK \
  -T SelectVariants \
  -R reference \
  --variant inVCF \
  --select_expressions "VQSLOD >= $lodCutoff" \
  --out outName.vcf

We don't use filter like "TruthSensitivityTranche99.00to99.90".

Fred

ADD COMMENTlink written 8.3 years ago by Frédéric Bigey290
0
gravatar for geneart$$
6 weeks ago by
geneart$$40
United States
geneart$$40 wrote:

I understand it has been a long time since this discussion here, but this has raised a similar question in me.

I have the tranches file generated from VQSR, however based on the tranch plot itself, I would pick 99.9 tranch to go to applyVQSR, as the TP is better in the 99.9 tranch based on the plot ,across cumulative and Tranch specific, TP. However the VQSLOD for that specific tranch is negative. So ,

  1. do I need to address this issue of picking a tranch cutoff based on tranch plot or VQSLOD score before applyVQSR? or do i just go ahead with apply VQSR and then filter based on VQSLOD? Lost here, Any suggestions appreciated!
  2. Also the TiTv is lower than what is suggested by GATK for a WES which is indicative of high FP.

targetTruthSensitivity,numKnown,numNovel,knownTiTv,novelTiTv,minVQSLod,filterName,model,accessibleTruthSites,callsAtTruthSites,truthSensitivity 90.00,1763732,507621,2.1238,1.0777,19.6606,VQSRTrancheSNP0.00to90.00,SNP,957871,862083,0.9000 99.00,2245012,1280942,2.1201,1.0087,4.9358,VQSRTrancheSNP90.00to99.00,SNP,957871,948292,0.9900 99.90,2423827,1416814,2.1208,1.0424,-1.6519,VQSRTrancheSNP99.00to99.90,SNP,957871,956913,0.9990 100.00,2492280,1489960,2.1202,1.0545,-22554.8245,VQSRTrancheSNP99.90to100.00,SNP,957871,957871,1.0000

ADD COMMENTlink written 6 weeks ago by geneart$$40
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1871 users visited in the last hour