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

161 cnv benchmark subworkflow #163

Merged
merged 10 commits into from
Mar 7, 2025
Merged

Conversation

kubranarci
Copy link
Contributor

  • new analysis type especially for CNV benchmarking. Enables intersection of BED regions using bedtools.
  • truvari works for cnv benchmarking only if pctseq = 0 --> means no sequence resolve needed.
  • new test case for cnv benchmarking and snapshots
  • intersect added to the metromap
  • intersection will be possible my target_regions input, if not provided VCF can be converted to BED for intersection analysis.

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/variantbenchmarking branch on the nf-core/test-datasets repository.
  • Make sure your code lints (nf-core pipelines lint).
  • Ensure the test suite passes (nextflow run . -profile test,docker --outdir <OUTDIR>).
  • Check for unexpected warnings in debug mode (nextflow run . -profile debug,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).

@kubranarci kubranarci linked an issue Mar 6, 2025 that may be closed by this pull request
Comment on lines +19 to +25
test.branch{
def vcf_file = it[1]
def regions_file = it[2]
vcf : vcf_file
regions : regions_file
other: false}
.set{test_samples}
Copy link
Contributor

Choose a reason for hiding this comment

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

What's this layout? 😁

Copy link
Contributor Author

Choose a reason for hiding this comment

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

I have two types of inputs, one is vcf the other is regions bed file. This stream should branch them accordingly. First check assures if vcf file is given, second is to check if regions exist :D it is a bit messy I know. But let me know if you have something better on your mind

Copy link
Contributor

@nvnieuwk nvnieuwk left a comment

Choose a reason for hiding this comment

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

Some minor comments, you can merge this in once you've addressed these :)

@kubranarci kubranarci merged commit 2d2773f into dev Mar 7, 2025
13 checks passed
@kubranarci kubranarci deleted the 161-cnv-benchmark-subworkflow branch March 7, 2025 10:18
@kubranarci kubranarci linked an issue Mar 7, 2025 that may be closed by this pull request
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.

CNV benchmark subworkflow fix broken link describing truth sets
2 participants