Commit 3b37a335 authored by Mathieu Giraud's avatar Mathieu Giraud
Browse files

germline.cpp: seed for incomplete DD2-DD3 rearrangements has to be short

The sequence TRDD2-01 is only 9 bp.
parent d2a0dbd5
......@@ -179,13 +179,13 @@ void MultiGermline::build_incomplete_set(string path)
// DD2-DD3
germline = new Germline("TRD+", 'd', path + "/TRDD2-01.fa", "", path + "/TRDJ.fa", -10, 20);
germline->new_index("#####-#####");
germline->new_index("#########");
germlines.push_back(germline);
germline = new Germline("TRD+", 'd', path + "/TRDV.fa", "", path + "/TRDD3-01.fa", -10, 20);
germline->new_index("#####-#####");
germline->new_index("#########");
germlines.push_back(germline);
germline = new Germline("TRD+", 'd', path + "/TRDD2-01.fa", "", path + "/TRDD3-01.fa", -10, 20);
germline->new_index("#####-#####");
germline->new_index("#########");
germlines.push_back(germline);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment