Hi, not reporting a bug, just a question. I have two genomes, but one is a large g

Different genomes scales about gggenomes HOT 4 CLOSED

tgodoy commented on September 26, 2024

Different genomes scales

from gggenomes.

Comments (4)

thackl commented on September 26, 2024

You mean zoom in on the second genome so it becomes larger on the plot? Hm, I can't think of anything that would be very easy.

Of the top of my head I can think of two maybe options:

What would be possible is to break up the long contig, i.e. only show the regions with some space around that actually align. focus() on links can do that automatically (let me know if you want try something like that but don't know how).
Manually scale the coordinates of one of the genomes. Something like multiplying every length, start and end column in for Xtrop_Chr1 in gAll, sAll and tblastx_tibble by 10 or so.

from gggenomes.

tgodoy commented on September 26, 2024

Thanks, I think that I will probe the the second solution. But In case I need to do a zoom in some regions, How I can do it?

from gggenomes.

thackl commented on September 26, 2024

Here's a minimal example that you might be able to adapt for your case.

library(gggenomes)

g0 <- tribble(
  ~seq_id, ~start, ~end,
  "A", 100, 200,
  "A", 10000, 10100,
  "B", 100, 200,
  "B", 300, 500
)

# clusters are just a faster way for manually creating links
c0 <- tribble(
  ~cluster_id, ~feat_id,
  "cls1", "f1",
  "cls1", "f3",
  "cls2", "f2",
  "cls2", "f4")

p1 <- gggenomes(g0) %>%
  add_clusters(links=c0) +
  geom_seq() +
  geom_gene() +
  geom_link() +
  geom_gene_tag(aes(label=feat_id)) 
p1

# zoom in on loci with links, group links less than 200 apart, show 200 flank
# you may have to play with those values
p1 %>% 
  focus(.track_id=links, .max_dist = 200, .expand=200) + 
  geom_seq_label()

# and add nice locus labels that indicate start/end of regions
p1 %>% 
  focus(.track_id=links, .max_dist = 200, .expand=200,
    .locus_id=str_glue("{seq_id}:{start}-{end}")) + 
  geom_seq_label()

# scaling up B 25 times - not sure how much sense this representation makes, though...
g_scaled <- g0 %>%
  mutate(across(c(start, end), ~ifelse(str_detect(seq_id, "B"), .x * 25, .x)))

gggenomes(g_scaled) %>% add_clusters(c0) + geom_seq() + geom_gene() + geom_link()