search

zhengxwen / SeqArray

Data management of large-scale whole-genome sequence variant calls (Development version only)
http://www.bioconductor.org/packages/SeqArray
43 stars 12 forks source link

The created BGZF VCF file is truncated #47

Closed zhengxwen closed 5 years ago

zhengxwen commented 5 years ago 
library(SeqArray)

f <- seqOpen(seqExampleFileName("gds"))
seqGDS2VCF(f, "tmp.vcf.gz")

q("no")

Run the following shell script, and the output VCF file is not completed (SeqArray <= v1.23.7), because the file is not closed appropriately:

gunzip -c tmp.vcf.gz
zhengxwen commented 5 years ago

It is fixed in the commit (SeqArray_1.23.8).