-
Notifications
You must be signed in to change notification settings - Fork 33
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
problem with expand_no #83
Comments
Hi, Yes, you're correct. This sounds like it might be a bug. I will try and take a look this week. |
Thanks very much - I can provide a sample dataset for diagnostics if helpful. |
The problem appears to happen on one side of the gene - so the script behaves as expected to the left of the gene but not the right: e.g. where expand_no -= 5 and gene = aac6IIc~~~aac61bcr~~~~~~aacA4 |
Hi, Apologies, this took me a bit longer to get to than expected. I am having trouble reproducing the issue. It would be great if you could provide the sample dataset that leads to the problem. |
My understanding of expan_no in panatoo-gene-neighbourhood is that it should give N genes up and downstream of the target gene X. In my version of panaroo (from conda - version 1.2.3) - this works for --expand_no 1 but not for --expand_no 2 where some have many more genes than 2 up or downstream
Do you have any idea why this might be?
The text was updated successfully, but these errors were encountered: