search

Bioconductor / VariantAnnotation

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

Unable to Export VCF After rbind #30

Open DarioS opened 5 years ago

DarioS commented 5 years ago

A minimal example using one file, although I noticed the problem when using two distinct files.

vcf <- readVcf(system.file("extdata", "ex2.vcf", package="VariantAnnotation") , "hg19",
               param = ScanVcfParam(geno = NA))
vcf <- rbind(vcf, vcf)
writeVcf(vcf, file = "test.vcf")

Error is "all 'geno()' fields must be named".