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

GenomicsDBImport error from # fields in AS_UNIQ_ALT_READ_COUNT

Answered
0

3 comments

  • Avatar
    Genevieve Brandt (she/her)

    Hi Jeff Gaither,

    This issue is coming from a bug with the AS_UNIQ_ALT_READ_COUNT. The annotation is not meant for germline use, so you won't need to have it in your file. 

    Here is another issue with that annotation in github, it's currently still open: https://github.com/broadinstitute/gatk/issues/7298. I commented with your example.

    For the rest of your files, I would recommend not adding this annotation. You can just stick to the standard annotations with -G StandardAnnotation and -G AS_StandardAnnotation in HaplotypeCaller.

    Please let me know if you have any more questions.

    Best,

    Genevieve

    1
    Comment actions Permalink
  • Avatar
    Jeff Gaither

    Hi Genevieve, 

    Cool, then I'll just cut the AS_UNIQ_ALT_READ_COUNT from future runs. Thanks for your help!

    To anyone coming across this issue, removing the annotation from your VCF is as simple as 

    bcftools annotate -x INFO/AS_UNIQ_ALT_READ_COUNT your.vcf.gz

    Jeff

    0
    Comment actions Permalink
  • Avatar
    Genevieve Brandt (she/her)

    No problem! Glad we could figure out the issue.

    0
    Comment actions Permalink

Please sign in to leave a comment.

Powered by Zendesk