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

SNPgenie for one bacterial strain #40

Open
SolayMane opened this issue Nov 9, 2020 · 1 comment
Open

SNPgenie for one bacterial strain #40

SolayMane opened this issue Nov 9, 2020 · 1 comment

Comments

@SolayMane
Copy link

Dear Chase,
I have a vcf file for snps of one strain against a ref. Can I use SNPgenie to calculate Pi trough the whole genome using sliding window?
Sincerely

@singing-scientist
Copy link
Contributor

Greetings @SolayMane ! Yes you can. If you want to do it yourself for a certain window, πN can be calculated as sum(N_diffs) / sum(N_sites), where N=nonsynonymous, and πS in the same way using S (for synonymous). But make sure the windows are not too small, or they will be subject to wild random fluctuations! Let me know if that helps.

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

2 participants