This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Problem with vcf format to run gatk FastaAlternateReferenceMaker

Hello,

I'm trying to run gatk FastaAlternateReferenceMaker to get the FASTA files from my vcf file.

I run : gatk FastaAlternateReferenceMaker -R my_reference.fa -O my_output.fasta -V my_file.vcf

I get that error message :

Using GATK jar /opt/apps/gcc-8.1.0/gatk-4.1.0.0/gatk-package-4.1.0.0-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 /opt/apps/gcc-8.1.0/gatk-4.1.0.0/gatk-package-4.1.0.0-local.jar FastaAlternateReferenceMaker -R my_reference.fa -O my_output.fasta -V my_file.vcf    14:15:23.712 INFO  NativeLibraryLoader - Loading libgkl_compression.so from jar:file:/opt/apps/gcc-8.1.0/gatk-4.1.0.0/gatk-package-4.1.0.0-local.jar!/com/intel/gkl/native/libgkl_compression.so
14:15:30.426 INFO  FastaAlternateReferenceMaker - ------------------------------------------------------------
14:15:30.427 INFO  FastaAlternateReferenceMaker - The Genome Analysis Toolkit (GATK) v4.1.0.0
14:15:30.427 INFO  FastaAlternateReferenceMaker - For support and documentation go to https://software.broadinstitute.org/gatk/
14:15:30.430 INFO  FastaAlternateReferenceMaker - Executing as *** on Linux v3.10.0-1127.18.2.el7.x86_64 amd64
14:15:30.430 INFO  FastaAlternateReferenceMaker - Java runtime: Java HotSpot(TM) 64-Bit Server VM v1.8.0_45-b14
14:15:30.430 INFO  FastaAlternateReferenceMaker - Start Date/Time: November 3, 2020 2:15:23 PM CET
14:15:30.430 INFO  FastaAlternateReferenceMaker - ------------------------------------------------------------
14:15:30.430 INFO  FastaAlternateReferenceMaker - ------------------------------------------------------------
14:15:30.431 INFO  FastaAlternateReferenceMaker - HTSJDK Version: 2.18.2
14:15:30.432 INFO  FastaAlternateReferenceMaker - Picard Version: 2.18.25
14:15:30.432 INFO  FastaAlternateReferenceMaker - HTSJDK Defaults.COMPRESSION_LEVEL : 2
14:15:30.432 INFO  FastaAlternateReferenceMaker - HTSJDK Defaults.USE_ASYNC_IO_READ_FOR_SAMTOOLS : false
14:15:30.432 INFO  FastaAlternateReferenceMaker - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_SAMTOOLS : true
14:15:30.432 INFO  FastaAlternateReferenceMaker - HTSJDK Defaults.USE_ASYNC_IO_WRITE_FOR_TRIBBLE : false
14:15:30.432 INFO  FastaAlternateReferenceMaker - Deflater: IntelDeflater
14:15:30.432 INFO  FastaAlternateReferenceMaker - Inflater: IntelInflater
14:15:30.433 INFO  FastaAlternateReferenceMaker - GCS max retries/reopens: 20
14:15:30.433 INFO  FastaAlternateReferenceMaker - Requester pays: disabled
14:15:30.433 INFO  FastaAlternateReferenceMaker - Initializing engine
14:15:30.853 INFO  FeatureManager - Using codec VCFCodec to read file file:///my_file.vcf
14:15:30.859 INFO  FastaAlternateReferenceMaker - Shutting down engine
[November 3, 2020 2:15:30 PM CET] org.broadinstitute.hellbender.tools.walkers.fasta.FastaAlternateReferenceMaker done. Elapsed time: 0.12 minutes.
Runtime.totalMemory()=623378432
org.broadinstitute.hellbender.exceptions.GATKException: Error initializing feature reader for path my_file.vcf
        at org.broadinstitute.hellbender.engine.FeatureDataSource.getTribbleFeatureReader(FeatureDataSource.java:353)
        at org.broadinstitute.hellbender.engine.FeatureDataSource.getFeatureReader(FeatureDataSource.java:305)
        at org.broadinstitute.hellbender.engine.FeatureDataSource.<init>(FeatureDataSource.java:256)
        at org.broadinstitute.hellbender.engine.FeatureManager.addToFeatureSources(FeatureManager.java:234)
        at org.broadinstitute.hellbender.engine.FeatureManager.initializeFeatureSources(FeatureManager.java:208)
        at org.broadinstitute.hellbender.engine.FeatureManager.<init>(FeatureManager.java:155)
        at org.broadinstitute.hellbender.engine.GATKTool.initializeFeatures(GATKTool.java:417)
        at org.broadinstitute.hellbender.engine.GATKTool.onStartup(GATKTool.java:638)
        at org.broadinstitute.hellbender.engine.ReferenceWalker.onStartup(ReferenceWalker.java:36)
        at org.broadinstitute.hellbender.cmdline.CommandLineProgram.runTool(CommandLineProgram.java:136)
        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)
