test

broadinstitute · Feb 4, 2024 · 2f10191 · 2f10191
1 parent a213288
commit 2f10191
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 4 deletions.
diff --git a/tools/TagSort/src/metricgatherer.cpp b/tools/TagSort/src/metricgatherer.cpp
@@ -202,7 +202,7 @@ CellMetricGatherer::CellMetricGatherer(std::string metric_output_file,
 {
   // write metrics csv header
   std::string s;
-  for (int i=0; i<25; i++)
+  for (int i=0; i<26; i++)
     metrics_csv_outfile_ << "," << kCommonHeaders[i]; // TODO ok to start with ,?
   for (int i=0; i<11; i++)
     metrics_csv_outfile_ << "," << cell_specific_headers[i];
@@ -324,7 +324,7 @@ GeneMetricGatherer::GeneMetricGatherer(std::string metric_output_file,
 {
   // write metrics csv header
   std::string s;
-  for (int i=0; i<25; i++)
+  for (int i=0; i<26; i++)
     metrics_csv_outfile_ << "," << kCommonHeaders[i]; // TODO ok to start with ,?
   for (int i=0; i<2; i++)
     metrics_csv_outfile_ << "," << gene_specific_headers[i];

diff --git a/tools/TagSort/src/metricgatherer.h b/tools/TagSort/src/metricgatherer.h
@@ -81,7 +81,7 @@ class MetricGatherer
   void clearCellAndGeneCommon();
   bool isMitochondrial(LineFields const& fields) const;
 
-  const std::string kCommonHeaders[25] =
+  const std::string kCommonHeaders[26] =
   {
     "n_reads",
     "noise_reads",
@@ -107,7 +107,8 @@ class MetricGatherer
     "reads_per_fragment",
     "fragments_per_molecule",
     "fragments_with_single_read_evidence",
-    "molecules_with_single_read_evidence"
+    "molecules_with_single_read_evidence",
+    "n_mitochondrial_reads"
   };
 
   void parseAlignedReadFields(LineFields const& fields, std::string hyphenated_tags);