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

Funcotator datasources.v1.7 gencode raise error

0

3 comments

  • Official comment
    Avatar
    Jonn Smith

    Tjitske - I responded to your github post, but in case you didn't see it:

     

    @Tjitskedv Typically we ask people make a new issue in github or create a forum post, rather than commenting on open pull requests.

    To answer your question directly - do not use v1.7 datasources with GATK 4.1.7.X or 4.1.8.X. Those data sources are not yet compatible with the Funcotator code in the Master branch. When we release GATK 4.1.9.0 they will work with Funcotator. This is not a bug in the code - it is an artifact of how we have to do data source releases.

    See issue #6708 for all the gory details.

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

    Hi Tjitske de Vries, it looks like there is an issue with compatibility between your files. Where did you get your GTF file? Also, are any of your data sources using hg38 instead of hg19? This will lead to errors.

    0
    Comment actions Permalink
  • Avatar
    Tjitske de Vries

    Hi Genevieve Brandt, thank you for your response. I downloaded the GTF file as part of the bundle  funcotator_dataSources.v1.7.20200521s.tar.gz from https://gatk.broadinstitute.org/hc/en-us/articles/360042912011-Funcotator

    al my datasources are hg19, I did not make any new data sources, only use those of the provided bundle.

    I have used a Mutect2 vcf file which was filtered by FilterMutectCalls, I use the exact same reference file for all steps. 

    Is there something that I can test/check/change for the compatibility? (removing gencode datasource is not an option since this is the only one that is required.) - do I need another version? 

    Hope you can help

    Update: I have downloaded the v1.6 data source bundle which contained the genecode V19 file now I used that gencode file instead of the one from v1.7, (copied it to the v1.7 folder so did not change any other file) and now it works. So I think the genecode-backmap file is not really compatible with funcotator for hg19, was this tested? or is this problem specific for my files? it seems to be a more general problem, so maybe good to check and update the files of the data source bundle?

    kind regards Tjitske

    0
    Comment actions Permalink

Please sign in to leave a comment.

Powered by Zendesk