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

How to change mapping parameters when dealing with highly repetitive region #82

Open
ghost opened this issue Nov 28, 2017 · 0 comments
Open

Comments

@ghost
Copy link

ghost commented Nov 28, 2017

We have done a genome assembly of a bacterium based on Oxford Nanopore data. The assembly looks very good overall, but I have one region where I want to confirm the assembly via mapping the longest 20% of the ONP reads back to the assembly.

The problem is that this region is very repetitive. I guess you could call it long tandem repeats of several hundred base pairs.

I think the assembly is correct, but graphmap is probably having trouble there. To me it looks like it is finding a seed on one particular repeat region of the refercence but the correct position would actually be shifted one/some repeat(s) left/right. Then the extension works but is "out of phase" in terms of repeats. When I have a look at the repeats on IGV I have insertions an deletions all about the same size (probably corresponding to the length of the tandem repeat).

The reason why I think that the assembly is correct is because I have about 30 to 40% of the long reads spanning the region correctly without any long indel.

How can I tweak the parameters of graphmap to avoid this problem?

Or do you think my hypothesis is wrong?

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

No branches or pull requests

0 participants