KeyError: 'V_SEQ_LENGTH'

Hello, I'm getting the following error when running MakeDb.py. Interestingly, this seems only to occur when running on T-cell receptor aligned output. Will continue digging...

(env2)dkoppstein@master:/dkoppstein/150521SG_v1.8_round2$ source /dkoppstein/150521SG_v1.8_round2/env/env2/bin/activate && python -m pdb `which MakeDb.py` igblast -i intermediate/150521SG_150520-1-P2_S2_L001_TR.fmt7 -s intermediate/150521SG_150520-1-P2_S2_L001_assembled_atleast-2.fasta --fail --outname 150521SG_150520-1-P2_S2_L001_TR --outdir intermediate
> /dkoppstein/src/bitbucket.org/javh/changeo/MakeDb.py(4)<module>()
-> """
(Pdb) c
        START> MakeDB
      ALIGNER> IgBlast
ALIGN_RESULTS> 150521SG_150520-1-P2_S2_L001_TR.fmt7
     SEQ_FILE> 150521SG_150520-1-P2_S2_L001_assembled_atleast-2.fasta
     NO_PARSE> False

PROGRESS> 00:16:37 [#                   ]   5% ( 3,440) 0.3 minTraceback (most recent call last):
  File "/usr/lib/python2.7/pdb.py", line 1314, in main
    pdb._runscript(mainpyfile)
  File "/usr/lib/python2.7/pdb.py", line 1233, in _runscript
    self.run(statement)
  File "/usr/lib/python2.7/bdb.py", line 400, in run
    exec cmd in globals, locals
  File "<string>", line 1, in <module>
  File "/dkoppstein/src/bitbucket.org/javh/changeo/MakeDb.py", line 4, in <module>
    """
  File "/dkoppstein/src/bitbucket.org/javh/changeo/MakeDb.py", line 516, in parseIgBlast
    no_parse=no_parse, out_args=out_args)
  File "/dkoppstein/src/bitbucket.org/javh/changeo/MakeDb.py", line 421, in writeDb
    for record in db_gen:
  File "/dkoppstein/src/bitbucket.org/javh/changeo/MakeDb.py", line 223, in readIgBlast
    db_gen['N1_LENGTH'] = int(d_align[8]) - int(db_gen['V_SEQ_LENGTH']) - int(db_gen['V_SEQ_START'])
KeyError: 'V_SEQ_LENGTH'
 Uncaught exception. Entering post mortem debugging
Running 'cont' or 'step' will restart the program
> /dkoppstein/src/bitbucket.org/javh/changeo/MakeDb.py(223)readIgBlast()
-> db_gen['N1_LENGTH'] = int(d_align[8]) - int(db_gen['V_SEQ_LENGTH']) - int(db_gen['V_SEQ_START'])
(Pdb) db_gen
{'JUNCTION': 'GGTCACCGTCTCCTCAGGG', 'IN_FRAME': 'F', 'D_CALL': 'TRBD1*01,TRBD2*02', 'J_CALL': 'None', 'STOP': 'T', 'FUNCTIONAL': 'F', 'SEQUENCE_INPUT': 'AGCTCTGGGAGAGGAGCCCAGGACTGGGATTCCGAGGTGTTTCCATTCAGTGATCTGCACTGAACACAGAGGACTCGCCATGGAGTTTGGGCTGAGCTGGGTTTTCCGTGTTGCTATTTTAAAAGGTGTCCAGTTTGAGGTGCGTCTGGTGGAGTCGGGGGGAGACTTGGTCCTGAGAGGGTCCCAGGGACGTTGCGGGGCGCCACCTCTTGCCTGGGGTCCTGGCACTGTTGTCACAATGTGACAGCTGGTTCGAGTCCTGGGGCCATGGAACCCTGGTCACCGTCTCCTCAGGGAGTGCATCCGCCC', 'SEQUENCE_ID': 'GTCATCCTTTGTTGTTT_GTTCGAGT|CONSCOUNT=2|PRCONS=0814-bs0-PCR2-IgM', 'D_SEQ_START': '292', 'V_CALL': 'None', 'JUNCTION_LENGTH': 19}
(Pdb) d_align
0                                                     D
1     GTCATCCTTTGTTGTTT_GTTCGAGT|CONSCOUNT=2|PRCONS=...
2                                              TRBD1*01
3                                                100.00
4                                                     5
5                                                     0
6                                                     0
7                                                     0
8                                                   292
9                                                   296
10                                                    5
11                                                    9
12                                                  8.4
13                                                 10.4
Name: 2, dtype: object
(Pdb) db_gen
{'JUNCTION': 'GGTCACCGTCTCCTCAGGG', 'IN_FRAME': 'F', 'D_CALL': 'TRBD1*01,TRBD2*02', 'J_CALL': 'None', 'STOP': 'T', 'FUNCTIONAL': 'F', 'SEQUENCE_INPUT': 'AGCTCTGGGAGAGGAGCCCAGGACTGGGATTCCGAGGTGTTTCCATTCAGTGATCTGCACTGAACACAGAGGACTCGCCATGGAGTTTGGGCTGAGCTGGGTTTTCCGTGTTGCTATTTTAAAAGGTGTCCAGTTTGAGGTGCGTCTGGTGGAGTCGGGGGGAGACTTGGTCCTGAGAGGGTCCCAGGGACGTTGCGGGGCGCCACCTCTTGCCTGGGGTCCTGGCACTGTTGTCACAATGTGACAGCTGGTTCGAGTCCTGGGGCCATGGAACCCTGGTCACCGTCTCCTCAGGGAGTGCATCCGCCC', 'SEQUENCE_ID': 'GTCATCCTTTGTTGTTT_GTTCGAGT|CONSCOUNT=2|PRCONS=0814-bs0-PCR2-IgM', 'D_SEQ_START': '292', 'V_CALL': 'None', 'JUNCTION_LENGTH': 19}

