GATK for Rnor6 - Funcotator error
AnsweredHello, I'm trying to find variant from rat sequences. I used Rnor6 ref but I got this error A USER ERROR has occurred: File file://rat_mutation/funcotator_dataSources.v1.7.20200521g/gencode/rn6/Rattus_norvegicus.Rnor_6.0.104.sort.gff3.gz is of the wrong type. It should contain Features of type GencodeGtfFeature, but instead contains Features of type(s): [Gff3Feature]
REQUIRED for all errors and issues:
a) GATK version used: v4.2.6.1
b) Exact command used:
c) Entire program log:
Using GATK jar /miniconda3/envs/ratmu/share/gatk4-4.2.6.1-1/gatk-package-4.2.6.1-local.jar
Running:
java -Dsamjdk.use_async_io_read_samtools=false -Dsamjdk.use_async_io_write_samtools=true -Dsamjdk.use_async_io_write_tribble=false -Dsamjdk.compression_level=2 -jar miniconda3/envs/ratmu/share/gatk4-4.2.6.1-1/gatk-package-4.2.6.1-local.jar Funcotator -V funcotator_dataSources.v1.7.20200521g/lmm_known/rn6/Rattus_norvegicus.vcf.gz --data-sources-path funcotator_dataSources.v1.7.20200521g -R funcotator_dataSources.v1.7.20200521g/gencode/rn6/genome.fa --ref-version rn6 --output-file-format VCF -O rat_mutation/variant/Rnor6_variants.funcotated.vcf --disable-sequence-dictionary-validation
11:05:56.589 INFO NativeLibraryLoader - Loading libgkl_compression.so from jar:file:miniconda3/envs/ratmu/share/gatk4-4.2.6.1-1/gatk-package-4.2.6.1-local.jar!/com/intel/gkl/native/libgkl_compression.so
11:05:56.724 INFO Funcotator - ------------------------------------------------------------
11:05:56.724 INFO Funcotator - The Genome Analysis Toolkit (GATK) v4.2.6.1
11:05:56.724 INFO Funcotator - For support and documentation go to https://software.broadinstitute.org/gatk/
11:05:56.724 INFO Funcotator - Executing as nkdang@cdragon054 on Linux v3.10.0-1160.15.2.el7.x86_64 amd64
11:05:56.724 INFO Funcotator - Java runtime: OpenJDK 64-Bit Server VM v1.8.0_332-b09
11:05:56.725 INFO Funcotator - Start Date/Time: August 3, 2022 11:05:56 AM CDT
11:05:56.725 INFO Funcotator - ------------------------------------------------------------
11:05:56.725 INFO Funcotator - ------------------------------------------------------------
11:05:56.725 INFO Funcotator - HTSJDK Version: 2.24.1
11:05:56.725 INFO Funcotator - Picard Version: 2.27.1
11:05:56.725 INFO Funcotator - Built for Spark Version: 2.4.5
11:05:56.725 INFO Funcotator - HTSJDK Defaults.COMPRESSION_LEVEL : 2
11:05:56.725 INFO Funcotator - HTSJDK Defaults.USE_ASYNC_IO_READ_FOR_SAMTOOLS : false
11:05:56.725 INFO Funcotator - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_SAMTOOLS : true
11:05:56.725 INFO Funcotator - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_TRIBBLE : false
11:05:56.725 INFO Funcotator - Deflater: IntelDeflater
11:05:56.725 INFO Funcotator - Inflater: IntelInflater
11:05:56.725 INFO Funcotator - GCS max retries/reopens: 20
11:05:56.725 INFO Funcotator - Requester pays: disabled
11:05:56.726 INFO Funcotator - Initializing engine
11:05:57.029 INFO FeatureManager - Using codec VCFCodec to read file file://rat_mutation/funcotator_dataSources.v1.7.20200521g/lmm_known/rn6/Rattus_norvegicus.vcf.gz
11:05:57.175 INFO Funcotator - Done initializing engine
11:05:57.175 INFO Funcotator - Skipping sequence dictionary validation.
11:05:57.175 INFO Funcotator - Processing user transcripts/defaults/overrides...
11:05:57.175 INFO Funcotator - Initializing data sources...
11:05:57.177 INFO DataSourceUtils - Initializing data sources from directory: funcotator_dataSources.v1.7.20200521g
11:05:57.177 INFO DataSourceUtils - Data sources version: 1.7.2020521g
11:05:57.178 INFO DataSourceUtils - Data sources source: ftp://gsapubftp-anonymous@ftp.broadinstitute.org/bundle/funcotator/funcotator_dataSources.v1.7.20200521g.tar.gz
11:05:57.178 INFO DataSourceUtils - Data sources alternate source: gs://broad-public-datasets/funcotator/funcotator_dataSources.v1.7.20200521.tar.gz
11:05:57.183 INFO DataSourceUtils - Resolved data source file path: file:///rat_mutation/Rattus_norvegicus.Rnor_6.0.104.sort.gff3.gz -> file:///rat_mutation/funcotator_dataSources.v1.7.20200521g/gencode/rn6/Rattus_norvegicus.Rnor_6.0.104.sort.gff3.gz
11:05:57.183 INFO DataSourceUtils - Resolved data source file path: file:///rat_mutation/genome.fa -> file://rat_mutation/funcotator_dataSources.v1.7.20200521g/gencode/rn6/genome.fa
11:05:57.186 INFO DataSourceUtils - Resolved data source file path: file://rat_mutation/Rattus_norvegicus.vcf.gz -> file://rat_mutation/funcotator_dataSources.v1.7.20200521g/lmm_known/rn6/Rattus_norvegicus.vcf.gz
11:05:57.186 INFO Funcotator - Finalizing data sources (this step can be long if data sources are cloud-based)...
11:05:57.186 INFO DataSourceUtils - Resolved data source file path: file://rat_mutation/Rattus_norvegicus.Rnor_6.0.104.sort.gff3.gz -> file:///rat_mutation/funcotator_dataSources.v1.7.20200521g/gencode/rn6/Rattus_norvegicus.Rnor_6.0.104.sort.gff3.gz
11:05:57.186 INFO DataSourceUtils - Setting lookahead cache for data source: Gencode : 100000
11:05:57.230 INFO Funcotator - Shutting down engine
[August 3, 2022 11:05:57 AM CDT] org.broadinstitute.hellbender.tools.funcotator.Funcotator done. Elapsed time: 0.01 minutes.
Runtime.totalMemory()=2120744960
***********************************************************************
A USER ERROR has occurred: File file://rat_mutation/funcotator_dataSources.v1.7.20200521g/gencode/rn6/Rattus_norvegicus.Rnor_6.0.104.sort.gff3.gz is of the wrong type. It should contain Features of type GencodeGtfFeature, but instead contains Features of type(s): [Gff3Feature]
***********************************************************************
Set the system property GATK_STACKTRACE_ON_USER_EXCEPTION (--java-options '-DGATK_STACKTRACE_ON_USER_EXCEPTION=true') to print the stack trace.
I appreciate all your help in advance.
Thanks
-
Hi Ngoc Khoi Dang,
Thank you for writing to the GATK forum! I hope that we can help you sort this out.
Unfortunately, Funcotator does not support non-human genomes at this point. That said, we are actively working on expanding access to other species. This issue is definitely on our radar; we hope to add this additional support in a few more months.
Please stay tuned! If you have any further questions, please feel free to reach out at any time.
Best,
Anthony
Please sign in to leave a comment.
1 comment