Commit 4ac85f93 authored by Mikaël Salson's avatar Mikaël Salson

germline.cpp: The expected recombination is a TRDD2/TRDJ (see BIOMED-2 paper)

parent 1253a4f8
......@@ -201,7 +201,7 @@ void MultiGermline::build_incomplete_set(string path, int max_indexing)
add_germline(new Germline("TRA+D", 'a', path + "/TRDD_upstream.fa", "", path + "/TRAJ.fa", -10, 80, max_indexing), SEED_S13);
// DD-JD + DD2-DD3
add_germline(new Germline("TRD+", 'd', path + "/TRDD_upstream.fa", "", path + "/TRDJ.fa", -10, 60, max_indexing), SEED_9);
add_germline(new Germline("TRD+", 'd', path + "/TRDD2_upstream.fa", "", path + "/TRDJ.fa", -10, 60, max_indexing), SEED_9);
add_germline(new Germline("TRD+", 'd', path + "/TRDV.fa", "", path + "/TRDD3_downstream.fa", -10, 50, max_indexing), SEED_9);
add_germline(new Germline("TRD+", 'd', path + "/TRDD2_upstream.fa", "", path + "/TRDD3_downstream.fa", -10, 50, max_indexing), SEED_9);
......
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