Skip to content
Snippets Groups Projects

Cutadapt set -m=1 to avoid empty reads

Merged CJHerrmann requested to merge cutadapt_min_len into dev

closes #161 (closed)
Cutadapt's default behavior is to keep all reads, even if they are empty after adapter removal. However, STAR will fail if encountering empty reads.
Solution:
Add -m=1 to all cutadapt rules, before the additional parameters from rule_config.yaml are inserted. Like this, if no rule_config.yaml is specified, -m=1 will be used, otherwise it will be overwritten by the value from rule_config.

The following tests have been performed locally:

  • confirm that cutadapt uses the second parameter value if a parameter is specified twice (tested with toy file in command line cutadapt call)
  • Specify the whole sequence of first read from input file as 3' adapter in the samples.tsv and run zarp with and without rule_config (This failed prior to adding -m=1)

Not tested:

  • cluster execution
  • whether the second value of -m is also used when cutadapt is called from snakemake (and not directly from command line)

Merge request reports

Loading
Loading

Activity

Filter activity
  • Approvals
  • Assignees & reviewers
  • Comments (from bots)
  • Comments (from users)
  • Commits & branches
  • Edits
  • Labels
  • Lock status
  • Mentions
  • Merge request status
  • Tracking
  • Loading
  • Loading
  • Loading
  • Loading
  • Loading
  • Loading
  • Loading
  • Loading
  • Loading
  • Loading
Please register or sign in to reply
Loading