-
Notifications
You must be signed in to change notification settings - Fork 53
Vcf2cytosure #456
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
Vcf2cytosure #456
Conversation
|
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Great to see that it has been added! :D
Just one small thing, I would change the name of the switch to skip and set it to true by default so that it fits better the rest of the code.
nextflow_schema.json
Outdated
}, | ||
"vcf2cytosure_switch": { | ||
"type": "boolean", | ||
"description": "Turn on the vcf2cytosure subworkflow", | ||
"fa_icon": "fas fa-toggle-on", | ||
"help_text": "vcf2cytosure can generate CGH files from a structural variant VCF file that can be analysed in the CytoSure interpretation software. Cut offs for allele frequencies and bin sizes can be modified in the config file. Turned off by default." |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Instead of switch, I would call it something along the lines of skip_vcf2cytosure and have it set to true by default, that way it would be similar to hoe the rest of the switches are build.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I know what you mean, it would be great to keep the naming consistent in the pipeline. I had it exactly like you suggested at first but ran into issues with turning turning the sub-workflow on using the command line. Since the flag was set to true by default, giving the option --skip_vcf2cytosure
will just skip it again.
I will try to find a solution
Co-authored-by: Lucía Peña-Pérez <lucia.pena.perez@scilifelab.se>
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good! :D
|
||
withName: '.*GENERATE_CYTOSURE_FILES:VCF2CYTOSURE' { | ||
ext.args = { [ | ||
meta.sex.equals('1') ? '--sex male' : '--sex female', |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
If meta.sex
is set to anything else than 1 it will default to female. I believe that is the general default of the program.
|
||
withName: '.*GENERATE_CYTOSURE_FILES:VCF2CYTOSURE' { | ||
ext.args = { [ | ||
meta.sex.equals('1') ? '--sex male' : '--sex female', |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
If meta.sex
is set to anything else than 1 it will default to female. I believe that is the general default of the program.
vcf2cytosure is a tool to convert SV vcf:s to a format that can be loaded into the CytoSure interpretation software. This PR includes a subworkflow for that tool. Since the workflow and tool is fairly specific to our setup it is turned off by default. It is turned on by supplying the option
--vcf2cytosure_switch
PR checklist
nf-core lint
).nextflow run . -profile test,docker --outdir <OUTDIR>
).nextflow run . -profile test_one_sample,docker --outdir <OUTDIR>
).docs/usage.md
is updated.docs/output.md
is updated.CHANGELOG.md
is updated.README.md
is updated (including new tool citations and authors/contributors).