Entering edit mode
I am seeing an error in the devel version of VariantAnnotation for VCF to VRanges conversion after this commit.
> library(VariantAnnotation) > fl <- system.file("extdata", "chr22.vcf.gz", package="VariantAnnotation") > vcf <- readVcf(fl, "hg19") > as(vcf, "VRanges") Error in if allis.na(x)) || all(x == "PASS")) return(FilterMatrix(matrix(nrow = length(x), : missing value where TRUE/FALSE needed > traceback() 3: parseFilterStrings(rd$FILTER) 2: asMethod(object) 1: as(vcf, "VRanges") > sessionInfo() R version 3.3.0 (2016-05-03) Platform: x86_64-apple-darwin13.4.0 (64-bit) Running under: OS X 10.11.5 (El Capitan) locale: [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8 attached base packages: [1] stats4 parallel stats graphics grDevices utils datasets [8] methods base other attached packages: [1] VariantAnnotation_1.19.4 Rsamtools_1.25.0 [3] Biostrings_2.41.4 XVector_0.13.2 [5] SummarizedExperiment_1.3.6 Biobase_2.33.0 [7] GenomicRanges_1.25.8 GenomeInfoDb_1.9.1 [9] IRanges_2.7.11 S4Vectors_0.11.7 [11] BiocGenerics_0.19.1 loaded via a namespace (and not attached): [1] XML_3.98-1.4 GenomicAlignments_1.9.4 bitops_1.0-6 [4] DBI_0.4-1 RSQLite_1.0.0 GenomicFeatures_1.25.14 [7] zlibbioc_1.19.0 BiocParallel_1.7.4 tools_3.3.0 [10] BSgenome_1.41.2 biomaRt_2.29.2 RCurl_1.95-4.8 [13] rtracklayer_1.33.8 AnnotationDbi_1.35.3