Entering edit mode
                    I used annotatr to annotate my methylation data and got four transcripts that I can't find any information for, no ensembl gene id translation or gene symbol translation: ENST00000602953.1 ENST00000552238.1 ENST00000451960.1 Any idea where these are coming from?
## create Granges object
regions <- GRanges(
  seqnames = sig_windows_filtered$chrom,
  ranges = IRanges(sig_windows_filtered$start, sig_windows_filtered$end),
  mcols = list(
    diff = sig_windows_filtered$diff, 
    type = sig_windows_filtered$type
  )
)
# Select annotations for intersection with regions
annots = c('hg38_cpgs', 'hg38_basicgenes', 'hg38_enhancers_fantom')
# Build the annotations (a single GRanges object)
annotations = build_annotations(genome = 'hg38', annotations = annots)
# Intersect the regions read in with the annotations
dm_annotated = annotate_regions(
  regions = regions,
  annotations = annotations,
  ignore.strand = TRUE,
  quiet = FALSE)
# Coerce to a data.frame
df_dm_annotated = data.frame(dm_annotated)
