Skip to content

Commit

Permalink
stop throwing exception on weird indel length thing
Browse files Browse the repository at this point in the history
since apparently it can still happen
  • Loading branch information
psathyrella committed Dec 9, 2020
1 parent b5436d7 commit f0c7f1a
Showing 1 changed file with 2 additions and 1 deletion.
3 changes: 2 additions & 1 deletion python/indelutils.py
Original file line number Diff line number Diff line change
Expand Up @@ -642,7 +642,8 @@ def combine_indels(regional_indelfos, full_qrseq, qrbounds, uid=None, debug=Fals
assert has_indels(rfo) # calling fcn needs to not add it if it doesn't have indels
joint_indelfo['genes'][region] = rfo['genes'][region]
if utils.non_gap_len(rfo['qr_gap_seq']) != qrbounds[region][1] - qrbounds[region][0]: # should be fixed by overlapping boundary shifter
raise Exception('%sqr_gap_seq non-gap length %d not the same as qrbound length %d in %s region indelfo' % ('%s: ' % uid if uid is not None else '', utils.non_gap_len(rfo['qr_gap_seq']), qrbounds[region][1] - qrbounds[region][0], region))
return None # UPDATE eh screw it this managed to happen *again* (see issue #310)
# raise Exception('%sqr_gap_seq non-gap length %d not the same as qrbound length %d in %s region indelfo' % ('%s: ' % uid if uid is not None else '', utils.non_gap_len(rfo['qr_gap_seq']), qrbounds[region][1] - qrbounds[region][0], region))
qr_gap_seq += [rfo['qr_gap_seq']]
gl_gap_seq += [rfo['gl_gap_seq']]

Expand Down

0 comments on commit f0c7f1a

Please sign in to comment.