diff --git a/smwrappers/explore_mapping_parameters/wrapper.py b/smwrappers/explore_mapping_parameters/wrapper.py index 7f94718afbd403d149dc16c02413d7e45c1b3bc7..0cebcfcea98bf04bc43e5247ba69c81ce06f5a60 100644 --- a/smwrappers/explore_mapping_parameters/wrapper.py +++ b/smwrappers/explore_mapping_parameters/wrapper.py @@ -27,7 +27,7 @@ def mapping_command(aligner): cmd="niceload --noswap -q hisat2 -p {snakemake.threads} --dta --seed 123 -t {snakemake.params.these_settings} --mm -x {snakemake.params.index} -U {snakemake.params.fastq} --no-unal -S /dev/null" echo ${{cmd}} 1> {snakemake.log.log} eval ${{cmd}} 1>> {snakemake.log.log} 2> {snakemake.log.err} -rate=$(grep 'overall alignment rate' {snakemake.log.err} | tail -1 | tee -a {snakemake.log.log} | cut -d"%%" -f1) +rate=$(grep 'overall alignment rate' {snakemake.log.err} | tail -1 | tee -a {snakemake.log.log} | cut -d"%" -f1) echo -e "${{rate}}\\t{snakemake.params.these_settings}" >> {snakemake.params.rates_record} """ elif aligner == "bowtie2": @@ -36,7 +36,7 @@ echo -e "${{rate}}\\t{snakemake.params.these_settings}" >> {snakemake.params.rat cmd="niceload --noswap -q bowtie2 -p {snakemake.threads} --seed 123 -t {snakemake.params.these_settings} --mm -x {snakemake.params.index} -U {snakemake.params.fastq} --no-unal -S /dev/null" echo ${{cmd}} > {snakemake.log.log} eval ${{cmd}} 1>> {snakemake.log.log} 2> {snakemake.log.err} -rate=$(grep 'overall alignment rate' {snakemake.log.err} | tail -1 | tee -a {snakemake.log.log} | cut -d"%%" -f1) +rate=$(grep 'overall alignment rate' {snakemake.log.err} | tail -1 | tee -a {snakemake.log.log} | cut -d"%" -f1) echo -e "${{rate}}\\t{snakemake.params.these_settings}" >> {snakemake.params.rates_record} """ # elif aligner == "crac":