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: gcnv contig wrapper scripts do not clean their output folders (#443) #468

Conversation

Nicolai-vKuegelgen
Copy link
Contributor

Some of the gcnv step_parts and their associated wrappers create collections of output folders that are currently not tracked by snakemake. Additionally gcnv does only partially overwrite these output files & folders, especially if it is run with fewer samples than in a previous run.

This changes the case_mode wrappers for contig_ploidy and call_cnvs step-parts to always remove the directories they will write into to ensure all output files belong to the current snappy run.

This could be extended to the cohort_mode wrappers, however they might be used in cases where removal of previous samples is not intended, since they could still be relevant (whereas in case mode this is not the case).

Fixes #443

…etermineGermlineContigPloidy runs (instead of skipping them)

- always purge previous output of GCNV call_cnvs_case_mode step when it's run
…results. This is the cleanest (and least error prone) way of stopping any intermixing with previous results.
@coveralls
Copy link

coveralls commented Nov 21, 2023

Coverage Status

coverage: 85.915%. remained the same
when pulling 0feeb78 on 443-gcnv_contig_ploidy-wrapper-script-crahses-because-gatk-determinegermlinecontigploidy-output-is-not-cleaned-but-always-fully-read
into 50d0781 on main.

Copy link
Member

@holtgrewe holtgrewe left a comment

Choose a reason for hiding this comment

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

LGTM

@Nicolai-vKuegelgen Nicolai-vKuegelgen merged commit 618f10d into main Nov 29, 2023
8 checks passed
@Nicolai-vKuegelgen Nicolai-vKuegelgen deleted the 443-gcnv_contig_ploidy-wrapper-script-crahses-because-gatk-determinegermlinecontigploidy-output-is-not-cleaned-but-always-fully-read branch November 29, 2023 15:21
@tedil tedil mentioned this pull request Jun 28, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
3 participants