-
Notifications
You must be signed in to change notification settings - Fork 6
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
Generate PoN for somatic SVs? #125
Comments
It's something to try, but I'm not sure if we have a very good sense of this: gnomAD-SV paper: https://www.biorxiv.org/content/10.1101/578674v1 My current thought is that this will require a great deal of investigation by CCS before we implement something. EDIT: The larger question is the following: how do you precisely implement such a filter? Let's consider SNVs. Let's say we're confident (based on the panel of normals) that a germline SNP at position X is a germline SNP, and should not be considered a somatic SNP. Ergo, it should be filtered. That's a pretty easy problem. Now let's consider SVs. SVs are denoted by two breakpoints and an event annotation. Based on a hypothetical "SV PoN", let's say we were 100% confident that an SV which is an INS at The outputs from SV callers will almost never show that. The variant could be mislabeled, or the breakpoints at are But then the question is, how do you incorporate the inherent "awfulness" of standard short-read SV callers? If you create filters which filter exactly for the events such as "INS at bkpt1=X and bkpt2=Y", I'm pretty sure you'll find nothing. Best case scenario could be a few correctly filtered event, worse case scenario would be a massive mess. I'm also not convinced the world has analyzed a high enough sample size to say that a given SV is a "standard" SV (with the exception of a handful of cases), but that's another kettle of fish**. **kettle of fish for the fiskgryta and fisksoppa, of course |
We'll need to devote science time to this @kpjonsson |
See #126, but for SVs. This is in addition to filtering against gnomAD.
The text was updated successfully, but these errors were encountered: