The current GATK version is 3.6-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Powered by Vanilla. Made with Bootstrap.
Last chance to register for the GATK workshop next week in Basel, Switzerland!

how to increase our call confidence

naarkhoonaarkhoo Posts: 38Member

Our call variants analysis, has more far more variants than we expect, including some low-call-confidences and non-coding variants on whole exome. Is there any suggestion on how to increase our call confidence and decrease the total number of variants. Specifically, I had question about the out_mode option.

below is a sample code of our UnifiedGenotyper,

java -Xmx8g \
-jar /medpop/mpg-psrl/Parabase/Tools/GenomeAnalysisTK-2.3-9-ge5ebf34/GenomeAnalysisTK.jar \
-T UnifiedGenotyper \
-R /medpop/mpg-psrl/Parabase/hg19/ucsc.hg19.fasta \
-I /medpop/mpg-psrl/Parabase/Projects/final/final.list \
-o /medpop/mpg-psrl/Parabase/Projects/final/S284.snp.small.vcf \
-stand_call_conf 50.0 \
-stand_emit_conf 10.0 \
--dbsnp /medpop/mpg-psrl/Parabase/hg19/dbsnp_137.hg19.vcf \
-l INFO \
-A DepthOfCoverage \
-A HaplotypeScore \
-A InbreedingCoeff \
-A LowMQ \
-glm SNP \
-dcov 200 \
-nt 4 \
-nct 3 \



  • CarneiroCarneiro Posts: 274Administrator, Dev admin

    are they passing filters? your emit_conf is pretty low, so you'll probably see a lot of variants between 10-50 that will be emitted but not pass your Q50 threshold.

    what is your question about the -out_mode?

  • naarkhoonaarkhoo Posts: 38Member

    Thanks; I will try a higher value for emit_conf,

    Regarding -out_mode, I expect to see less variants by setting it to EMIT_ALL_CONFIDENT_SITES, than EMIT_VARIANTS_ONLY, but I got more ! can you explain a little bit the differences ?!

    Thank you

  • CarneiroCarneiro Posts: 274Administrator, Dev admin

    emit all confident sites will emit all sites including reference that pass your Q10 confidence threshold for emitting.

Sign In or Register to comment.