[go: nahoru, domu]

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

Ns in core gene alignment #366

Closed
karchern opened this issue Nov 29, 2017 · 6 comments
Closed

Ns in core gene alignment #366

karchern opened this issue Nov 29, 2017 · 6 comments
Labels

Comments

@karchern
Copy link

Hello,

I occasionally have a stretch of Ns in my core gene alignment. I don't think this makes any sense. Has this issue occured before? I tried searching past issues but didn't find anything.

Thank you.

@andrewjpage
Copy link
Member

Hi,
If a core gene is absent you will get a stretch of N's. This is because by default the definition of a core gene is that it must be in 99% of isolates, so for that 1% you need to fill in the gap with something to make it all align.

@karchern
Copy link
Author

Hi Andrew,
thanks for the answer. I always assumed roary filled up missing sites with gaps (- character) by default. I assume it uses "-" when part of a gene is missing and Ns otherwise?

@andrewjpage
Copy link
Member

Sorry I forgot I changed it recently, its the other way around, a dash is used for a missing gene.

@karchern
Copy link
Author

I'm using roary 3.8.2. I assume this is not the most recent version. Do you know in which version you flipped? From visual inspection of my core gene alignment, I strongly believe that Ns correspond to entire missing genes and gaps correspond to missing bits of a gene.

Thank you for the clarification!

@andrewjpage
Copy link
Member

It was changed in the last release:
https://github.com/sanger-pathogens/Roary/releases

@karchern
Copy link
Author

Great, thanks for the quick response!

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

No branches or pull requests

2 participants