[BioC] readVcf bgzip error

Davis, Brian Brian.Davis at uth.tmc.edu
Fri Oct 12 22:33:21 CEST 2012


Valerie,

Thanks for looking into this.  I'm new to Biocondoctor, but not R.  I hadn't realized update.package didn't update Bioconductor.  Anyways, this did indeed fix my problem.

Again, thanks.

Brian

-----Original Message-----
From: Valerie Obenchain [mailto:vobencha at fhcrc.org] 
Sent: Friday, October 12, 2012 3:07 PM
To: Davis, Brian
Cc: bioconductor at r-project.org
Subject: Re: [BioC] readVcf bgzip error

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



More information about the Bioconductor mailing list