This is a test version of Biostars. For the public version, visit https://www.biostars.org.
VCF annotations -> readable info using PyVCF?

Hey, so I am working on a script that will run a VCF through several different databases for annotations, and then take the resulting VCF's and extract all of the new information each database added and put it into json files and tsv files. The issue I'm having is that once I load the ouput file into PyVCF, I am not sure how to extract information from the comments (which define what each of the annotations mean) and match that with annotations in the info fields. Is there a way to do this, or should I just be comparing every database manually in PyVCF instead of outsourcing the annotation?

pyvcf

1 answer

Have a look at this:

http://pyvcf.readthedocs.org/en/latest/API.html#vcf.Reader.infos

Oh, of course. Thanks so much!

Log in to answer this question.