Entering edit mode
Hi, Val.
Is there in interest in simply ignoring unknown INFO and GENOTYPE
fields
when parsing VCF files, perhaps by issuing a warning instead of an
error?
There are LOTS of malformed VCF files out there. In some cases, they
are
not useable, but in this case, they can be perfectly useable if these
unknown fields are simply ignored.
> dat = readVcf('tmp.gatk.vcf',genome='hg19')
Error: scanVcf: record 22 INFO 'KGPilot123' not found
path:
/Volumes/CCRBioinfo/projects/RosenbergImmuneStudy/staging/tmp.gatk.vcf
Thanks,
Sean
[[alternative HTML version deleted]]