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

Create a line plot for synteny #74

Closed
GallVp opened this issue Aug 9, 2023 · 1 comment · Fixed by #89
Closed

Create a line plot for synteny #74

GallVp opened this issue Aug 9, 2023 · 1 comment · Fixed by #89
Assignees
Milestone

Comments

@GallVp
Copy link
Member

GallVp commented Aug 9, 2023

No description provided.

@GallVp
Copy link
Member Author

GallVp commented Aug 11, 2023

From Ignacio Carvajal

It is apparently possible to run Dgenies locally on linux. I have not played around with it too much because the nice thing about it is the interactivity of the results. I'll experiment running this on powerplant at some point and opening the results on powerplant from a browser. It would be nice to skip downloading on the local machine step which can be a pain, especially for larger geneomes. There are other dotplot visualization packages out there if you run into trouble with dgenies.

I think this is what you mean by filtering of the matches based on identity or length. See here for an example:
https://dgenies.toulouse.inra.fr/result/chimp_vs_human

Generating the PAF file that is requried as an input is as easy as:
'minimap2 -t 10 -cx asm10 ${ANITRA_ASSEMBLY} ${WA_ASSEMBLYA}> AN_WA23_asm10.paf'

This is comparing two raspberry assemblies using 10 threads. It takes 150 sec to run. The other inputs for Dgeneies are the two assemblies (${ANITRA_ASSEMBLY} and ${WA_ASSEMBLYA} in this case).

One important thing that I have not tested is the minimap2 parameters. I've been using asm10, which they say is for "intra-species asm-to-asm alignments where the expected divergence is 1%". You can also specify asm5 and asm20 which correspond to .5% and 2% respectively. 1% seems to work well in my cases but it may not be a one size fits all.

Please keep me in the loop for your ideas and if you need help/testing. Keen to contribute.

@GallVp GallVp added the v2 label Aug 13, 2023
@GallVp GallVp added this to the v1.3.1 milestone Apr 16, 2024
@GallVp GallVp self-assigned this Apr 16, 2024
@GallVp GallVp added feature request and removed v2 labels Apr 16, 2024
@GallVp GallVp modified the milestones: v1.3.1, 2.0 Apr 16, 2024
@GallVp GallVp modified the milestones: v1.3.2, 2.0.0 May 3, 2024
@GallVp GallVp mentioned this issue May 30, 2024
10 tasks
@GallVp GallVp closed this as completed in #89 Jun 5, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging a pull request may close this issue.

1 participant