[BioC] readVcf bgzip error

[Valerie Obenchain]

Brian,

Thanks for sending the sample data. I've tested the file on Windows and 
am not able to reproduce the error. It looks like you've got some old 
packages in your installation. A new version of Bioconductor (2.11) was 
just released last week. To update your packages you can follow the 
instructions here,

http://bioconductor.org/install/

Let me know if you still get the error after updating.

Valerie

Output for testing in devel :

## Read in the uncompressed file
 > library(VariantAnnotation)
 > fl <- "vcftest.vcf"
 > vcf <- readVcf(fl, "hg19")
 > vcf
class: VCF
dim: 700 997
genome: hg19
exptData(1): header
fixed(4): REF ALT QUAL FILTER
info(39): NS DP ... HD HP
geno(6): GT VR ... GQ FT
rownames(700): 1:6711146 1:6711190 ... 1:9662200 1:9662234
rowData values names(1): paramRangeID
colnames(997): A00003 A00057 ... A16420 A16457
colData names(1): Samples

## Compress and read in
 > compressVcf <- bgzip(fl, tempfile())
 > idx <- indexTabix(compressVcf, "vcf")
 > tab <- TabixFile(compressVcf, idx)
 > cmp <- readVcf(tab, "hg19")
 > cmp
class: VCF
dim: 700 997
genome: hg19
exptData(1): header
fixed(4): REF ALT QUAL FILTER
info(39): NS DP ... HD HP
geno(6): GT VR ... GQ FT
rownames(700): 1:6711146 1:6711190 ... 1:9662200 1:9662234
rowData values names(1): paramRangeID
colnames(997): A00003 A00057 ... A16420 A16457
colData names(1): Samples

 > sessionInfo()
R Under development (unstable) (2012-10-05 r60879)
Platform: i386-w64-mingw32/i386 (32-bit)

locale:
[1] LC_COLLATE=English_United States.1252
[2] LC_CTYPE=English_United States.1252
[3] LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C
[5] LC_TIME=English_United States.1252

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

other attached packages:
[1] VariantAnnotation_1.5.3 Rsamtools_1.11.1        Biostrings_2.27.2
[4] GenomicRanges_1.11.2    IRanges_1.17.2          BiocGenerics_0.5.

loaded via a namespace (and not attached):
  [1] AnnotationDbi_1.21.2   Biobase_2.19.0         biomaRt_2.15.0
  [4] bitops_1.0-4.1         BSgenome_1.27.1        DBI_0.2-5
  [7] GenomicFeatures_1.11.1 parallel_2.16.0        RCurl_1.95-1.1
[10] RSQLite_0.11.2         rtracklayer_1.19.0     stats4_2.16.0
[13] tools_2.16.0           XML_3.95-0.1           zlibbioc_1.5.0




On 10/10/2012 06:33 AM, Davis, Brian wrote:
> Valerie,
>
> This is human subject data so I'll have to work on getting permissions to share on my end.  In the mean time I'll try to reproduce with 1000 genomes data.
>
>
> Brain
>
> -----Original Message-----
> From: Valerie Obenchain [mailto:vobencha at fhcrc.org]
> Sent: Tuesday, October 09, 2012 6:19 PM
> To: Davis, Brian
> Cc: bioconductor at r-project.org
> Subject: Re: [BioC] readVcf bgzip error
>
> Hi Brian,
>
> I'm not sure what's going on here. Can you point me to where you got this file or is it small enough to send?
>
> Valerie
>
>
>
>
> On 10/09/2012 12:06 PM, Davis, Brian wrote:
>> I'm seeing an error when I read in a compress vcf, but not when I read in the uncompressed vcf.  Can anyone point me in the right direction to figure out what I'm doing wrong?  I've tried this on 3 different vcfs with the same error (different record fails).
>>
>>> # read in a complete file
>>> fl<- "first10K.vcf"
>>> vcf<- readVcf(fl, "hg19")
>>> vcf
>> class: VCF
>> dim: 9934 998
>> genome: hg19
>> exptData(1): header
>> fixed(4): REF ALT QUAL FILTER
>> info(39): NS DP ... HD HP
>> geno(6): GT VR ... GQ FT
>> rownames(9934): 1:69270 1:69360 ... 1:19597392 1:19597396 rowData
>> values names(1): paramRangeID
>> colnames(998): A00003 A00057 ... A16457 ''
>> colData names(1): Samples
>>> # now try again but compress it first
>>> fl<- "first10K.vcf"
>>> compressVcf<- bgzip(fl, tempfile())
>>> idx<- indexTabix(compressVcf, "vcf")
>>> tab<- TabixFile(compressVcf, idx)
>>> vcf<- readVcf(tab, "hg19")
>> Error: scanVcf: record 4370 INFO '0/0:.:130:131:.:.' not found
>>     path:
>> C:\Users\bdavis2\AppData\Local\Temp\RtmpwrXcST\file1dc84cff4177
>>> sessionInfo()
>> R version 2.15.1 (2012-06-22)
>> Platform: x86_64-pc-mingw32/x64 (64-bit)
>>
>> locale:
>> [1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United
>> States.1252 [3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C
>> [5] LC_TIME=English_United States.1252
>>
>> attached base packages:
>> [1] stats     graphics  grDevices utils     datasets  methods   base
>>
>> other attached packages:
>> [1] VariantAnnotation_1.2.10 Rsamtools_1.8.6          Biostrings_2.24.1        GenomicRanges_1.8.13
>> [5] IRanges_1.14.4           BiocGenerics_0.2.0
>>
>> loaded via a namespace (and not attached):
>> [1] AnnotationDbi_1.18.1  Biobase_2.16.0        biomaRt_2.12.0        bitops_1.0-4.1
>>    [5] BSgenome_1.24.0       DBI_0.2-5             GenomicFeatures_1.8.3 grid_2.15.1
>>    [9] lattice_0.20-10       Matrix_1.0-9          RCurl_1.95-1.1        RSQLite_0.11.2
>> [13] rtracklayer_1.16.3    snpStats_1.6.0        splines_2.15.1        stats4_2.15.1
>> [17] survival_2.36-14      tools_2.15.1          XML_3.95-0.1          zlibbioc_1.2.0
>>
>> Brian
>>
>>
>> 	[[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