A question regarding synteny

[Rahimlou]

Just for my understanding, is it possible to show synteny for certain genes (i.e. some part of the genomes)? What function can be used to zoom a specific region in several genomes?

[thackl]

Yes, that is possible. The function to do that is focus. You can do something like this:

# full genomes
p1 <- gggenomes(emale_genes, emale_seqs) %>%
 add_sublinks(emale_prot_ava) +
 geom_seq() + geom_gene() + geom_gene_tag(aes(label=name)) + geom_link()
p1

# zoom in on major capsid protein
p1 %>% focus(name=="MCP", .expand=2000)