diff --git a/workflow/rna-seq.nf b/workflow/rna-seq.nf index c7b9cb1ecaf6ce85a641248e706a5c3394336e56..a94077891768522aac7320ab35922c7ab3361bd4 100644 --- a/workflow/rna-seq.nf +++ b/workflow/rna-seq.nf @@ -599,12 +599,12 @@ process inferMetadata { echo -e "LOG: inference of strandedness results is: \${spike}" >> ${repRID}.inferMetadata.log # determine species - if [ 1 -eq \$(echo \$(expr \${align_hu} ">=" 25)) ] && [ 1 -eq \$(echo \$(expr \${align_mo} "<" 25)) ] + if [ 1 -eq \$(echo \$(expr \${align_hu} ">=" 40)) ] && [ 1 -eq \$(echo \$(expr \${align_mo} "<" 40)) ] then species="Homo sapiens" bam="GRCh.sampled.sorted.bam" bed="./GRCh/bed/genome.bed" - elif [ 1 -eq \$(echo \$(expr \${align_mo} ">=" 25)) ] && [ 1 -eq \$(echo \$(expr \${align_hu} "<" 25)) ] + elif [ 1 -eq \$(echo \$(expr \${align_mo} ">=" 40)) ] && [ 1 -eq \$(echo \$(expr \${align_hu} "<" 40)) ] then species="Mus musculus" bam="GRCm.sampled.sorted.bam"