长基因组区域的gggenes：有可能分裂它吗？

Question

我正在使用 gggenome 来描绘一个包含许多基因的长区域。我想用几条不同的线来显示该区域。有谁知道可以吗？

这是现在的样子：

ggplot(genes, aes(xmin = start, xmax = end, y = molecule, label = gene)) + #, fill = gene
  geom_gene_arrow()

Answer 1

没有您的数据，但根据

gggenes

包中的示例，尝试

ggplot(genes, aes(xmin = start, xmax = end, y = molecule, label = gene)) +
facet_wrap(~ molecule, scales = "free")

Answer 2

这是我的数据框基因：

molecule    gene    start   end strand  direction   unknown
Seq1    unknown1    1   1158    +   forward Other
Seq1    unknown3    1155    2039    +   forward Other
Seq1    unknown167  2030    2485    +   forward HTJ
Seq1    unknown168  2600    3592    +   forward Uncharacterized
Seq1    unknown169  3769    3599    -   reverse Uncharacterized
Seq1    unknown170  3985    3770    -   reverse Uncharacterized
Seq1    unknown171  4191    3982    -   reverse Uncharacterized
Seq1    unknown18015    4697    4188    -   reverse Uncharacterized
Seq1    unknown172  5230    4694    -   reverse Uncharacterized
...