diff --git a/CHANGELOG.md b/CHANGELOG.md index f0d2e01..b130143 100644 --- a/CHANGELOG.md +++ b/CHANGELOG.md @@ -8,20 +8,22 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0 ### `Added` - Functionality to create DROP databases and to add samples to existing ones [#147](https://github.com/genomic-medicine-sweden/tomte/pull/147) +- A switch `--skip_variant_calling` for variant calling [#168](https://github.com/genomic-medicine-sweden/tomte/pull/168) ### `Fixed` ### `Parameters` -| Old parameter | New parameter | -| -------------------------------- | ------------------- | -| `--max_cpus` | | -| `--max_memory` | | -| `--max_time` | | -| `--validationShowHiddenParams` | | -| `--validationSkipDuplicateCheck` | | -| `--validationS3PathCheck` | | -| `--monochromeLogs` | `--monochrome_logs` | +| Old parameter | New parameter | +| -------------------------------- | ------------------------ | +| `--max_cpus` | | +| `--max_memory` | | +| `--max_time` | | +| `--validationShowHiddenParams` | | +| `--validationSkipDuplicateCheck` | | +| `--validationS3PathCheck` | | +| `--monochromeLogs` | `--monochrome_logs` | +| | `--skip_variant_calling` | ### `Changed` diff --git a/conf/modules/call_variants.config b/conf/modules/call_variants.config index 5ae2895..f2ff28e 100644 --- a/conf/modules/call_variants.config +++ b/conf/modules/call_variants.config @@ -16,6 +16,7 @@ process { withName: '.*CALL_VARIANTS.*' { + ext.when = { (!params.skip_variant_calling) } publishDir = [ enabled: false, ] @@ -60,7 +61,7 @@ process { saveAs: { filename -> filename.equals('versions.yml') ? null : filename }, ] } - withName: '.*CALL_VARIANTS:SPLIT_MULTIALLELICS' { + withName: '.*CALL_VARIANTS:SPLIT_MULTIALLELICS' { ext.prefix = { "${meta.id}_norm" } ext.args = { '--output-type z --multiallelics -both --write-index=tbi' } } diff --git a/conf/test.config b/conf/test.config index a96e542..b64930b 100644 --- a/conf/test.config +++ b/conf/test.config @@ -26,6 +26,7 @@ params { // Skip when GITHUB actions skip_drop_ae = System.getenv("GITHUB_ACTIONS").equals(null) ? false : true skip_drop_as = System.getenv("GITHUB_ACTIONS").equals(null) ? false : true + skip_variant_calling = true // Input data input = "${projectDir}/test_data/samplesheet_chr21.csv" diff --git a/docs/usage.md b/docs/usage.md index 2d157fe..1c305de 100644 --- a/docs/usage.md +++ b/docs/usage.md @@ -190,9 +190,11 @@ The mandatory and optional parameters for each category are tabulated below. | --------- | -------------------------------- | | | variant_caller1 | | | bcftools_caller_mode2 | +| | skip_variant_calling3 | 1 If it is not provided by the user, the default value is bcftools
-2 If it is not provided by the user, the default value is multiallelic +2 If it is not provided by the user, the default value is multiallelic
+3 If it is not provided by the user, the default value is false #### 5. SNV annotation (ensembl VEP) @@ -284,7 +286,7 @@ To build your own database you will need at least 50 for aberrant expression, if - `--skip_drop_ae false` if you want to get a database for aberrant expression - `--skip_subsample_region false` if you have sequenced any material with overrepresented regions (such as hemoglobin in whole blood) we recommend to remove it by setting this parameter to false and providing a bed with the overrepresented region with `--subsample_bed` - `--skip_downsample false` if you have very deeply sequenced samples, we recommend to downsample, the default is 60M read pairs -- `--skip_build_tracks true`, `--skip_stringtie true`, `--skip_vep true` as most users will be interested in getting the database rather than other downstream results +- `--skip_build_tracks true`, `--skip_stringtie true`, `--skip_variant_calling`, `--skip_vep true` as most users will be interested in getting the database rather than other downstream results Running DROP with many samples requires a lot of time and a lot of memory, that is why we recommend to subsample overrepresented regions and downsample if you have deeply sequenced samples. If your run fails for either of this reasons, try to relaunch it from the work directory where DROP was run so that DROP continues from the point where it failed (if you restart the pipeline with `-resume` it will begin from the start and it will likely fail in the same way). diff --git a/nextflow.config b/nextflow.config index 1f738bd..2825015 100644 --- a/nextflow.config +++ b/nextflow.config @@ -44,6 +44,7 @@ params { save_mapped_as_cram = true // Variant calling + skip_variant_calling = false variant_caller = 'bcftools' bcftools_caller_mode = 'multiallelic' skip_build_tracks = false diff --git a/nextflow_schema.json b/nextflow_schema.json index 03b80e7..00499ae 100644 --- a/nextflow_schema.json +++ b/nextflow_schema.json @@ -273,6 +273,12 @@ "default": "", "fa_icon": "fas fa-filter", "properties": { + "skip_variant_calling": { + "type": "boolean", + "default": false, + "description": "Skip variant calling for all samples.", + "fa_icon": "fas fa-toggle-off" + }, "variant_caller": { "type": "string", "default": "bcftools",