diff --git a/dge_workflow/dge_utils.sh b/dge_workflow/dge_utils.sh
index 9cc7225539668586f3d6f39385d9d1f9eda86119..516d6095ec24a8a88f1fb103e1d718b23da3c4f5 100755
--- a/dge_workflow/dge_utils.sh
+++ b/dge_workflow/dge_utils.sh
@@ -19,6 +19,10 @@ export PATH=/home/brandl/bin/deepTools/bin:$PATH
 export PATH=/projects/bioinfo/holger/bin/FastQC_0.11.2:$PATH
 export PATH=/projects/bioinfo/holger/bin/bedtools-2.23.0/bin/:$PATH
 export PATH=/home/brandl/bin/spinr:$PATH
+export PATH=/home/brandl/bin/subread-1.4.6-p3-Linux-x86_64/bin:$PATH
+
+
+
 # which tophat;  which bowtie2; which cuffdiff
 
 
diff --git a/dge_workflow/star_align.sh b/dge_workflow/star_align.sh
index ce937ad8afdaefcbe31f2526e25804b06c3fec45..927db0354f7491bbbb3ecaf5e014be85fc5dbdd7 100755
--- a/dge_workflow/star_align.sh
+++ b/dge_workflow/star_align.sh
@@ -92,9 +92,18 @@ done
 
 wait4jobs .tophatjobs
 
+dge_bam_correlate . &
+
+
+## estimate expresssion with http://bioinf.wehi.edu.au/featureCounts/
+##Summarize multiple datasets at the same time:
+#featureCounts -t exon -g gene_id -a annotation.gtf -o counts.txt library1.bam library2.bam library3.bam
+
+#bamFile=control_3.bam
+#featureCounts -t exon -g gene_id -a ${gtfFile} -o feature_counts.txt ${bamFile} -T 5
+mysub "${project}__feature_counts" "featureCounts -t exon -g gene_id -a ${gtfFile} -o gene_counts.txt $(ls *bam) -T 5" | blockScript
 
 
-dge_bam_correlate . &
 
 ## create tophat mapping report
 ## todo adjust report to STAR