Hello
I am studying genotyping by RAD-seq.
I have interests in population genetic structure within wild plant species.
So I would like to use GATK for genotyping / variant calling.
However,
I felt in a problem.
My working sentence is
"java -jar gatk-package-4.1.2.0-local.jar HaplotypeCaller -R /mnt/c/Popgen -I /mnt/c/Popgen/samtools-1.9/sample18.bam -O /mnt/c/Popgen/GATK/sample18.vcf"
and then,
I got this result.
-------------------------------------------------------------------------------------------------------------------------------------------
11:45:13.056 INFO NativeLibraryLoader - Loading libgkl_compression.so from jar:file:/mnt/c/popgen/GATK/gatk-package-4.1.2.0-local.jar!/com/intel/gkl/native/libgkl_compression.so
May 07, 2020 11:45:14 AM shaded.cloud_nio.com.google.auth.oauth2.ComputeEngineCredentials runningOnComputeEngine
INFO: Failed to detect whether we are running on Google Compute Engine.
11:45:14.793 INFO HaplotypeCaller - ------------------------------------------------------------
11:45:14.793 INFO HaplotypeCaller - The Genome Analysis Toolkit (GATK) v4.1.2.0
11:45:14.794 INFO HaplotypeCaller - For support and documentation go to https://software.broadinstitute.org/gatk/
11:45:14.795 INFO HaplotypeCaller - Executing as bwook1243@bwook1243 on Linux v4.4.0-18362-Microsoft amd64
11:45:14.795 INFO HaplotypeCaller - Java runtime: OpenJDK 64-Bit Server VM v1.8.0_242-8u242-b08-0ubuntu3~18.04-b08
11:45:14.796 INFO HaplotypeCaller - Start Date/Time: May 7, 2020 11:45:13 AM KST
11:45:14.797 INFO HaplotypeCaller - ------------------------------------------------------------
11:45:14.798 INFO HaplotypeCaller - ------------------------------------------------------------
11:45:14.799 INFO HaplotypeCaller - HTSJDK Version: 2.19.0
11:45:14.799 INFO HaplotypeCaller - Picard Version: 2.19.0
11:45:14.800 INFO HaplotypeCaller - HTSJDK Defaults.COMPRESSION_LEVEL : 2
11:45:14.801 INFO HaplotypeCaller - HTSJDK Defaults.USE_ASYNC_IO_READ_FOR_SAMTOOLS : false
11:45:14.801 INFO HaplotypeCaller - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_SAMTOOLS : true
11:45:14.802 INFO HaplotypeCaller - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_TRIBBLE : false
11:45:14.804 INFO HaplotypeCaller - Deflater: IntelDeflater
11:45:14.808 INFO HaplotypeCaller - Inflater: IntelInflater
11:45:14.809 INFO HaplotypeCaller - GCS max retries/reopens: 20
11:45:14.809 INFO HaplotypeCaller - Requester pays: disabled
11:45:14.810 INFO HaplotypeCaller - Initializing engine
11:45:14.815 INFO HaplotypeCaller - Shutting down engine
[May 7, 2020 11:45:14 AM KST] org.broadinstitute.hellbender.tools.walkers.haplotypecaller.HaplotypeCaller done. Elapsed time: 0.03 minutes.
Runtime.totalMemory()=288882688
java.lang.IllegalArgumentException: File is not a supported reference file type: /mnt/c/Popgen
at htsjdk.samtools.reference.ReferenceSequenceFileFactory.lambda$getFastaExtension$0(ReferenceSequenceFileFactory.java:252)
at java.util.Optional.orElseGet(Optional.java:267)
at htsjdk.samtools.reference.ReferenceSequenceFileFactory.getFastaExtension(ReferenceSequenceFileFactory.java:252)
at htsjdk.samtools.reference.ReferenceSequenceFileFactory.getDefaultDictionaryForReferenceSequence(ReferenceSequenceFileFactory.java:222)
at org.broadinstitute.hellbender.utils.fasta.CachingIndexedFastaSequenceFile.checkFastaPath(CachingIndexedFastaSequenceFile.java:166)
at org.broadinstitute.hellbender.utils.fasta.CachingIndexedFastaSequenceFile.(CachingIndexedFastaSequenceFile.java:129)
at org.broadinstitute.hellbender.utils.fasta.CachingIndexedFastaSequenceFile.(CachingIndexedFastaSequenceFile.java:111)
at org.broadinstitute.hellbender.utils.fasta.CachingIndexedFastaSequenceFile.(CachingIndexedFastaSequenceFile.java:96)
at org.broadinstitute.hellbender.engine.ReferenceFileSource.(ReferenceFileSource.java:35)
at org.broadinstitute.hellbender.engine.ReferenceDataSource.of(ReferenceDataSource.java:27)
at org.broadinstitute.hellbender.engine.GATKTool.initializeReference(GATKTool.java:423)
at org.broadinstitute.hellbender.engine.GATKTool.onStartup(GATKTool.java:693)
at org.broadinstitute.hellbender.engine.AssemblyRegionWalker.onStartup(AssemblyRegionWalker.java:161)
at org.broadinstitute.hellbender.cmdline.CommandLineProgram.runTool(CommandLineProgram.java:137)
at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMainPostParseArgs(CommandLineProgram.java:191)
at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:210)
at org.broadinstitute.hellbender.Main.runCommandLineProgram(Main.java:162)
at org.broadinstitute.hellbender.Main.mainEntry(Main.java:205)
at org.broadinstitute.hellbender.Main.main(Main.java:291)
-------------------------------------------------------------------------------------------------------------------------------------------
I think that the connection to Google compute engine is not working on my computer. Isn't it?
I prepared reference files (fasta (*.fa), *fa.fai, *.dict).
How can I solve this problem?
Please give your kind help.