Operon boxes about gggenomes HOT 2 CLOSED

library(tidyverse)
library(gggenomes)

genes <- tibble(
  seq_id = rep(c("s1", "s2"), each=5), # five genes per sequence
  start = c(50, 250, 500, 1000, 1300, 50, 450, 700, 1200, 1500),
  end = c(150, 450, 950, 1200, 1800, 350, 650, 1150, 1400, 2000) 
)
  
operons <-  tibble(
  seq_id = c("s1", "s2"),
  start = c(225, 425),
  end = c(1225,1425),
  name= c("operon X", "operon Y"),
)

gggenomes(genes, feats=operons) +
  geom_feat(size=5) +
  geom_feat_note(aes(label=name), nudge_y = -.1) +
  geom_gene()

# write example data to file for demonstation
# bed format: seq_id start end strand score name ...
write_tsv(operons, "operons.bed", col_names=FALSE)

# read a coordinate file (such as .bed, .gff, .gbk)
operons.bed <- read_feats("operons.bed")

gggenomes(genes, feats=operons.bed) +
  geom_feat(size=5) +
  geom_feat_note(aes(label=name), nudge_y = -.1) +
  geom_gene()