Corresponding entry in the .fmt7 file:

# IGBLASTN 2.2.29+
# Query: GTCATCCTTTGTTGTTT_GTTCGAGT|CONSCOUNT=2|PRCONS=0814-bs0-PCR2-IgM
# Database: db/TR.V.human.F+ORF+infrP.ungapped.fasta db/TR.D.human.F+ORF+infrP.ungapped.fasta db/TR.J.human.F+ORF+infrP.ungapped.fasta
# Domain classification requested: imgt

# V-(D)-J rearrangement summary for query sequence (Top V gene match, Top D gene match, Top J gene match, Chain type, stop codon, V-J frame, Productive, Strand).  Multiple equivalent top matches having the same score and percent identity, if present, are separated by a comma.
TRGVA*01        TRBD1*01,TRBD2*02       N/A     VB      Yes     N/A     No      +

# V-(D)-J junction details based on top germline gene matches (V end, V-D junction, D region, D-J junction, J start).  Note that possible overlapping nucleotides at VDJ junction (i.e, nucleotides that could be assigned to either rearranging gene) are indicated in parentheses (i.e., (TACT)) but are not included under the V, D, or J gene itself
GGTCA   CCGTCTCCT       CAGGG   N/A     N/A

# Hit table (the first field indicates the chain type of the hit)
# Fields: query id, subject id, % identity, alignment length, mismatches, gap opens, gaps, q. start, q. end, s. start, s. end, evalue, bit score
# 4 hits found
V       GTCATCCTTTGTTGTTT_GTTCGAGT|CONSCOUNT=2|PRCONS=0814-bs0-PCR2-IgM TRGVA*01        75.00   40      8       1
       2       245     282     16      55      0.17    26.8
V       GTCATCCTTTGTTGTTT_GTTCGAGT|CONSCOUNT=2|PRCONS=0814-bs0-PCR2-IgM TRBV6-8*01      65.22   46      16      0
       0       54      99      196     241     1.5     23.7
D       GTCATCCTTTGTTGTTT_GTTCGAGT|CONSCOUNT=2|PRCONS=0814-bs0-PCR2-IgM TRBD1*01        100.00  5       0       0
       0       292     296     5       9       8.4     10.4
D       GTCATCCTTTGTTGTTT_GTTCGAGT|CONSCOUNT=2|PRCONS=0814-bs0-PCR2-IgM TRBD2*02        100.00  5       0       0
       0       294     298     10      14      8.4     10.4

Comments (11)