Entering edit mode
Hello! Could you help me please.
I work with vcf file:
> vcf <- readVcf("test.vcf.gz", "hg19")
Then I use predictCoding() to get information about nucleotide substitutions and output looks like this:
> coding <- predictCoding(vcf_mod, txdb, Hsapiens) > coding[5:7] GRanges with 3 ranges and 17 metadata columns: 1:865692_C/G chr1 [865692, 865692] + | <NA>
I'd like to save this file on local disk (I mean coding).How can I do it?
Thank you!
With devel rtracklayer, you can:
This is stored and indexed in a manner that is efficient for restricted range queries:
Thank you very much!)