diff --git a/dge_workflow/featcounts_deseq.R b/dge_workflow/featcounts_deseq.R
index d83df30c548f87376bc280a9321b7ea6087b6bb8..2a85e4a30d2403b52cdd321cb8a4b03ee30c93c4 100755
--- a/dge_workflow/featcounts_deseq.R
+++ b/dge_workflow/featcounts_deseq.R
@@ -73,7 +73,7 @@ genesAfter <- nrow(countMatrix)
 # TODO add general info about the deseq process/model here
 
 
-get_sample_from_replicate <- function(repName) str_match(repName, "(.*)_[0-9]{1}$")[,2]
+get_sample_from_replicate <- function(repName) str_match(repName, "(.*)_[R]?[0-9]{1}$")[,2]
 
 # Define or load a contrasts matrix
 if(!is.null(contrasts_file)){