Opened 9 days ago

Closed 8 days 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).

Change History (1)

comment:1 by Nicklas Nordborg, 8 days ago

Resolution: fixed
Status: newclosed

In 7868:

Fixes #1614: WGS variant calling fails when creating SBS96 counts

We only accept ACGT in the TRINUCLEO annotation.

Note: See TracTickets for help on using tickets.