Cutadapt set -m=1 to avoid empty reads
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.tsvand run zarp with and withoutrule_config(This failed prior to adding-m=1)
Not tested:
- cluster execution
- whether the second value of
-mis also used when cutadapt is called from snakemake (and not directly from command line)