fix rc assignment bug (multi-sample bambu-sc)

GoekeLab · Feb 17, 2025 · 4725911 · 4725911
1 parent 824db18
commit 4725911
Showing 1 changed file with 5 additions and 5 deletions.
diff --git a/R/bambu-processReads.R b/R/bambu-processReads.R
@@ -84,14 +84,14 @@ bambu.processReads <- function(reads, annotations, genomeSequence,
 
             mcols(readGrgList[[i]])$BC <- as.factor(mcols(readGrgList[[i]])$BC)
 
-            if(!isFALSE(demultiplexed)){ 
-                mcols(readGrgList[[i]])$sampleID <- as.numeric(mcols(readGrgList[[i]])$BC)
-            } else {
-                mcols(readGrgList[[i]])$sampleID <- i
-            }
         }
         readGrgList <- do.call(c, readGrgList)    
         mcols(readGrgList)$id <- seq_along(readGrgList) 
+        if(!isFALSE(demultiplexed)){ 
+          mcols(readGrgList[[i]])$sampleID <- as.numeric(mcols(readGrgList[[i]])$BC)
+        } else {
+          mcols(readGrgList[[i]])$sampleID <- i
+        }
         readClassList <- constructReadClasses(readGrgList, genomeSequence = genomeSequence,annotations = annotations,
             stranded = stranded, min.readCount = min.readCount, 
             fitReadClassModel = fitReadClassModel, min.exonOverlap = min.exonOverlap,