Genome Analysis Toolkit

Variant Discovery in High-Throughput Sequencing Data

GATK process banner

Need Help?

Search our documentation

Community Forum

Hi, How can we help?

Developed in the Data Sciences Platform at the Broad Institute, the toolkit offers a wide variety of tools with a primary focus on variant discovery and genotyping. Its powerful processing engine and high-performance computing features make it capable of taking on projects of any size. Learn more

Spurios indel

0

1 comment

  • Avatar
    Gökalp Çelik

    Hi ulitskyi

    Looks like you are using amplicon data with very high depths. When I looked at the toy data you provided that seemingly HOMVAR SNP call is also called as HET by HaplotypeCaller which sure tells me that this kind of data is usually not very compatible with HaplotypeCaller's default assembly engine. 

    Here are a couple of parameters that I recommend you to test. You need to add these parameters to HaplotypeCaller. My initial test showed me that SNP is now called as HOMVAR and the spurious INDEL is not called at all. 

    --pileup-detection true 
    --use-pdhmm true
    --pileup-detection-enable-indel-pileup-calling true

    These parameters enable the partial hmms generated by the pileup detection mode to be added to the assembly engine to correct failed assemblies. Last parameter does not showup in the help text but believe us it is there. All of these parameters are also enabled by default when DRAGEN 3.7.8 compatibility mode is enabled. One thing to keep in mind is that when PDHMM is enabled GVCF output cannot be used (At least currently). 

    Normally when HaplotypeCaller is used without kmer and mbq parameters that indel site is not called at all and SNP site is called as HOMVAR correctly albeit AD values are not quite right. You might want to check your added parameters again if my additional parameters don't quite do justice for your calls. 

    I hope this helps.

    Regards. 

     

    0
    Comment actions Permalink

Please sign in to leave a comment.

Powered by Zendesk