Quantcast
Channel: Recent Discussions — GATK-Forum
Viewing all articles
Browse latest Browse all 12345

UnifiedGenotyper different results depending on the output mode

$
0
0

Hello

I'm using GATK UnifiedGenotyper (v.3.6) for variant calling (extremely relaxed settings, no filtration). I found a false-negative variant (chr7:117188684T>G) that is clearly visible in IGV. But if I run UnifiedGenotyper with the same settings but in EMIT_ALL_SITES mode, the variant is detected with high quality. Why the variant is not detected in EMIT_VARIANTS_ONLY mode?

Here is my commands:

java -jar GenomeAnalysisTK.jar -R /data/reference/hg19.fasta -T UnifiedGenotyper -I /data/sample_data.bam -stand_call_conf 1 -stand_emit_conf 1 -L /data/targets.bed --output_mode EMIT_VARIANTS_ONLY -o /data/emit_variants_only.vcf

java -jar GenomeAnalysisTK.jar -R /data/reference/hg19.fasta -T UnifiedGenotyper -I /data/sample_data.bam -stand_call_conf 1 -stand_emit_conf 1 -L /data/targets.bed --output_mode EMIT_ALL_SITES -o /data/emit_all_sites.vcf

Thanks,
Tamara


Viewing all articles
Browse latest Browse all 12345

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>