Commit 4e80e044 authored by Vidjil Team's avatar Vidjil Team

fuse.py: .clntab, CDR3 from JUNCTION

'sequence.JUNCTION.raw nt seq' is the JUNCTION, including V-Cys and J-Phe/J-Trp,
whereas the CDR3 does not include these amino acids.
parent dd18b7d4
......@@ -554,7 +554,7 @@ class ListWindows:
w.d["seg"]["5start"]=0
w.d["seg"]["4end"]=0
w.d["seg"]["4start"]=0
w.d["seg"]["cdr3"] = tab["sequence.JUNCTION.raw nt seq"]
w.d["seg"]["cdr3"] = tab["sequence.JUNCTION.raw nt seq"][3:-3]
listw.append((w , w.d["reads"][0]))
......
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