Mentions légales du service

Skip to content
Snippets Groups Projects
Commit 9845f58d authored by BOULLE Olivier's avatar BOULLE Olivier
Browse files

moved reverse complement to optimize

parent 68fd6549
No related branches found
No related tags found
No related merge requests found
...@@ -36,7 +36,6 @@ std::string get_cluster_id(const std::string& read, const std::string& seq_to_fi ...@@ -36,7 +36,6 @@ std::string get_cluster_id(const std::string& read, const std::string& seq_to_fi
*/ */
std::string read_rc = reverse_complement(read); // reverse complement of the read sequence
double score_fw, score_rc; // alignment score with the forward and reverse complement sequence double score_fw, score_rc; // alignment score with the forward and reverse complement sequence
int x_fw, y_fw, x_rc, y_rc; //x_fw and x_rc not used int x_fw, y_fw, x_rc, y_rc; //x_fw and x_rc not used
...@@ -49,6 +48,8 @@ std::string get_cluster_id(const std::string& read, const std::string& seq_to_fi ...@@ -49,6 +48,8 @@ std::string get_cluster_id(const std::string& read, const std::string& seq_to_fi
return read.substr(y_fw, split_level); return read.substr(y_fw, split_level);
} }
std::string read_rc = reverse_complement(read); // reverse complement of the read sequence
// test local alignement with the reverse complement sequence // test local alignement with the reverse complement sequence
SmithWaterman(read_rc, seq_to_find, score_rc, x_rc, y_rc); SmithWaterman(read_rc, seq_to_find, score_rc, x_rc, y_rc);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment