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

Put ./. genotype when there is no power to identify a variant #137

Closed
mfoll opened this issue Oct 20, 2016 · 1 comment
Closed

Put ./. genotype when there is no power to identify a variant #137

mfoll opened this issue Oct 20, 2016 · 1 comment
Assignees

Comments

@mfoll
Copy link
Member

mfoll commented Oct 20, 2016

At the moment we put 0/0 when QVAL<50, but this can be due to the absence of mutation but also absence of power to detect a mutation. This should now be easy now that we implemented #133.

We need to have a parameter to choose the allelic fraction we are trying to identify:

  • For germline variants we should use the same parameter use to classify somatic/germline status (20% at the moment).
  • For somatic mutation this will be a new parameter.
@tdelhomme
Copy link
Member

@nalcala nextflow part seems to be ok, I've also added lines in the R script to recover the values, keeping your previous synthax (afmin and sigma).

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

No branches or pull requests

3 participants