Caused by: htsjdk.tribble.TribbleException$MalformedFeatureFile: Unable to parse header with error: Your input file has a malformed header: The FORMAT field was provided but there is no genotype/sample data, for input source: my_file.vcf
        at htsjdk.tribble.TribbleIndexedFeatureReader.readHeader(TribbleIndexedFeatureReader.java:263)
        at htsjdk.tribble.TribbleIndexedFeatureReader.<init>(TribbleIndexedFeatureReader.java:102)
        at htsjdk.tribble.TribbleIndexedFeatureReader.<init>(TribbleIndexedFeatureReader.java:127)
        at htsjdk.tribble.AbstractFeatureReader.getFeatureReader(AbstractFeatureReader.java:120)
        at org.broadinstitute.hellbender.engine.FeatureDataSource.getTribbleFeatureReader(FeatureDataSource.java:350)
        ... 14 more
Caused by: htsjdk.tribble.TribbleException$InvalidHeader: Your input file has a malformed header: The FORMAT field was provided but there is no genotype/sample data
        at htsjdk.variant.vcf.AbstractVCFCodec.parseHeaderFromLines(AbstractVCFCodec.java:185)
        at htsjdk.variant.vcf.VCFCodec.readActualHeader(VCFCodec.java:111)
        at htsjdk.tribble.AsciiFeatureCodec.readHeader(AsciiFeatureCodec.java:79)
        at htsjdk.tribble.AsciiFeatureCodec.readHeader(AsciiFeatureCodec.java:37)
        at htsjdk.tribble.TribbleIndexedFeatureReader.readHeader(TribbleIndexedFeatureReader.java:261)
        ... 18 more

I show the header of my vcf file, which looks good :

##fileformat=VCFv4.2
##INFO=<ID=DP,Number=1,Type=Integer,Description="Total Depth">
##INFO=<ID=AF,Number=A,Type=Float,Description="Allele Frequency">
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
##FORMAT=<ID=HQ,Number=2,Type=Integer,Description="Haplotype Quality">
#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO    FORMAT
Super-Scaffold_100001   672149  .       C       T       .       PASS    DP=210;AF=0.51  GT:HQ
Super-Scaffold_100001   862122  .       A       T       .       PASS    DP=305;AF=0.5   GT:HQ
Super-Scaffold_100001   931168  .       C       A       .       PASS    DP=127;AF=0.5   GT:HQ
Super-Scaffold_100001   967240  .       C       T       .       PASS    DP=127;AF=0.5   GT:HQ

Any idea why I get Your input file has a malformed header: The FORMAT field was provided but there is no genotype/sample data ?

Bests

vcf gatk phasing
• 2,798 views • link
written by pablo pablo

I think the main problem is here:

Error initializing feature reader for path my_file.vcf

Is that vcf file in the right location?

• 0 views • link
written by GenoMax GenoMax

Yes, the vcf file is in the right location. Even if I use the full path, it does not work.

• 0 views • link
written by pablo pablo

The FORMAT field was provided but there is no genotype/sample data, for input source: my_file.vcf

