Skip to content
This repository has been archived by the owner on Mar 17, 2023. It is now read-only.

Consider using "sourmash search --containment" for "translate" #86

Open
olgabot opened this issue Aug 3, 2020 · 1 comment
Open

Consider using "sourmash search --containment" for "translate" #86

olgabot opened this issue Aug 3, 2020 · 1 comment

Comments

@olgabot
Copy link
Contributor

olgabot commented Aug 3, 2020

Currently, sencha translate uses a simple match of whether 100% of the k-mers from the reading frame, match the reference proteome. But a "Franken k-mer" situation can happen where the reading frame has 100% match, but the k-mers are all from different genes. Using sourmash search --containment would only search for "consecutive" k-mers that all appear in a single gene (or maybe family of genes??) and would be an improvement over the current method.

image

Thanks to @bluegenes for the idea!

@olgabot
Copy link
Contributor Author

olgabot commented Aug 3, 2020

Maybe use sourmash lca search --containment by hacking the Least Common Ancestor (LCA) code to use gene/protein families instead of individual genes, so that the matches could be on a per-gene family level.

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

No branches or pull requests

1 participant