Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Support Microsynteny plots with genome alignment data #547

Open
8 tasks
Adamtaranto opened this issue Feb 18, 2023 · 3 comments
Open
8 tasks

Support Microsynteny plots with genome alignment data #547

Adamtaranto opened this issue Feb 18, 2023 · 3 comments
Assignees

Comments

@Adamtaranto
Copy link
Collaborator

Adamtaranto commented Feb 18, 2023

jcvi.graphics.ribbon module is a modified version of the synteny module which allows plotting of arbitrary alignments and multiple genomic feature tracks. Ribbon plots are not based on ortholog linkages but instead on whole genome alignments.

This issue is for tracking progress on the ribbon module in preparation for merging it into the main project.

  • Modify LASTZ wrapper to do pairwise alignments between a defined set of contigs
  • Add filter for alignment length and quality
  • Convert LASTZ output to .bed and .blocks files with wga2ribbon app
  • Extend wga2ribbon to use Minimap2 asm2asm alignment format
  • Modify synteny module to plot arbitrary alignments
  • Add logging and comment new functions
  • Check that we are using pre-existing jcvi functions where possible
  • Document use cases in new wiki page
@AlbertoBernacchi
Copy link

Hello Adamataranto,
have you managed to modify the synteny module to plot arbitrary alignments? Can you point me to a bioinformatic package that is able to do so? I'm currently searching for it but with no luck.
Thank you.

@Adamtaranto
Copy link
Collaborator Author

I have a development branch that works for this. It is currently undocumented and you will need to provide your own alignments.

I can provide some support if you want to give it a go.

@AlbertoBernacchi
Copy link

Hi Adamtarando,
thank you for the quick response. I'll try to explain what I want to do and what I can currently do with the data I have. I have two mitochondrial genomes of which I want to analyze the collinearity. However, I have yet to find a tool that allows for the analysis of non-coding regions. Currently I have the .bed files, .gff, .fasta and .blast of both mitochondrial genomes.
Thank you again for your help.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

2 participants