Skip to content
Snippets Groups Projects
Commit 6d6ff73a authored by Gervaise Henry's avatar Gervaise Henry :cowboy:
Browse files

Fix log in infer species if ambiguous

parent 079c85ae
Branches
Tags
2 merge requests!58Develop,!56Resolve "Detect error in inferMetadata for tracking"
Pipeline #8858 canceled with stages
in 5 minutes and 44 seconds
......@@ -810,11 +810,13 @@ process inferMetadata {
species="Homo sapiens"
bam="GRCh.sampled.sorted.bam"
bed="./GRCh/genome.bed"
echo -e "LOG: inference of species results in: \${species}" >> ${repRID}.inferMetadata.log
elif [ 1 -eq \$(echo \$(expr \${align_mo} ">=" 40)) ] && [ 1 -eq \$(echo \$(expr \${align_hu} "<" 40)) ]
then
species="Mus musculus"
bam="GRCm.sampled.sorted.bam"
bed="./GRCm/genome.bed"
echo -e "LOG: inference of species results in: \${species}" >> ${repRID}.inferMetadata.log
else
echo -e "LOG: ERROR - inference of species returns an ambiguous result: hu=\${align_hu} mo=\${align_mo}" >> ${repRID}.inferMetadata.log
if [ "${speciesForce}" == "" ]
......@@ -836,8 +838,6 @@ process inferMetadata {
bam="GRCm.sampled.sorted.bam"
bed="./GRCm/genome.bed"
fi
else
echo -e "LOG: inference of species results in: \${species}" >> ${repRID}.inferMetadata.log
fi
# infer experimental setting from dedup bam
......
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