Commit 0e595eb9 authored by Mathieu Giraud's avatar Mathieu Giraud

core/segment.{h,cpp}: compute and store CDR3productive

parent 80825945
......@@ -1076,6 +1076,8 @@ void FineSegmenter::findCDR3(){
JUNCTIONaa = nuc_to_aa(subsequence(getSequence().sequence, JUNCTIONstart, CDR3start-1))
+ CDR3aa + nuc_to_aa(subsequence(getSequence().sequence, CDR3end+1, JUNCTIONend));
CDR3productive = (CDR3nuc.length() % 3 == 0) && (CDR3aa.find('*') == string::npos);
}
json FineSegmenter::toJson(){
......
......@@ -126,6 +126,7 @@ protected:
int CDR3start, CDR3end;
string CDR3nuc;
string CDR3aa;
bool CDR3productive;
bool reversed, segmented, dSegmented;
int because;
......
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