Add WDL implementation for Salmon - #326
Conversation
✅ Snyk checks have passed. No issues have been found so far.
💻 Catch issues earlier using the plugins for VS Code, JetBrains IDEs, Visual Studio, and Eclipse. |
|
This also needs a |
c2d9dde to
4f222d5
Compare
a-frantz
left a comment
There was a problem hiding this comment.
This is looking really great! Left some comments to address, but this is close to the finish line 🚀
_Describe the problem or feature in addition to a link to the issues._ This is needed for writing decent tests in #326 BCR and ABL1 were chosen somewhat arbitrarily. We do have other test fixtures referencing these genes specifically (tests are for Arriba), but this doesn't interact with those tests v50 release chosen as its latest. Other fixtures use v31; I could use that release if we want Before submitting this PR, please make sure: - [x] You have added a few sentences describing the PR here. - [x] The code passes all CI tests without any errors or warnings. - [ ] You have added tests (when appropriate). - [ ] You have added an entry in any relevant CHANGELOGs (when appropriate). - [ ] If you have made any changes to the `scripts/` or `docker/` directories, please ensure any image versions have been incremented accordingly! - [x] You have updated the README or other documentation to account for these changes (when appropriate).
…d FASTQ data, add SE test
…add decoy-aware indexing test
71e9c9a to
6bc4b66
Compare
|
since I iterated on the index fixture several times (custom → v1 → v2 → reverted → BCR/ABL1), my branch's LFS history now carries a few superseded large blobs (~1.3GB total across old versions) even though the current file is tiny. If you squash-merge, this shouldn't carry into main's history — let me know if you'd prefer I try to clean this up |
we'll do a squash-merge 👍 no action needed on your end |
a-frantz
left a comment
There was a problem hiding this comment.
Good news: we're at the nitpicking stage of the review 🤣 Left a bunch of comments, but all of them are about very minor things. This looks great 🚀
…upport, misc review fixes
|
Just a quick note - builds_index_with_decoys fails locally for me on Windows with Unsupported reparse point type when reading back a symlink created inside the container. Manually confirmed the symlink itself is created correctly on the Linux side (via docker run ... ln -sf ...), so this looks like a Windows-host/Docker-volume limitation on my end, not an actual bug — should be fine in your CI, but flagging in case you want to verify. |
|
@PriyankaaXD I just pushed some small nits I wanted in before merge but didn't want to force us through another round of back and forth for such minor points. Hope you don't mind! @adthrasher could you give this an additional review? |
|
Thanks for testing at production scale — good catch that decoys need to factor into the memory calc more heavily. Happy to update the formula to include decoys_fasta_size once you've got final numbers, or if you'd rather just push the tuned value yourself given you already have the real data, that works too! |
adthrasher
left a comment
There was a problem hiding this comment.
Finished quant runs for reference with and without decoys.
#0 completed wall 338m00s queued 30.6s 4 cpu, 31.7 GiB peak 10.9 GiB, cpu 1337m03s
#0 completed wall 360m59s queued 31.5s 4 cpu, 10.2 GiB peak 4.3 GiB, cpu 1432m16s
|
Thanks @PriyankaaXD for implementing this! |
Adds a WDL implementation for Salmon (mapping-mode only), per the "tool wishlist" issue #228.
Tasks added in
tools/salmon.wdl:index— wrapssalmon indexquant— wrapssalmon quantAll "important options" from Salmon's docs are exposed as inputs, with defaults verified against
salmon quant --help-readsoutput on Salmon 2.6.0. Parameter documentation is copied from Salmon's official docs, per guidance in the issue. Scoped to mapping-mode only (FASTQ input) — no BAM/alignment-mode support, as requested.Tests added in
tools/salmon.ymlusing the new Sprocket test framework, covering both tasks with real output assertions. Verified locally:sprocket lintpasses cleanly,sprocket dev testpasses both tests.This is my first contribution to this project — happy to make any adjustments you'd like!
Before submitting this PR, please make sure:
scripts/ordocker/directories, please ensure any image versions have been incremented accordingly!