# Import the anarci function.
from anarci import anarci
# Format the sequences that we want to number.
sequences = [ ("12e8:H","AVVLTQTASPVSAAVGGTVTINCQASQSVYKNNRLAWYQQKPGQPPKLLIYGASTLASGVPSRFKGSGSGTQFTLTISDLECDDAATYYCAGGYSGNINVFGGGTEVVVK")]
# Hand the list of sequences to the anarci function. Number them with the IMGT scheme
results = anarci(sequences, scheme="imgt",output=False,assign_germline=True,allowed_species=['human',
'mouse','rat','rabbit','rhesus','pig','alpaca'])
# Unpack the results. We get three lists
numbering, alignment_details, hit_tables = results
print ("numbering: %s alignment_details: %s"%(numbering, alignment_details))
@broncio123 @FBoyles my test code is bellow:
the result of the alignment_details
alignment_details: [[{'id': 'rabbit_K', 'description': '', 'evalue': 3.1e-60, 'bitscore': 192.9, 'bias': 4.8, 'query_start': 0, 'query_end': 110, 'species': 'rabbit', 'chain_type': 'K', 'scheme': 'imgt', 'query_name': '12e8:H', 'germlines': {}}