Skip to content

Prevent Pyrodigal from adding stop codon indicators #447

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

Merged
merged 7 commits into from
Feb 6, 2025
Merged
Show file tree
Hide file tree
Changes from 2 commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
1 change: 1 addition & 0 deletions CHANGELOG.md
Original file line number Diff line number Diff line change
Expand Up @@ -23,6 +23,7 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
- [#435](https://github.com/nf-core/funcscan/pull/435) Fixed dependency errors within taxonomy merging scripts, updated the code and output for all three workflows. Bumped to version 0.1.1. (by @darcy220606)
- [#437](https://github.com/nf-core/funcscan/pull/437) Fixed file name error when supplying already preprocessed CARD database for ARG workflow. (by @jasmezz)
- [#446](https://github.com/nf-core/funcscan/pull/446) Updated antiSMASH modules to fix apptainer execution. (by @jasmezz and @jfy133)
- [#447](https://github.com/nf-core/funcscan/pull/447) Added `--annotation_pyrodigal_stopcodon` parameter to improve AMPlify screening. (by @jasmezz)

### `Dependencies`

Expand Down
1 change: 1 addition & 0 deletions conf/modules.config
Original file line number Diff line number Diff line change
Expand Up @@ -201,6 +201,7 @@ process {
params.annotation_pyrodigal_singlemode ? "-p single" : "-p meta",
params.annotation_pyrodigal_closed ? "-c" : "",
params.annotation_pyrodigal_forcenonsd ? "-n" : "",
params.annotation_pyrodigal_stopcodon ? '--no-stop-codon' : '',
Copy link
Member

Choose a reason for hiding this comment

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

The pipeline flag should also be nostopcodon so it's clearer to the user

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

Okay, is more intuitive yes. On the other hand I wanted to stick to the "codex" to have params default to false in the nextflow.config (the renamed annotation_pyrodigal_nostopcodon would have to be true by default).

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

Oh I just saw that I flipped the logic by mistake 🙈 I'll update the above line to
params.annotation_pyrodigal_stopcodon ? '' : '--no-stop-codon',

"-g ${params.annotation_pyrodigal_transtable}",
].join(' ').trim()
}
Expand Down
1 change: 1 addition & 0 deletions nextflow.config
Original file line number Diff line number Diff line change
Expand Up @@ -66,6 +66,7 @@ params {
annotation_pyrodigal_closed = false
annotation_pyrodigal_transtable = 11
annotation_pyrodigal_forcenonsd = false
annotation_pyrodigal_stopcodon = false

annotation_bakta_db = null
annotation_bakta_db_downloadtype = 'full'
Expand Down
7 changes: 7 additions & 0 deletions nextflow_schema.json
Original file line number Diff line number Diff line change
Expand Up @@ -518,6 +518,13 @@
"fa_icon": "fas fa-barcode",
"description": "Forces Pyrodigal to scan for motifs.",
"help_text": "Forces Pyrodigal to a full scan for motifs rather than activating the Shine-Dalgarno RBS finder, the default scanner for Pyrodigal to train for motifs.\n\nFor more information check the Pyrodigal [documentation](https://pyrodigal.readthedocs.io).\n\n> Modifies tool parameter(s):\n> - PYRODIGAL: `-n`"
},
"annotation_pyrodigal_stopcodon": {
"type": "boolean",
"fa_icon": "fa fa-star-of-life",
"description": "This forces Pyrodigal to append asterisks (`*`) as stop codon indicators.",
"help_text": "Some downstream tools like AMPlify cannot process sequences containing non-sequence characters like the stop codon indicator `*`. Thus, this flag is deactivated by default. Activate this flag to revert the behaviour and have Pyrodigal append `*` as stop codon indicator to annotated sequences.\n\nFor more information check the Pyrodigal [documentation](https://pyrodigal.readthedocs.io).\n\n> Modifies tool parameter(s):\n> - PYRODIGAL: `--no-stop-codon`",
"default": false
}
},
"fa_icon": "fas fa-file-signature"
Expand Down