如何在R中组合VCF文件？

Question

我在使用rbind使用bioconductor的libraryAnnotation组合VCF文件时遇到问题。

我正在读取两个VCF文件，当我尝试将它们与'rbind'按一定顺序组合时，却出现了错误，但按不同的顺序可以工作。

> vcf.full <- rbind(vcf1, vcf2) ## Doesnt work 
Error in as(from, to_class, strict = FALSE) : 
  no method or default for coercing “CompressedCharacterList” to “DNAStringSetList”
> vcf.full <- rbind(vcf2, vcf1) ## Work 
> 
> header(vcf1)
class: VCFHeader 
samples(503): NA06984 NA06989 ... NA20828 NA20832
meta(2): META contig
fixed(2): FILTER ALT
info(27): CIEND CIPOS ... EX_TARGET MULTI_ALLELIC
geno(1): GT
> header(vcf2)
class: VCFHeader 
samples(503): NA06984 NA06989 ... NA20828 NA20832
meta(2): META contig
fixed(2): FILTER ALT
info(27): CIEND CIPOS ... EX_TARGET MULTI_ALLELIC
geno(1): GT
> 
> dim(vcf1)
[1]  30 503
> dim(vcf2)
[1] 149 503

两个VCF文件似乎相同。知道为什么会发生吗？或对使用R组合多个VCF文件的其他方法有何建议？