Merge pull request #428 from kopardev/activeDev

added missing tail -n1 bam2bw rule

Rules/initialqcrnaseq.snakefile

@@ -898,7 +898,7 @@ bamCoverage -b $bam -o {output.fbw} --filterRNAstrand forward --binSize 20 --smo
 bamCoverage -b $bam -o {output.rbw} --filterRNAstrand reverse --binSize 20 --smoothLength 40 -p 32
 
 # reverse files if method is not dUTP/NSR/NNSR ... ie, R1 in the direction of RNA strand.
-fp=`awk '{{if($NF > 0.75) print "0.0"; else if ($NF<0.25) print "1.0"; else print "0.5";}}' {input.strandinfo}`
+fp=`tail -n1 {input.strandinfo}|awk '{{if($NF > 0.75) print "0.0"; else if ($NF<0.25) print "1.0"; else print "0.5";}}'`
 if [ $fp -lt 0.25 ];then
 mv {output.fbw} {output.fbw}.tmp
 mv {output.rbw} {output.fbw}