Issue #1558: Fixed a bug with output of GX/GN BAM tags without CB/UB.…

… Bug-fix alpha patch: 2.7.10a_alpha_220517
alexdobin · May 17, 2022 · 071e989 · 071e989
1 parent 070b52e
commit 071e989
Show file tree

Hide file tree

Showing 3 changed files with 4 additions and 2 deletions.
diff --git a/CHANGES.md b/CHANGES.md
@@ -1,3 +1,4 @@
+* Issue #1558: Fixed a bug with output of GX/GN BAM tags without CB/UB.
 * Fixed a bug with --soloMultiMappers for small number of cells cases.
 * Issue #1513: If --soloMultiMappers options are not requested, output "NoMulti" in the "Reads Mapped to Gene: Unique+Multiple Gene" line of the Summary.csv file.
 * Fixed a problem with STARsolo CellReads.stats output for no-passlist runs.

diff --git a/source/ParametersSolo.cpp b/source/ParametersSolo.cpp
@@ -398,10 +398,11 @@ void ParametersSolo::initialize(Parameters *pPin)
     if (featureYes[SoloFeatureTypes::VelocytoSimple] || featureYes[SoloFeatureTypes::Velocyto]) {//turn readInfo on for Gene needed by VelocytoSimple
         readInfoYes[SoloFeatureTypes::Gene]=true;
     };
+    samAttrFeature = featureFirst;
     if (samAttrYes){//pSolo.samAttrFeature=0 by default, so need to check samAttrYes
         if (   featureFirst == SoloFeatureTypes::Gene || featureFirst == SoloFeatureTypes::GeneFull ||
                featureFirst == SoloFeatureTypes::GeneFull_Ex50pAS || featureFirst == SoloFeatureTypes::GeneFull_ExonOverIntron ) {
-            samAttrFeature = featureFirst;
+            //all good
         } else {
             ostringstream errOut;
             errOut << "EXITING because of fatal PARAMETERS error: CB and/or UB attributes in --outSAMattributes require --soloFeatures Gene OR/AND GeneFull OR/AND GeneFull_Ex50pAS.\n";

diff --git a/source/VERSION b/source/VERSION
@@ -1 +1 @@
-#define STAR_VERSION "2.7.10a_alpha_220506"
+#define STAR_VERSION "2.7.10a_alpha_220517"