diff --git a/src/RNAseq_cec.nf b/src/RNAseq_cec.nf
index 37f8cad09bb9a8f58fad9b8ab0d4bc42d7192a6b..2f19bcde2c6e002f729e5ac36e25a25ade8afc5a 100644
--- a/src/RNAseq_cec.nf
+++ b/src/RNAseq_cec.nf
@@ -26,6 +26,11 @@ include {
    star_mapping_fastq_out: "star_bam/"
 )
 
+include { index_bam
+  } from './nf_modules/samtools/main.nf' addParams(
+   htseq_out: "star_bam/"
+)
+
 include { htseq_count_with_gff
   } from './nf_modules/htseq/main.nf' addParams(
    htseq_out: "htseq_count/"
@@ -86,9 +91,12 @@ workflow {
       star_index.collect(),
       fastp.out.fastq,
   )
+
+  //// index bam
+  index_bam(star_mapping_fastq.out.bam)
   
   htseq_count_with_gff(
-        star_mapping_fastq.out.bam,
+        index_bam.out.bam_idx,
         genome_gff3_file
         )