there is a FORMAT column but there is no associated genotype in your vcf.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • index fasta file for HaplotypeCaller
    written by Mojtaba • Mojtaba
    I am going to work with HaplotypeCaller to realign indels after some processes. HaplotypeCaller asks me to define an indexed fasta reference <ref.fasta.fai> and offers …
  • How does one create a custom VCF file?
    written by manfi • manfi
    I'm running an analysis where I use Mutect2 to look for reads that code for a particular pathogenic variant. The variant is given in a …
  • GATK GetPileupSummaries Java heap space
    written by Manuel Sokolov Ravasqueira • Manuel Sokolov Ravasqueira
    I am using GATK GetPileupSummaries in the following way: GENOME="/FILES/HUMAN_REFERENCES/hg19.fa" RECBAM="/FILES/${patient_id}/${patient_id}.recalibrated.bam" intervals_list="/FILES/HUMAN_REFERENCES/wgs_calling_regions.v1.interval_list" GERM="/FILES/HUMAN_REFERENCES/small_exac_common_3-hg19.vcf" PON="/FILES/HUMAN_REFERENCES/Mutect2-WGS-panel- b37-hg19.vcf"export GERM="/FILES/HUMAN_REFERENCES/af-only-gnomad-hg19.raw.sites.vcf" VCF="/FILES/${patient_id}/${patient_id}.recalibrated.vcf" OUTPUT="/FILES/${patient_id}/${patient_id}.getpileupsummaries.table" srun /mnt/beegfs/apptainer/images/gatk4.sif gatk GetPileupSummaries \ -I $RECBAM …
  • Error in VQSR first step
    written by Yoosef • Yoosef
    Hello. After inserting following codes for VQSR first step: java -jar gatk-package-4.3.0.0-local.jar VariantRecalibrator -O /home/yousef/Desktop/haplotype_hg38_SNP_Recal.vcf --resource:hapmap,known=false,training=true,truth=true,prior=15.0 '/media/yousef/EEFC0BDBFC0B9CC9/Sequencing/Next_generation_sequencing/Index_files/VQSR_data/resources_broad_hg38_v0_hapmap_3.3.hg38.vcf' --resource:omni,known=false,training=true,truth=false,prior=12.0 '/media/yousef/EEFC0BDBFC0B9CC9/Sequencing/Next_generation_sequencing/Index_files/VQSR_data/resources_broad_hg38_v0_1000G_omni2.5.hg38.vcf' --resource:1000G,known=false,training=true,truth=false,prior=10.0 '/media/yousef/EEFC0BDBFC0B9CC9/Sequencing/Next_generation_sequencing/Index_files/VQSR_data/resources_broad_hg38_v0_1000G_phase1.snps.high_confidence.hg38.vcf' --resource:dbsnp,known=true,training=false,truth=false,prior=2.0 '/media/yousef/EEFC0BDBFC0B9CC9/Sequencing/Next_generation_sequencing/Index_files/VQSR_data/resources_broad_hg38_v0_Homo_sapiens_assembly38.dbsnp138.vcf' --tranches-file /home/yousef/Desktop/Tranches.txt …
  • While using GATK BaseRecalibration :No report file is created and i cannot figure out why is that…
    written by Phoebe Magdy • Phoebe Magdy
    I'm trying to make base recalibration using GATK for bam files and here is the command I used: ``` !gatk --java-options "-Xmx2G" BaseRecalibrator \ -R …
  • gatk: Invalid lookback distance found
    written by rk.khayami94 • rk.khayami94
    I tried to run the ```FilterVariantTranches``` tool from GATK 4.2.4.1 docker image and encountered the following error: ``` java.lang.RuntimeException: Invalid lookback distance found. at com.intel.gkl.compression.IntelInflater.inflateNative(Native …
  • GATK Mutect2 errors during basic variant calling
    written by michael.flower.14 • michael.flower.14
    I've just installed GATK and am trying to do some basic variant calling. However when I try and run this line gatk Mutect2 -R $REF …
  • GATK HaplotypeCaller Error
    written by sahar.salimi92 • sahar.salimi92
    Hi, I run all the steps below to run GATK gatk HaplotypeCaller -R chrX.fa -I ERR188044_chrX.bam -O ERR188044_chrX.vcf -ERC GVCF But I got the following …
  • GATK BaseRecalibrator error.. how do I solve this?
    written by kwanghoon • kwanghoon
    I downloaded reference genome and known site vcf at GATK bundle (https://console.cloud.google.com/storage/browser/genomics-public-data/resources/broad/hg38/v0;tab=objects?pli=1&prefix=&forceOnObjectsSortingFiltering=false) I sorted and duduplicated bam file then tried to use baserecalibrator. But I …
  • gatk haplotypecaller crashing
    written by dishasharma35 • dishasharma35
    I am trying to run GATK HAPLOTYPECALLER and dit is crashing. Could someone help? numactl -N 0 -m 0 /home/administrator/IGIB/gatk-4.1.4.0/gatk --java-options "-Xmx50g" HaplotypeCaller -R "/home/administrator/IGIB/samples/Homo_sapiens_assembly38.fasta" …