Closed DarioS closed 4 years ago
Hi,
Thanks for the interest in our work and the feedback. I developed sigLASSO really as a new statistical learning algorithm to better assign mutational signatures. The algorithm inputs are the "mutational context" and "signatures".
To help users to use the tool, we supplied a simple basic bash script (get_context.sh) to preprocess the vcf/fasta and did an R wrapper on top of it. We encourage users to use this script as a prototype and modify it as needed.
I apologize that you found it hard to use. I have added a quick patch for this.
I am running
vcf2spec
but I see lots of errors like "WARNING. chromosome (1) was not found in the FASTA file. Skipping." The FASTA and VCF use the same style of chromosome names.The VCF file has the same style of chromosome IDs.
It looks like sigLASSO removes the chr prefix from the chromosome ID without asking the user, which is bad form. Surely, the user should be allowed to use chromosome names such as chr1, chr2.