|
|
@@ -297,17 +297,23 @@ impl SequencesGraph {
|
|
|
}
|
|
|
}
|
|
|
});
|
|
|
+
|
|
|
let n_base: usize = n_atcg.iter().sum();
|
|
|
let max_base = n_atcg.iter().max().unwrap();
|
|
|
let mut base = b'N';
|
|
|
if *max_base as f32 / n_base as f32 > 0.5 {
|
|
|
n_atcg.iter().enumerate().for_each(|(pos, n)| if n == max_base { base = atcg[pos] });
|
|
|
}
|
|
|
+
|
|
|
consensus_sequence.push(base);
|
|
|
}
|
|
|
}
|
|
|
|
|
|
println!("RES {:?}", consensus_sequence);
|
|
|
+
|
|
|
+ if consensus_sequence[consensus_sequence.len() - 1] == b'N' {
|
|
|
+ println!("NNNNNNNNNNNNNNNNNNNnn!")
|
|
|
+ }
|
|
|
consensus_sequence
|
|
|
}
|
|
|
|
