Skip SEQ when writing out

BimberLab · Feb 16, 2024 · 936c359 · 936c359
1 parent 959c8b6
commit 936c359
Showing 1 changed file with 10 additions and 2 deletions.
diff --git a/src/process/bam.rs b/src/process/bam.rs
@@ -19,13 +19,21 @@ use std::io::{BufWriter, Error, Write};
 const MAX_UMIS_IN_CHANNEL: usize = 50;
 
 fn bam_data_values(bam_data: &Vec<String>) -> String {
-    bam_data.join("\t")
+    bam_data
+        .iter()
+        .enumerate() 
+        .filter(|&(index, _)| index != 15) 
+        .map(|(_, value)| value.as_str())
+        .collect::<Vec<&str>>()
+        .join("\t")
 }
 
 fn bam_data_header(prefix: &str) -> String {
     BAM_FIELDS_TO_REPORT
         .iter()
-        .map(|&field| format!("{}_{}", prefix, field))
+        .enumerate() 
+        .filter(|&(index, _)| index != 15) 
+        .map(|(_, &field)| format!("{}_{}", prefix, field))
         .collect::<Vec<String>>()
         .join("\t")
 }