[BioC] Feature request in readVcf
Martin Morgan
mtmorgan at fhcrc.org
Sun Sep 23 06:16:51 CEST 2012
On 09/22/2012 04:54 AM, Michael Lawrence wrote:
> Just in case you don't know, it is possible to work around this using the
> "info" slot of ScanVcfParam. But yes having a simple option to ignore
> unknown INFO and FORMAT fields would be convenient. The same should also
> apply to FILTER and ALT (in the case of structural variants).
unknown INFO and FORMAT fields now (v. 1.3.30) generate a warning and
don't get parsed. Seems like FILTER and ALT don't currently get checked;
are you looking for a warning on invalid value?
Martin
>
> Michael
>
> On Fri, Sep 21, 2012 at 10:51 AM, Sean Davis <sdavis2 at mail.nih.gov> wrote:
>
>> 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]]
>>
>> _______________________________________________
>> Bioconductor mailing list
>> Bioconductor at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>> Search the archives:
>> http://news.gmane.org/gmane.science.biology.informatics.conductor
>>
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
>
--
Computational Biology / Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N.
PO Box 19024 Seattle, WA 98109
Location: Arnold Building M1 B861
Phone: (206) 667-2793
More information about the Bioconductor
mailing list