Skip to content

Abnumber cannot renumber TCR with imgt #23

@ym-guan

Description

@ym-guan

Hi! I found that when i use abnumber to renumber TCR with imge scheme, it told me something error, like follow
`AssertionError Traceback (most recent call last)
Cell In[10], line 1
----> 1 abchain = abnumber.Chain(seq, 'imgt')

File ~/miniconda3/envs/diffab/lib/python3.8/site-packages/abnumber/chain.py:96, in Chain.init(self, sequence, scheme, cdr_definition, name, assign_germline, allowed_species, **kwargs)
94 if isinstance(sequence, Seq):
95 sequence = str(sequence)
---> 96 results = _anarci_align(sequence, scheme=scheme, allowed_species=allowed_species, assign_germline=assign_germline)
97 if len(results) > 1:
98 raise ChainParseError(f'Found {len(results)} antibody domains in sequence: "{sequence}"')

File ~/miniconda3/envs/diffab/lib/python3.8/site-packages/abnumber/common.py:43, in _anarci_align(sequence, scheme, allowed_species, assign_germline)
41 v_gene = ali['germlines']['v_gene'][0][1] if assign_germline else None
42 j_gene = ali['germlines']['j_gene'][0][1] if assign_germline else None
---> 43 aa_dict = {Position(chain_type=chain_type, number=num, letter=letter, scheme=scheme): aa
44 for (num, letter), aa in positions if aa != '-'}
45 tail = sequence[end+1:]
46 results.append((aa_dict, chain_type, tail, species, v_gene, j_gene))

File ~/miniconda3/envs/diffab/lib/python3.8/site-packages/abnumber/common.py:43, in (.0)
41 v_gene = ali['germlines']['v_gene'][0][1] if assign_germline else None
42 j_gene = ali['germlines']['j_gene'][0][1] if assign_germline else None
---> 43 aa_dict = {Position(chain_type=chain_type, number=num, letter=letter, scheme=scheme): aa
44 for (num, letter), aa in positions if aa != '-'}
...
18 def _validate_chain_type(chain_type):
---> 19 assert chain_type in ['H', 'L', 'K'],
20 f'Invalid chain type "{chain_type}", it should be "H" (heavy), "L" (lambda light chian) or "K" (kappa light chain)'

AssertionError: Invalid chain type "A", it should be "H" (heavy), "L" (lambda light chian) or "K" (kappa light chain)`

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions