Open cmdcolin opened 3 years ago
carried over from https://github.com/GMOD/jbrowse-components/issues/1133
some R code to diagonalize here http://jmonlong.github.io/Hippocamplus/2017/09/19/mummerplots-with-ggplot2/
@cmdcolin Thanks for your suggestions and sharing the R code link. I tried the R code above. Unfortunately it seems not working with multiple chromosomes. Do you have code that can generate diagonal plot between genome reference of multiple chromosomes and assembled contigs of same species? Thanks!
here's the example code for mummerplot that takes a single sequence from each reference and target and require the start and end range of that sequences
mummerplot -x "[0,275287]" -y "[0,265111]" -postscript -p mummer mummer.mums
See for example
This is a human-to-human alignment loaded as a PAF file
Mummerplot automatically diagonalizes the plot and flips backwards scaffolds to create a by-default better display of the same data
We may be able to accomplish something similar, would require some computational tricks that probably exist in the mummer codebase
Being able to allow the user to manually flip a region, or re-order the regions, may be interesting as well