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

SALMON_INDEX runs when using --aligner star_rsem even if samples have explicit strandedness #975

Closed
drpatelh opened this issue Mar 24, 2023 · 1 comment
Labels
bug Something isn't working
Milestone

Comments

@drpatelh
Copy link
Member

drpatelh commented Mar 24, 2023

Description of the bug

nextflow run main.nf \
    --input https://raw.githubusercontent.com/nf-core/test-datasets/rnaseq/samplesheet/samplesheet_test.csv \
    --fasta 'https://github.com/nf-core/test-datasets/raw/rnaseq/reference/genome.fasta' \
    --gtf 'https://github.com/nf-core/test-datasets/raw/rnaseq/reference/genes.gtf' \
    -profile docker \
    --aligner star_rsem \
    --outdir results \

Even though we have specified reverse explicitly for all samples in the samplesheet the following process is still run to create the Salmon index in order to subsample the reads to infer the strandedness:

[3c/a91c08] process > NFCORE_RNASEQ:RNASEQ:FASTQ_SUBSAMPLE_FQ_SALMON:SALMON_INDEX (genome.transcripts.fa)             [100%] 1 of 1, cached: 1 ✔

This process shouldn't be run, however, it is being triggered since this logic evaluates to true:

!params.salmon_index && !('salmon' in prepareToolIndices)

We need to find a way to add to that logic somehow to also check that these channels are empty:

ch_strand_fastq.auto_strand,

If you have no intention of running pseudo-alignment, for now, the easiest workaround is to provide a dummy path that exists to --salmon_index /my/random/existing/folder/ as this won't actually be used or validated elsewhere by the pipeline and will save creating the index.

@drpatelh drpatelh added the bug Something isn't working label Mar 24, 2023
@drpatelh drpatelh added this to the 3.12 milestone Mar 24, 2023
drpatelh added a commit to drpatelh/nf-core-rnaseq that referenced this issue Mar 27, 2023
@drpatelh drpatelh modified the milestones: 3.12, 3.11 Mar 28, 2023
@drpatelh
Copy link
Member Author

Fixed in #978

drpatelh added a commit that referenced this issue Mar 28, 2023
Fix #975 and replace System.exit(1) calls in template
robsyme pushed a commit that referenced this issue Mar 29, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

1 participant