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

Fix input channels for variantrecalibrator #1022

Merged
merged 5 commits into from
May 16, 2023

Conversation

asp8200
Copy link
Contributor

@asp8200 asp8200 commented May 14, 2023

Fix input channels for variantrecalibrator in the joint-germline workflow.

There is currently no pytest testing that this works. I tested it manually with downsampled versions of the giab-samples HG002 and HG003.

PR checklist

  • This comment contains a description of changes (with reason).
  • If you've fixed a bug or added code that should be tested, add tests!
  • If you've added a new tool - have you followed the pipeline conventions in the contribution docs- [ ] If necessary, also make a PR on the nf-core/sarek branch on the nf-core/test-datasets repository.
  • Make sure your code lints (nf-core lint).
  • Ensure the test suite passes (nextflow run . -profile test,docker --outdir <OUTDIR>).
  • Usage Documentation in docs/usage.md is updated.
  • Output Documentation in docs/output.md is updated.
  • CHANGELOG.md is updated.
  • README.md is updated (including new tool citations and authors/contributors).

@asp8200 asp8200 requested a review from maxulysse as a code owner May 14, 2023 18:50
@FriederikeHanssen
Copy link
Contributor

Changelog :)

@github-actions
Copy link

github-actions bot commented May 14, 2023

nf-core lint overall result: Passed ✅

Posted for pipeline commit c026ee3

+| ✅ 153 tests passed       |+
#| ❔   9 tests were ignored |#

❔ Tests ignored:

✅ Tests passed:

Run details

  • nf-core/tools version 2.8
  • Run at 2023-05-16 08:48:05

@asp8200
Copy link
Contributor Author

asp8200 commented May 14, 2023

"Fix" doesn't work when either known_snps_vqsr, known_indels_vqsr or dbsnp_vqsr isn't set.

dbsnp_vqsr = params.dbsnp_vqsr ?: Channel.empty()
known_indels_vqsr = params.known_indels_vqsr ?: Channel.empty()
known_snps_vqsr = params.known_snps_vqsr ?: Channel.empty()
dbsnp_vqsr = params.dbsnp_vqsr ? Channel.value(params.known_snps_vqsr) : Channel.empty()
Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@FriederikeHanssen and @maxulysse : The joint-germline genotyping now seems to work both when known_snps_vqsr etc. are set and when they are not.

Should we do a similar initialisation of the other channels in this section?

@amizeranschi
Copy link
Contributor

Is there any way I could test the changes from this PR, to see if they solve #1016 ?

@maxulysse
Copy link
Member

Yes, you can run:

nextflow run asp8200/sarek -r fix_variantrecalibrator ...

@FriederikeHanssen
Copy link
Contributor

that would be great @amizeranschi to test that it all works now :)

@amizeranschi
Copy link
Contributor

I tested this now and can confirm that this PR fixes #1016.

@maxulysse
Copy link
Member

Changelog :)

Just need the changelog and we're good to go

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

Successfully merging this pull request may close these issues.

Error running VQSR in joint VC: nextflow.exception.MissingProcessException: Missing process or function mix
4 participants