This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Why is PyVCF's Reader failing to read VCF file?

Hello,

Does anyone know why I am getting this error with the Reader function in the PyVCF module. The VCF file used here was read normally using a Perl VCF module.

>>> fl = vcf.Reader(filename='sample.vcf.gz', compressed=True)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File ".local/lib/python2.7/site-packages/PyVCF-0.6.7-py2.7-linux-x86_64.egg/vcf/parser.py", line 288, in __init__
    self._parse_metainfo()
  File ".local/lib/python2.7/site-packages/PyVCF-0.6.7-py2.7-linux-x86_64.egg/vcf/parser.py", line 309, in _parse_metainfo
    key, val = parser.read_info(line)
  File ".local/lib/python2.7/site-packages/PyVCF-0.6.7-py2.7-linux-x86_64.egg/vcf/parser.py", line 123, in read_info
    "One of the INFO lines is malformed: %s" % info_string)
SyntaxError: One of the INFO lines is malformed: ##INFO=<ID=SOMATIC,Number=.,Description="Somatic status of existing variant (from ensembl VEP) associated with the most severe consequence in the CQ field">

Many thanks

software-error snp pyvcf
• 4,402 views • link
updated by Ram Ram • written by wrel2015 • wrel2015

1 answer

Hi there,

I think PyVCF requires that you add a "Type" descriptor to your INFO lines. The valid options are "Integer", "Float", "Flag", "Character", or "String".

So if you change the line above to look like this:

##INFO=<ID=SOMATIC,Number=.,Type=String,Description="Somatic status of existing variant (from ensembl VEP) associated with the most severe consequence in the CQ field">

It should work. Hope that helps!

• 0 views • link
updated by Ram Ram • written by brian@codenhance.com brian@codenhance.com

thank you Brian. This did quite down the INFO error but now there is a "contig lines malformation" error, I guess because the contig lines come after the INFO lines in the header and they too lack a "Type=String" attribute.

Is there an option in the Reader function of PyVCF that forces header lines to be read as strings or avoids me having to manipulate the header lines to make it work? For example if I were to read thousands of VCF files.

• 0 views • link
written by wrel2015 • wrel2015

Log in to answer this question.

More posts like this
  • MSIPred
    written by a.basitkhan1990 • a.basitkhan1990
    Hello! I am new to python (first time) I am trying to use this on my MAF data: https://github.com/wangc29/MSIpred I attempt to run the first …
  • pVAC-Seq (personalized Variant Antigens by Cancer Sequencing) with qsub
    written by soribaby88 • soribaby88
    Hi. I'm trying to do a job with pVAC-seq. But, a error occured and I don't understand the error message. My command /opt/sge/bin/lx-amd64/qsub /src/qsub.8 /data/D161740/Research/DongJin/01.programs/Python-3.5.6/python …
  • error while Kissplice run
    written by KG • KG
    I am using kissplice to identify SNPs from the transcriptome data. I am getting an error. This is how my terminal looks like: kissplice2refgenome -o …
  • kover learn issue
    written by alhicks32188 • alhicks32188
    Hi, I have run kover (http://aldro61.github.io/kover/) successfully previously, but in attempting to run in on a new dataset, I was able to create the dataset …
  • problems installing HTSeq
    written by Lila M Lila M
    Hi everybody I'm analyzing some data using DEXSEq, for that proposal I've used python2.7 home/my_user/R/x86_64-pc-linux-gnu-library/3.3/DEXSeq/python_scripts/dexseq_prepare_annotation.py Homo_sapiens.GRCh37.75.gtf Homo_sapiens.GRC37.gff3 And I get the next error mesage Could …
  • Error in running kissplice2reftranscriptome
    written by c.evang • c.evang
    Dear All and kissplice-users, I'm run into an error when running kissplice2reftranscriptome; the error was: Opening all files and formatting tables... Collecting data... Traceback (most …
  • Oncotator v1.8.0.0 usage, download database.
    written by Peavil • Peavil
    Hi. I used oncotator v1.8.0.0 with old-version database, It was running well. But, when I downloaded new database(oncotator_v1_ds_April052016.tar.gz. ) at (http://gatkforums.broadinstitute.org/gatk/discussion/4154/howto-install-and-run-oncotator-for-the-first-time), It doesn't run, and …
  • Error while running MACS2 ATAC-seq data
    written by Benn Benn
    Hi Guys, I am trying to run MACS2 for ATAC-seq data, it worked with a part of my bam files, but with another set of …
  • python's HGVS module UTA connection error
    written by wrel2015 • wrel2015
    Hi, I am getting an error whenever I try to connect to the public UTA database hosted by Invitae as part of the HGVS module …
  • Troubleshooting MACS bdgcmp
    written by mparker2 • mparker2
    I am trying to generate a log FE Bedgraph file from some chip seq data. I already have some bedgraph tracks for the control and …