diff --git a/this/src/fast_mapper.rs b/this/src/fast_mapper.rs index 1e8c6b2..5c46320 100644 --- a/this/src/fast_mapper.rs +++ b/this/src/fast_mapper.rs @@ -384,11 +384,9 @@ impl FastMapper{ let mut good = Vec::::with_capacity( genes.len() ); let most_matches = genes.values().max().unwrap_or(&0); let mut best_gene = 0; //otherwise this throws an error later - let mut genes_with_best_matches = 0; for ((gene_id, level ), matches) in genes{ if level == &prob_level && matches == most_matches { good.push( *gene_id ); - genes_with_best_matches +=1; } } if good.len() ==1 { @@ -398,13 +396,6 @@ impl FastMapper{ } return true } - if genes_with_best_matches == 1{ - ret.push( best_gene ); - if report { - eprintln!("3 This was selected as good: {} or {}\n",best_gene, self.names_store[best_gene] ) - } - return true - } // now we have one possibility left -> As many int reads as there are normal reads. // that is exactly what one would expect if it is actually a normal read