gggenes of a long genome region: is it possible to split it?

70 Views Asked by At

I am using gggenome to depict a long region with many genes. I would like to show the region split in several different lines. Does anyone know if it is possible?

This is the way now appear:

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

enter image description here

2

There are 2 best solutions below

0
markhogue On

Don't have your data, but based on the example in the gggenes pacakge, try

ggplot(genes, aes(xmin = start, xmax = end, y = molecule, label = gene)) +
facet_wrap(~ molecule, scales = "free")
1
Antonio J. Pérez On

This is my dataframe genes:

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
...