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

Calculating unmutated frequency in IgDiscover output #4

Open
scharch opened this issue Feb 14, 2024 · 0 comments
Open

Calculating unmutated frequency in IgDiscover output #4

scharch opened this issue Feb 14, 2024 · 0 comments

Comments

@scharch
Copy link

scharch commented Feb 14, 2024

According to the docs here, orgdbstats should be run on the unzipped filtered.tsv. However, the sequence_alignment field output by IgDiscover is ungapped, which results in ogrdbstats assigning 0 unmutated_sequences for all alleles, even though IgDiscover itself thinks that about 20% of the sequences are unmutated. Any way around this?

This is a non-model species. I've manually gapped the "reference" (such as it is) database, but I can't do that for the bulk sequencing data and I'm not finding an annotation tool that does gaps and can handle non-model species.

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

1 participant