-
Notifications
You must be signed in to change notification settings - Fork 3
Description
Hi,
I used HMMSTR v 1.0.5 with coordinates on CANVAS disease involving RFC1 gene with differents motifs
chr4 39348424 39348479 AAAAG
chr4 39348424 39348479 AAGGG
chr4 39348424 39348479 AAAGG
The problem is that I got the exact same results for each motif which is quite impossible,
After looking at the result file (_genotype_calls.tsv) I saw that first column should be the ID but without filled the 5th column (Target name), every motif were named chr4:39348424-39348479,
Then I tried to set a uniq name for each in this column and now results are different for each CANVAS motif
You may either raise an error for multi-motif if target column is not filled or add motif as key to avoid this behaviour,
Best,
JB