Opened 7 weeks ago
Closed 7 weeks ago
#1614 closed defect (fixed)
WGS variant calling fails when creating SBS96 counts
Reported by: | Nicklas Nordborg | Owned by: | Nicklas Nordborg |
---|---|---|---|
Priority: | major | Milestone: | Reggie v5.3.1 |
Component: | net.sf.basedb.reggie | Keywords: | |
Cc: |
Description
The error:
java.lang.RuntimeException: Traceback (most recent call last): File "/home/scanb/lorry/jobs/WGSVariantCall-240930-1051-4345-4df359/sbs96.py", line 81, in <module> sbsCounts[sbs] = sbsCounts[sbs]+1 KeyError: 'C[C>T]N'
The cause seems to be a variant that is located where the neighbor position has a 'N' in the reference. The context becomes C[C>T]N
which is not one of the 96 SBS context. The script needs to be updated to handle this situation (eg. ignore the variant).
Note:
See TracTickets
for help on using tickets.
In 7868: