!REQUIRES: python $VIDJIL_DIR/tools/check_python_version.py
!LAUNCH: $VIDJIL_DIR/vidjil $VIDJIL_DEFAULT_OPTIONS -3 -z 1 -G $VIDJIL_DIR/germline/IGH -w 60 -r 5 -e 10 -b data  $VIDJIL_DIR/data/Stanford_S22.fasta > /dev/null ; cat out/data.vidjil | python $VIDJIL_DIR/tools/format_json.py -1

$ Number of reads
e1:"total": [13153]

$ Number of segmented reads
e1:"segmented": [13153]

$ Most abundant window
1:"id": "CCACCTATTACTGTACCCGGGAGGAACAATATAGCAGCTGGTACTTTGACTTCTGGGGCC".*"reads": \\[8\\]

$ Affect values are over all the sequence
1: "affectValues": .[^}]*"start": 1, "stop": 128

$ Segmentation
1:"name": "IGHV3-23.05 6/ACCCGGGAGGAACAATAT/9 IGHD6-13.01 0//5 IGHJ4.02"

#                                                     Cys-T-=R=-E-=E=-Q-=Y=-S-=S=-W-=Y=-F-=D=-F-Trp
#          1         2         3         4         5  | *|  6         7   *     8    **   9    |  |10        11        12
# 12345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678
# 5-----------------------------------------------------5N----------------ND----------DJ-----------------------------------------J
# GCTGTACCTGCAAATGAACAGCCTGCGAGCCGAGGACACGGCCACCTATTACTGTACCCGGGAGGAACAATATAGCAGCTGGTACTTTGACTTCTGGGGCCAGGGGATCCTGGTCACCGTCTCCTCAG
#                                                        ACCCGGGAGGAACAATAT
#                                                                          agcagctggtac
# 5(V): 1 -> 55
# 4(D): 74 -> 85
# 3(J): 86 -> 128
# CDR3: 56 -> 94
# JUNCTION: 53 -> 97

$ Segmentation details - V
1:"5": ."delRight": 6, "name": "IGHV3-23.05", "stop": 55.

$ Segmentation details - D
1:"4": ."delLeft": 9, "delRight": 0, "name": "IGHD6-13.01", "start": 74, "stop": 85.

$ Segmentation details - J
1:"3": ."delLeft": 5, "name": "IGHJ4.02", "start": 86.

$ Segmentation details - N1, N2
1:"N1": 18,
1:"N2": 0,

$ Segmentation details - CDR3, JUNCTION
1:"cdr3": ."aa": "TREEQYSSWYFDF", "start": 56, "stop": 94.
1:"junction": ."aa": "CTREEQYSSWYFDFW", .* "start": 53, "stop": 97.

$ Second sequence has a DNA sequence provided
1:"id": "CTGTGCGAGAGGTTACTATGATAGTAGTGGTTATTACGGGGTAGGGCAGTACTACTACTA".*"sequence": "[ACGT]+",

$ Second sequence also has evalues
1:"id": "CTGTGCGAGAGGTTACTATGATAGTAGTGGTTATTACGGGGTAGGGCAGTACTACTACTA".*"evalue": ."val": "[0-9\.e-]+"

$ All 'start' fields are 1-based, they never equal to zero
0: "start": 0