[BioC] VariantAnnotation: Specifying 'seqinfo' at import with 'readVcf'

Julian Gehring julian.gehring at embl.de
Tue Sep 24 11:29:18 CEST 2013


Hi,

Is there a direct way to specifiy the 'seqinfo' of a genome for the 
import of a VCF file using 'readVcf'?  I'm aware that one can change it 
with the 'seqinfo' method afterwards, but for large VCF files this can 
take a significant amount of time.

An alternative would be to sneak it in by the 'which' arguments, such as:

readVcf(file, genome, ScanVcfParam(which = as(seq_info, "GRanges")))

but this requires the file to be indexed beforehand.

Best wishes
Julian



More information about the Bioconductor mailing list