Commit 826a4bbd authored by Mathieu Giraud's avatar Mathieu Giraud

core/segment.cpp: toJson, remove spurious output of D gene

score_D may be not properly initialized, leading to spurious output
of D genes in the .vidjil file (as for TRB-VJ recombinations in bug20150909.fa)

The correct way to check for D is to use isDSegmented.
parent b502ac3a
......@@ -999,7 +999,7 @@ json FineSegmenter::toJson(){
seg["5start"] = 0;
seg["5end"] = Vend;
if (score_D.size()>0){
if (isDSegmented()) {
seg["4"] = segmented_germline->rep_4.label(best_D);
seg["4start"] = Dstart;
seg["4end"] = Dend;
......
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