Fix some gene names Seurat mangles

code/run-seurat.R

@@ -61,6 +61,10 @@ run_seurat <- function(sce, pars) {
 
   result <- tibble::as_tibble(result)
 
+  # Fix some mangled gene names.
+  result$gene <- ifelse(result$gene == "Hba-a2-loc1", "Hba-a2_loc1", result$gene)
+  result$gene <- ifelse(result$gene == "Hba-a2-loc2", "Hba-a2_loc2", result$gene)
+
   validate_mgs_result(result)
 
   # Get the elapsed time.