SeqAnのチュートリアル

gistfile1.cpp

#include <iostream>
#include <string>

#include <seqan/sequence.h>

using namespace std;

int main () {
	seqan::String<char> text = "this is test using by seqan";
	seqan::String<char> pattern = "test";
	

	
	seqan::String<int> score;
	resize(score, seqan::length(text));
	
	for (unsigned i=0; i<seqan::length(text)-seqan::length(pattern)+1; i++) {
		int localScore = 0;
		for ( unsigned j = 0; j<seqan::length(pattern); j++) {
			if ( text[j+i] == pattern[j] ){
				localScore++;
			}
			score[i] = localScore;
		}
		
		for ( unsigned i = 0; i < seqan::length(score); i++){
			std::cout << score[i] << "hit count";
			std::cout << std::endl;
			
			return 0;
		}
	}
}