search

Bioconductor / VariantAnnotation

Annotation of Genetic Variants
https://bioconductor.org/packages/VariantAnnotation
23 stars 20 forks source link

scanBcfHeader _hts_rewind() error in 1.29.19 #22

Closed smgogarten closed 5 years ago

smgogarten commented 5 years ago

One of my unit tests in SeqArray is now failing to read an example file:

vcffile <- SeqArray::seqExampleFileName("vcf")
vcf <- readVcf(vcffile, genome="hg19")
Error in scanBcfHeader(bf) : [internal] _hts_rewind() failed
zhengxwen commented 5 years ago

The VCF file in SeqArray has been recompressed by bgzip, and it works with the latest VariantAnnotation.