Skip to content

Commit ea2f0ba

Browse files
committed
Additionally "per-read" collapse immunoseq dataset, probably fixes a reported bug
1 parent 32856dd commit ea2f0ba

File tree

2 files changed

+7
-2
lines changed

2 files changed

+7
-2
lines changed

src/main/groovy/com/antigenomics/vdjtools/misc/CommonUtil.groovy

Lines changed: 5 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -385,6 +385,11 @@ public class CommonUtil {
385385

386386
final static String OOF_SYMBOLS_POSSIBLE = /([atgc#~_\?])+/, OOF_CHAR = "_", STOP_CHAR = "*"
387387

388+
/**
389+
* Translates CDR3 sequence in V><J directions simultaneously, supports non-coding sequences
390+
* @param seq sequence to translate
391+
* @return amino acid sequence, as String object
392+
*/
388393
static String translate(String seq) {
389394
if (seq.length() == 0)
390395
return ""

src/main/groovy/com/antigenomics/vdjtools/misc/Software.groovy

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -49,13 +49,13 @@ enum Software {
4949
"cdr1nt", "cdr2nt", "cdr3nt", "cdr1aa", "cdr2aa", "cdr3aa",
5050
"inFrame", "noStop", "complete",
5151
4.times { _ }].flatten()),
52-
ImmunoSeq("\t", null, 1, true, false, ["cdr3nt", "cdr3aa", "count", "freq", "cdr3Length",
52+
ImmunoSeq("\t", null, 1, true, true, ["cdr3nt", "cdr3aa", "count", "freq", "cdr3Length",
5353
_, "v", 5.times { _ },
5454
_, "d", 5.times { _ },
5555
_, "j", 5.times { _ },
5656
6.times { _ },
5757
_, "VEnd", "DStart", "DEnd", "JStart"].flatten()),
58-
ImmunoSeqV3("\t", null, 1, true, false, ["cdr3nt", "cdr3aa", "count", "freq", "cdr3Length",
58+
ImmunoSeqV3("\t", null, 1, true, true, ["cdr3nt", "cdr3aa", "count", "freq", "cdr3Length",
5959
_, "v", 5.times { _ },
6060
_, "d", 5.times { _ },
6161
_, "j", 5.times { _ },

0 commit comments

Comments
 (0)