# NOT RUN {
# GENOME.class <- readVCF("chr1.vcf.gz",1000,"1",1,100000)
# GENOME.class.split <- split_data_into_GFF_features(GENOME.class,"Homo_sapiens.GRCh37.73.gtf",
# "1", "gene")
# GENOME.class.split@region.names
# }
Run the code above in your browser using DataLab