Commit 9e2b476c authored by nchernia's avatar nchernia
Browse files

some systems do not have zcat. also grep now handles more complicated ligation patterns

parent 94e871e8
Loading
Loading
Loading
Loading
+3 −3
Changes for SLURM/scripts/countligations.sh: 3 added lines, 3 removed lines.
Original line number Diff line number Diff line
@@ -27,10 +27,10 @@
# Juicer version 1.5
if [ "$usegzip" -eq 1 ]
then 
    num1=$(paste <(zcat $name1$ext) <(zcat $name2$ext) | grep -c $ligation)
    num2=$(zcat ${name1}${ext} | wc -l | awk '{print $1}')
    num1=$(paste <(gunzip -c $name1$ext) <(gunzip -c $name2$ext) | grep -cE $ligation)
    num2=$(gunzip -c ${name1}${ext} | wc -l | awk '{print $1}')
else
    num1=$(paste $name1$ext $name2$ext | grep -c $ligation)
    num1=$(paste $name1$ext $name2$ext | grep -cE $ligation)
    num2=$(wc -l ${name1}${ext} | awk '{print $1}')
fi
echo -ne "$num1 " > ${name}${ext}_norm.txt.res.txt