Attention:
The front line support team will be unavailable to answer questions until May 27th 2019 as we are celebrating Memorial Day. We will be back soon after. Thank you for your patience and we apologize for any inconvenience!

call variants for a single sample using HC

I am using HC and the following command to call variants for a single sample using the genotyping_mode DISCOVERY. This is still running.
java -Xmx6g -jar $GATK/GenomeAnalysisTK.jar \
-T HaplotypeCaller \
-R ${reference} \
-I ${result}/reduced_reads.bam \
-L chr21 \
--genotyping_mode DISCOVERY \
-stand_emit_conf 10 \
-stand_call_conf 30 \
-o ${result}/raw_variants.vcf


Came across another command in Best Practices document.
java–jarGenomeAnalysisTK.jar
–T HaplotypeCaller\
–R human.fasta \
–I sample1.bam \
–o sample1.g.vcf \
–ERC GVCF\
-­‐-­‐variant_index_typeLINEAR\
-­‐-­‐variant_index_parameter 128000

Can you please let me know the difference and what command do I need to call variants for a single sample?
Output is a gvcf which is an intermediate step followed by VCF file.

Best Answers

Answers

Sign In or Register to comment.