Error: ##fileformat=VCFv4.2 does not exist
4
0
Entering edit mode
16 months ago

Hello everybody,

I am using Pharmcat to preprocess my vcf file, and for this I am running this command

python3 pharmcat_vcf_preprocessor.py -vcf NA12801.VCF

But I am getting this error

Error: ##fileformat=VCFv4.2 does not exist

I have generated my vcf file by using gatk Haplotypecaller

./gatk --java-options "-Xmx4g" HaplotypeCaller -R GCF_000001405.40_GRCh38.p14_genomic.fna -I sample.bam -O sample.VCF

Kindly help me out in resolving this issue.

Thanks in advance

Pharmcat VCF • 1.1k views
ADD COMMENT
0
Entering edit mode
16 months ago

Have a look in your VCF file if that line exists.

eg

grep fileformat *.VCF

It might be VCF 4.1 or 4.3 etc.

ADD COMMENT
0
Entering edit mode

Yes this line exist in my vcf file

##fileformat=VCFv4.2
ADD REPLY
0
Entering edit mode
16 months ago
Lian ▴ 40

Look at the pharmcat_vcf_preprocessor.py program, post the code near the error Error: ##fileformat=VCFv4.2 does not existis.

ADD COMMENT
0
Entering edit mode
16 months ago
barslmn ★ 2.2k

Error: ##fileformat=VCFv4.2 does not exist

Looks like a file does not exist error and checking source it looks like its trying to read a txt file with paths to VCF files. So instead of giving the VCF directly give a txt file with a path to VCF file.

ADD COMMENT
0
Entering edit mode
16 months ago
sbstevenlee ▴ 480

I would also recommend checking out a different star allele-calling tool named PyPGx (see docs). It supports VCF files produced from HaplotypeCaller. Disclaimer: I'm the developer of PyPGx.

ADD COMMENT

Login before adding your answer.

Traffic: 1383 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6