1 #include "ebwt_search_util.h"
2 #include "seqan/file.h"
8 * Print a hit along with information about the backtracking
9 * regions constraining the hit.
11 void printHit(const vector<String<Dna5> >& os,
13 const String<Dna5>& qry,
21 // Print pattern sequence
22 cout << " Pat: " << qry << endl;
23 // Print text sequence
26 for(size_t i = 0; i < qlen; i++) {
27 cout << os[h.h.first][h.h.second + i];
30 for(int i = (int)qlen-1; i >= 0; i--) {
31 cout << os[h.h.first][h.h.second + i];
36 for(int i = (int)qlen-1; i >= 0; i--) {
37 if (i < (int)unrevOff) cout << "0";
38 else if(i < (int)oneRevOff) cout << "1";
39 else if(i < (int)twoRevOff) cout << "2";
40 else if(i < (int)threeRevOff) cout << "3";