Merge pull request #8 from klarman-cell-observatory/boli

Fixed a bug for detecting totalA or totalBC chemistry
lilab-bcb · May 17, 2021 · 0756c33 · 0756c33
2 parents dcc9a9c + cfca35b
commit 0756c33
Showing 1 changed file with 5 additions and 3 deletions.
diff --git a/generate_count_matrix_ADTs.cpp b/generate_count_matrix_ADTs.cpp
@@ -231,9 +231,11 @@ void detect_totalseq_type(string& extra_info) {
 			feature_iter = feature_index.find(binary_feature);
 			ntotA += (feature_iter != feature_index.end() && feature_iter->second.item_id >= 0);
 
-			binary_feature = barcode_to_binary(safe_substr(read2.seq, totalseq_BC_pos, feature_blen));
-			feature_iter = feature_index.find(binary_feature);
-			ntotBC += (feature_iter != feature_index.end() && feature_iter->second.item_id >= 0);
+			if (read2.seq.length() >= totalseq_BC_pos + feature_blen) {
+				binary_feature = barcode_to_binary(safe_substr(read2.seq, totalseq_BC_pos, feature_blen));
+				feature_iter = feature_index.find(binary_feature);
+				ntotBC += (feature_iter != feature_index.end() && feature_iter->second.item_id >= 0);				
+			}
 
 			++cnt;
 		}