GangSTR merge error

[Zhiqin]

Hi Helia

Thank you for developing EnsembleTR. When I merge vcfs from HipSTR, GangSTR and ExpansionHunter, I found that all but the vcf from B worked successfully. Then I use EnsembleTR to merge GangSTR vcf only, and I get the following error:

File "/data2/share/home/zhangzhiqin/.local/bin/EnsembleTR", line 33, in <module>
    sys.exit(load_entry_point('EnsembleTR==1.0.0', 'console_scripts', 'EnsembleTR')())
  File "/data2/share/home/zhangzhiqin/.local/lib/python3.9/site-packages/EnsembleTR-1.0.0-py3.9.egg/ensembletr/main.py", line 80, in run
  File "/data2/share/home/zhangzhiqin/.local/lib/python3.9/site-packages/EnsembleTR-1.0.0-py3.9.egg/ensembletr/main.py", line 41, in main
  File "/data2/share/home/zhangzhiqin/.local/lib/python3.9/site-packages/EnsembleTR-1.0.0-py3.9.egg/ensembletr/vcfio.py", line 228, in getMergableCalls
  File "/data2/share/home/zhangzhiqin/.local/lib/python3.9/site-packages/EnsembleTR-1.0.0-py3.9.egg/ensembletr/recordcluster.py", line 42, in __init__
  File "/data2/share/home/zhangzhiqin/.local/lib/python3.9/site-packages/trtools-5.0.2-py3.9.egg/trtools/utils/utils.py", line 363, in GetCanonicalMotif
    repseq_f = GetCanonicalOneStrand(repseq)
  File "/data2/share/home/zhangzhiqin/.local/lib/python3.9/site-packages/trtools-5.0.2-py3.9.egg/trtools/utils/utils.py", line 400, in GetCanonicalOneStrand
    if nucToNumber[newseq[j]] < nucToNumber[canonical[j]]:
KeyError: '1'

When running GangSTR, I ran all the samples together, and I don't know if that's the cause of the error.

Best,
Zhiqin