plot_arcs_compare
Opened this issue · 1 comments
region argument in plot_arcs_compare does not behave as expected. This argument still displays the entire sequence x-axis but only plots interactions that begin prior to the region cutoff. Even if the interaction extends beyond the region cutoff, the entire arc is displayed. For example:
plot = rnav.plot_arcs_compare(
samples=[6kb_RNA1, 6kb_RNA2],
sequence="ss",
plot_kwargs={"dpi": 50},
interactions="pairprob",
region=[1,2000]
)
returns a 6kb x axis with the full sequence displayed. Arcs exist for the first 2kb. Some arcs begin before nucleotide 2000 but extend to nucleotide 2200 and are still fully displayed.
On another note - it would also be very useful to be able to plot annotations on these graphs. They currently don't accept that argument.