This is a test version of Biostars. For the public version, visit https://www.biostars.org.
vcfbreakmulti Broken VCF header, no column names?

I'm trying to use vcfbreakmulti on some vcf file:

vcfbreakmulti fileIN.vcf > fileOUT.vcf

I'm getting:

Broken VCF header, no column names?
 at /home/krzysiek/Programy/vcftools-vcftools-1d27c24/src/perl/Vcf.pm line 172
        Vcf::throw('Vcf4_2=HASH(0x2074ef8)', 'Broken VCF header, no column names?') called at /home/vcftools/src/perl/Vcf.pm line 867

My vcftools version is 0.1.15. fileIN.vcf is vcf version 4.1.

grep "#CHROM" fileIN.vcf
#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO    FORMAT  sampleName

Any ideas?

Edit: I've traced the error, it was not related to vcfbreakmulti. I was running it together with vcf-annotate in a pipe and vcf-annotate was trying to save to the same file as it was reading. Thanks everybody for suggestions.

vcftools
• 3,446 views • link
written by vpnfatrsqlbahzpchbvi • vpnfatrsqlbahzpchbvi

please, show us 

head -n1 fileIN.vcf

and

grep -B 5 "#CHROM" fileIN.vcf
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

head -n1 fileIN.vcf

##fileformat=VCFv4.1

grep -B 5 "#CHROM" fileIN.vcf

##INFO=<ID=OPOS,Number=.,Type=Integer,Description="List of original allele positions">
##INFO=<ID=OREF,Number=.,Type=String,Description="List of original reference bases">
##INFO=<ID=OALT,Number=.,Type=String,Description="List of original variant bases">
##INFO=<ID=OMAPALT,Number=.,Type=String,Description="Maps OID,OPOS,OREF,OALT entries to specific ALT alleles">
##deamination_metric=0.162637
#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO    FORMAT  sampleName
• 0 views • link
written by vpnfatrsqlbahzpchbvi • vpnfatrsqlbahzpchbvi

One suggestion would be to use vcf-validator from vcftools to check for any specific issues with the vcf.

• 0 views • link
written by microfuge microfuge 
 vcf-validator fileIN.vcf
    The header tag 'contig' not present for CHROM=chr1. (Not required but highly recommended.)
    The header tag 'contig' not present for CHROM=chr2. (Not required but highly recommended.)
    The header tag 'contig' not present for CHROM=chr3. (Not required but highly recommended.)
...
    The header tag 'contig' not present for CHROM=chr21. (Not required but highly recommended.)
    The header tag 'contig' not present for CHROM=chr22. (Not required but highly recommended.)
    The header tag 'contig' not present for CHROM=chrX. (Not required but highly recommended.)

I've got the same results for files that are working with vcfbreakmulti with no problem.

• 0 views • link
written by vpnfatrsqlbahzpchbvi • vpnfatrsqlbahzpchbvi

0 answers

No answers yet.

Log in to answer this question.

More posts like this