|
|
@@ -303,6 +303,8 @@ impl SequencesGraph {
|
|
|
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] });
|
|
|
+ } else {
|
|
|
+
|
|
|
}
|
|
|
|
|
|
consensus_sequence.push(base);
|
|
|
@@ -311,6 +313,7 @@ impl SequencesGraph {
|
|
|
|
|
|
println!("RES {:?}", consensus_sequence);
|
|
|
println!("last = {:?}", consensus_sequence.last().unwrap().to_owned());
|
|
|
+ println!("last = {:?}", consensus_sequence.last().unwrap().to_owned() == b'N');
|
|
|
|
|
|
|
|
|
if consensus_sequence.last().unwrap().to_owned() == b'N' {
|
