do not index starbam

Author: kopardev

workflow/rules/rsem.smk

@@ -59,7 +59,7 @@ rule get_rsem_counts:
 set -exuf -o pipefail
 cd {params.rsemdir}
 # inter strandedness
-samtools index -@{threads} {input.starbam}
+# samtools index -@{threads} {input.starbam}
 infer_experiment.py -r {input.bed12} -i {input.starbam} -s 1000000 > {output.strandinfo}
 # Get strandedness to calculate Forward Probability
 fp=$(tail -n1 {output.strandinfo} | awk '{{if($NF > 0.75) print "0.0"; else if ($NF<0.25) print "1.0"; else print "0.5";}}')