This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Unable to generate VCF using bcftools
Exception in thread "main" htsjdk.tribble.TribbleException$MalformedFeatureFile: Unable to parse header with error: Input stream does not contain a BCF encoded file; BCF magic header info not found, at record 0 with position 0:, for input source: file:///home/arshiya/efs/FastQC/ERR101899_sorted.bcf
    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 htsjdk.tribble.AbstractFeatureReader.getFeatureReader(AbstractFeatureReader.java:80)
    at htsjdk.variant.vcf.VCFFileReader.<init>(VCFFileReader.java:140)
    at htsjdk.variant.vcf.VCFFileReader.<init>(VCFFileReader.java:92)
    at picard.vcf.VcfFormatConverter.doWork(VcfFormatConverter.java:97)
    at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:305)
    at picard.cmdline.PicardCommandLine.instanceMain(PicardCommandLine.java:103)
    at picard.cmdline.PicardCommandLine.main(PicardCommandLine.java:113)
Caused by: htsjdk.tribble.TribbleException: Input stream does not contain a BCF encoded file; BCF magic header info not found, at record 0 with position 0:
    at htsjdk.variant.bcf2.BCF2Codec.error(BCF2Codec.java:482)
    at htsjdk.variant.bcf2.BCF2Codec.readHeader(BCF2Codec.java:153)
    at htsjdk.variant.bcf2.BCF2Codec.readHeader(BCF2Codec.java:58)
    at htsjdk.tribble.TribbleIndexedFeatureReader.readHeader(TribbleIndexedFeatureReader.java:261)
    ... 10 more
next-gen
• 2,175 views • link
updated by andrew.j.skelton73 andrew.j.skelton73 • written by arshiya.g02 • arshiya.g02

Thank you.. What is the command for converting the BAM file to BCF

• 0 views • link
written by arshiya.g02 • arshiya.g02

1 answer

your file 'ERR101899_sorted.bcf' is not a BCF file

check this with:

file ERR101899_sorted.bcf

output should be something like " Binary Call Format (BCF) version 2.1" or "compressed gzip file"

otherwise show us the output of

head ERR101899_sorted.bcf
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

This is what I get when I use the header command 

chr1    633613  A   1   ^!. 2
chr1    633614  A   1   .   3
chr1    633615  A   1   .   9
chr1    633616  A   1   .   =
chr1    633617  T   1   .   M
chr1    633618  A   1   .   [
chr1    633619  T   1   .   [
chr1    633620  T   1   .   g
chr1    633621  A   1   .   ]
chr1    633622  A   1   .   g
• 0 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by arshiya.g02 • arshiya.g02

it looks like the output of samtools mpileup. There is something wrong in your worklow.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this
  • 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 …
  • 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 …
  • No header in VCF file
    written by Begonia_pavonina • Begonia_pavonina
    Hello everyone, I am working with a specific variant calling pipeline, and the output is a VCF file missing headers. It seems there are no …
  • Visualizing VCF file on IGV
    written by Buffo Buffo
    Hi all, This is my first time trying to understand and visualize a VCF file on IGV, I make my vcf as follows: samtools mpileup …
  • 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 …
  • Problem with vcf format to run gatk FastaAlternateReferenceMaker
    written by pablo pablo
    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 …
  • Changing BCF genome coordinates (Picard LiftOverVcf error)
    written by wpierrick • wpierrick
    Hello all, I am trying to convert my BCF files from GrCH38 to GrCH37. I am using Picard LiftOverVcf as it can work with BCF …