brennanpincardiff / drawProteins

Creating package to draw proteins from Uniprot API
Other
33 stars 19 forks source link

Error in draw_chains function, incomplete representation #33

Open MauriAndresMU1313 opened 2 years ago

MauriAndresMU1313 commented 2 years ago

I´m trying to show 7 proteins of different CoVs, 4 of these proteins (red arrows) show only part of the protein S2 (subunit 2). I don´t understand why this happens. I reviewed if the protein in Uniprot is complete and it looks like that all is fine. Here is my script:

data <- drawProteins::get_features("P0DTC2 A0A6B9WHD3 A0A6M3G9R1 P59594 E0XIZ3 K9N5Q8 U5NJG5")

p <- draw_canvas(real_data)
p <- draw_chains(p, real_data, label_size =  2.5)
p <- draw_domains(p, real_data, show.legend = F, label_domains = F)
p <- draw_motif(p, real_data)

p <- p + theme_bw(base_size = 10) +
  theme(panel.grid.minor = element_blank(),
        panel.grid.major = element_blank()) +
  theme(axis.ticks = element_blank(),
        axis.text.y = element_blank()) +
  theme(panel.border = element_blank())

So, clearly, there is a issue in these protein draws, that only shows the S2. There is a way to force the code to show the complete protein?

Here is the plot

Screenshot_1