Skip to content
Snippets Groups Projects
Commit e108c736 authored by Jonathan Gesell's avatar Jonathan Gesell
Browse files

Merge branch '44-stranded.logic' into '45-parallelizeTin'

Fix infer logic

See merge request !25
parents af133ded 28850994
Branches
Tags
4 merge requests!37v0.0.1,!27Develop,!26Resolve "Chunk bam for parallel tin calculation",!25Fix infer logic
Pipeline #6339 failed with stages
in 1 hour, 24 minutes, and 31 seconds
......@@ -431,7 +431,7 @@ process inferMetadata {
set path (inBam), path (inBai), path (inBamChr) from dedupBam_inferMeta
output:
path "infer.csv" into inferedMetadata
path "infer.tsv" into inferedMetadata
path "${inBam.baseName}.tin.xls" into tin
path "${repRID}.insertSize.inner_distance_freq.txt" optional true into innerDistance
......@@ -456,7 +456,7 @@ process inferMetadata {
percentF=`bash inferMeta.sh sef ${repRID}.rseqc.log`
percentR=`bash inferMeta.sh ser ${repRID}.rseqc.log`
fi
if [ \$percentF -gt 0.25 ] && [ \$percentR -lt 0.25 ]
if [ 1 -eq \$(echo \$(expr \$percentF ">" 0.25)) ] && [ 1 -eq \$(echo \$(expr \$percentR "<" 0.25)) ]
then
stranded="forward"
if [ \$endness == "PairEnd" ]
......@@ -465,7 +465,7 @@ process inferMetadata {
else
strategy="++,--"
fi
elif [ \$percentR -gt 0.25 ] && [ \$percentF -lt 0.25 ]
elif [ 1 -eq \$(echo \$(expr \$percentR ">" 0.25)) ] && [ 1 -eq \$(echo \$(expr \$percentF "<" 0.25)) ]
then
stranded="reverse"
if [ \$endness == "PairEnd" ]
......@@ -487,6 +487,6 @@ process inferMetadata {
# write infered metadata to file
echo \${endness},\${stranded},\${strategy},\${percentF},\${percentR},\${fail} > infer.csv
echo -e \${endness}'\\t'\${stranded}'\\t'\${strategy}'\\t'\${percentF}'\\t'\${percentR}'\\t'\${fail} > infer.tsv
"""
}
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment