modules icon indicating copy to clipboard operation
modules copied to clipboard

SPADES changed flag for single-read in illumina reads

Open aidaanva opened this issue 2 years ago • 2 comments

PR checklist

Changed flag for single-read illumina reads to --s, since -s is used in spades for singletons from a pair-end library (see Input data in https://github.com/ablab/spades/blob/spades_3.15.5/README.md)

  • [x] This comment contains a description of changes (with reason).
  • [ ] If you've fixed a bug or added code that should be tested, add tests!
  • [ ] If you've added a new tool - have you followed the module conventions in the contribution docs
  • [ ] If necessary, include test data in your PR.
  • [ ] Remove all TODO statements.
  • [x] Emit the versions.yml file.
  • [ ] Follow the naming conventions.
  • [ ] Follow the parameters requirements.
  • [ ] Follow the input/output options guidelines.
  • [ ] Add a resource label
  • [x] Use BioConda and BioContainers if possible to fulfil software requirements.
  • Ensure that the test works with either Docker / Singularity. Conda CI tests can be quite flaky:
    • [x] PROFILE=docker pytest --tag <MODULE> --symlink --keep-workflow-wd --git-aware
    • [ ] PROFILE=singularity pytest --tag <MODULE> --symlink --keep-workflow-wd --git-aware
    • [ ] PROFILE=conda pytest --tag <MODULE> --symlink --keep-workflow-wd --git-aware

aidaanva avatar May 25 '23 12:05 aidaanva

Is this actually correct? --s comes under the 'multiple libraries' section, where my understanding of the implementation of this module is that it only ever accepts a single library (i.e. -1 vs -2, rather than --pe-1 -0pe-2 etc.?

The docs below would imply to me that for single-end libraries you would just give -1?

image

I tried to add the -1 flag as suggested by @jfy133, however when I used the -1 flag for the single-end alone the following error occurs:

== Error ==  the number of files with left paired reads is not equal to the number of
         files with right paired reads (library number: 1, library type: paired-end)! 

When I used the `--s1`` the input is interpreted correctly:

Dataset parameters:
  RNA virus assembly mode
  Metagenomic mode
  Reads:
    Library number: 1, library type: single
      left reads: not specified
      right reads: not specified
      interlaced reads: not specified
      single reads: ['/private/var/folders/1c/72tb52hs181_9xn2l2s9xt4m0000gp/T/tmp51itgo6g/63/1b394bce5e84e0aff90f28e1af303c/test_2.fastq.gz']
      merged reads: not specified

But I agree that the documentation seems to indicate otherwise

aidaanva avatar Jul 04 '23 07:07 aidaanva

@aidaanva , do you need some help to get this finished off?

SPPearce avatar May 13 '24 14:05 SPPearce

Hi @SPPearce, I just returned from maternity leave, I will try to get to this within the next month.

aidaanva avatar May 24 '24 09:05 aidaanva

@aidaanva We are doing some nf-core cleaning and trying to close PRs when not used anymore (they can be reopened in the future ofc). Just to know, are you planning to work on this one? :)

luisas avatar Mar 11 '25 12:03 luisas

I donor have the time to work on this one right now. I will reopen it when I have the time. Thank you for asking.

aidaanva avatar Mar 11 '25 14:03 aidaanva

@aidaanva All cool! Thanks for replying.

I will close it then :)

luisas avatar Mar 11 '25 14:03 luisas