Service Notice: Due to the blizzard currently hammering the US Northeast, the Broad is shut down and the GATK forum will be mostly unattended while we hunker down and sip hot cocoa with marshmallows. Assuming the power stays on and we're able to dig ourselves out of the snow when it's all over, normal service should resume Wednesday or Thursday.

Why do you have to provide -I and -L if GENOTYPE_GIVEN_ALLELES is used: UnifiedGenotyper

hintzenhintzen Posts: 4Member
edited February 2013 in Ask the GATK team

Previously I have been running a command like this:

  java -jar /path/GenomeAnalysisTK.jar \ 

  -T UnifiedGenotyper \

  -R /path/human_g1k_v37.fasta \

  -et NO_ET \

  -K /path/key \

  -out_mode EMIT_ALL_SITES \

  --input_file /path/bam \

  -L /path/intervals \

  -gt_mode GENOTYPE_GIVEN_ALLELES \

  --alleles /path/vcf \

  --dbsnp /path/dbsnp_135.b37.vcf \

  -o /path/my.vcf

But I was reading the documentation again and I read this statement: GENOTYPE_GIVEN_ALLELES only the alleles passed in from a VCF rod bound to the -alleles argument will be used for genotyping

Which lead me to believe that there wasn't a need to include the lines: --input_file /path/bam \ -L /path/intervals \

because it would be redundant. But when I try to run without those line I get back an error message: Walker requires reads but none were provided.

Can you give an explaination as to why both of those lines AND GENOTYPE_GIVEN_ALLELES would be needed?

Post edited by hintzen on

Best Answer

Answers

Sign In or Register to comment.