correct bug in Assembler::correct_contigs()

zhixingfeng · Feb 14, 2020 · dadf534 · dadf534
1 parent 5456cf5
commit dadf534
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/...roj/project.xcworkspace/xcuserdata/zhixingfeng.xcuserdatad/UserInterfaceState.xcuserstate b/...roj/project.xcworkspace/xcuserdata/zhixingfeng.xcuserdatad/UserInterfaceState.xcuserstate
diff --git a/src/modules/assemble/assembler.cpp b/src/modules/assemble/assembler.cpp
@@ -1934,7 +1934,9 @@ void Assembler::correct_contigs(string ann_file, string out_file, double min_ove
                 is_corrected[i] = true;
                 for (auto k = 0; k < this->rl_ann_clust[j].cons_seq.size(); ++k){
                     if (cons_seq_hash[i].find(rl_ann_clust[j].cons_seq[k]) == cons_seq_hash[i].end() &&
-                        tested_loci_hash[i].find(rl_ann_clust[j].cons_seq[k]/4) == tested_loci_hash[i].end()){
+                        tested_loci_hash[i].find(rl_ann_clust[j].cons_seq[k]/4) == tested_loci_hash[i].end() &&
+                        int(rl_ann_clust[j].cons_seq[k]/4) >= rl_ann_clust[i].start &&
+                        int(rl_ann_clust[j].cons_seq[k]/4) <= rl_ann_clust[i].end){
 
                         cons_seq_corrected[i].insert(rl_ann_clust[j].cons_seq[k]);
                         tested_loci_corrected[i].insert(rl_ann_clust[j].cons_seq[k]/4);