public static void Test()
{
var bandOptions = new BandingOptions(4, 50);
var recursor = new SimpleQvRecursor((int)Move.ALL_MOVES, bandOptions);
var strandTpl = "ACGTACGTACGTACGT";
QvModelParams modelParams = QvModelParams.Default();
QvSequenceFeatures features = new QvSequenceFeatures("ACGTACGTCGT");
var evaluator = new QvEvaluator(features, strandTpl, modelParams);
var mutationEvaluator = new SimpleQvMutationScorer(evaluator, recursor);
var score = mutationEvaluator.Score();
Console.WriteLine(score);
}
public void Recursion()
{
var bandOptions = new BandingOptions(4, 50);
var recursor = new SimpleQvRecursor((int)(Move.BASIC_MOVES | Move.MERGE), bandOptions);
var strandTpl = "ACGTACGTACGTACGT";
QvModelParams modelParams = QvModelParams.Default();
QvSequenceFeatures features = new QvSequenceFeatures("ACGTACGTCGT");
var evaluator = new QvEvaluator(features, strandTpl, modelParams);
var mutationEvaluator = new SimpleQvMutationScorer(evaluator, recursor);
var score = mutationEvaluator.Score();
Console.WriteLine(score);
}
