The current GATK version is 3.7-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.
GATK 3.7 is here! Be sure to read the Version Highlights and optionally the full Release Notes.
Register now for the upcoming GATK Best Practices workshop, Feb 20-22 in Leuven, Belgium. Open to all comers! More info and signup at http://bit.ly/2i4mGxz

how to increase our call confidence

naarkhoonaarkhoo Member Posts: 38

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 -Djava.io.tmpdir=/tmp/evm-gatk/ \
-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 \
-out_mode EMIT_VARIANTS_ONLY \
--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 \

Thanks

Answers

Sign In or Register to comment.