Entering edit mode
I have a VCF created by Bis-SNP. One of the INFO fields is called CS but readVcf is ignoring it and I can't figure out why. A tiny example VCF can be downloaded from https://www.dropbox.com/sh/bmbyjgts26req8k/AAA4UkTTUU8IzNxkcy4ZP14Ua?dl=0 that can be used to reproduce the problem.
This is what I've tried:
> library(VariantAnnotation) > x <- readVcf('~/tmp/ex.vcf.gz', 'hg19') # Missing CS field. Should appear before Context. > info(x) DataFrame with 2 rows and 9 columns Context DB DP HQ MQ0 NS QD REF SB <CharacterList> <logical> <integer> <numeric> <integer> <integer> <numeric> <CharacterList> <numeric> rs55998931 YH TRUE 8 NA 0 1 NA CH -0.1412 chr1:10774_G/C SG FALSE 4 NA 0 1 NA CH -0.6254 > sessionInfo() R version 3.1.1 (2014-07-10) Platform: x86_64-apple-darwin10.8.0 (64-bit) locale: [1] en_AU.UTF-8/en_AU.UTF-8/en_AU.UTF-8/C/en_AU.UTF-8/en_AU.UTF-8 attached base packages: [1] stats4 parallel stats graphics grDevices utils datasets methods base other attached packages: [1] VariantAnnotation_1.12.2 Rsamtools_1.18.1 Biostrings_2.34.0 XVector_0.6.0 GenomicRanges_1.18.1 [6] GenomeInfoDb_1.2.2 IRanges_2.0.0 S4Vectors_0.4.0 BiocGenerics_0.12.0 loaded via a namespace (and not attached): [1] AnnotationDbi_1.28.1 base64enc_0.1-2 BatchJobs_1.4 BBmisc_1.7 Biobase_2.26.0 [6] BiocParallel_1.0.0 biomaRt_2.22.0 bitops_1.0-6 brew_1.0-6 BSgenome_1.34.0 [11] checkmate_1.5.0 codetools_0.2-9 DBI_0.3.1 digest_0.6.4 fail_1.2 [16] foreach_1.4.2 GenomicAlignments_1.2.0 GenomicFeatures_1.18.2 iterators_1.0.7 RCurl_1.95-4.3 [21] RSQLite_1.0.0 rtracklayer_1.26.1 sendmailR_1.2-1 stringr_0.6.2 tools_3.1.1 [26] XML_3.98-1.1 zlibbioc_1.12.0
Thanks for your help,
Pete
Thanks Pete. I'll have a look.
Valerie