Open SHuang-Broad opened 6 years ago
Sure will upload, was actually delayed because ultimately I would like to have this functionality integrated into the main of Mobster. Regarding the first issue, this may be difficult, need to take into account TSD and aswell as non-preicise predictions. Anyway, first things, first will upload the code. Thanks for contributing!
Thanks, Jayne!
If I understand correctly, the REF allele problem is not too difficult to solve, because per the VCF spec, the REF allele must be at least one base long, and start with the base pointed to by POS
.
And since we are dealing with insertions here, the REF allele should always be 1 base long, right?
Yes this is true, was thinking about TSD's and events with no supporting split reads. But most I think you idea is the most practical solution. (Have you seen the new code in devlop?)
Yes, I just saw that after you mentioned. Let me know if there's anything I can help with that (test run, code review etc.)
Hi,
No REF allele issue was fixed or not? I could not visualize my mobster vcf files in IGV.
Thanks.
I've found (actually IGV did) two issues with the output from this tool:
REF
allele is outputINFO
annotation value, "no tsd", that contains a space. VCF spec does not allow space inINFO
.The second problem is easier to fix, while the first one needs more work (it needs access to reference fasta at the least).
It would be nice if the source code for that could be made public and I am happy to issue PRs for these two.
Thanks!