-
Notifications
You must be signed in to change notification settings - Fork 3
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
How to construct a complex indel from SNVs and simple indels #1
Comments
Thanks for your interest. This library expects you to write your own scripts to construct complex indels from simple indel inputs. For example, please see: Currently, this library constructs a complex indel from a simple indel, not from SNV. Please let me know if you need more info. Or, please contact me if you need a sample script with a sample dataset. |
Hi @rawagiha |
No, this is my fault!! Please apply phase(how="complex") to valn as fixed in the doc. My apologies about this. By the way, inputting SNV won't throw an error. Applying"phase" to VariantAlignment object representing a SNV returns a NullVariant object (https://indelpost.readthedocs.io/en/latest/api.html#indelpost.NullVariant). |
Additional note on constructing complex indels: In actual practice, we do not know in advance which SNVs and indels to combine. Or, the variant caller may not report all members of the target complex indel. Further, SNVs and Indels may be reported in two separate VCF files. For these reasons, indelPost only requires one indel to construct the target complex indel. It will search other members in the BAM file. In your case, the deletion at 55242469 is enough for construction. The other SNVs will be automatically searched. If the input indel is simple, the phase method returns the input indel as it is. |
Hi @rawagiha Thank you for your information. I got it working. I noticed a small typo in the doc at https://indelpost.readthedocs.io/en/latest/examples.html#annotating-complex-indels-from-simple-indels The line of code |
Wow! Thanks a lot for finding another typo (will fix it!). |
Hi, I read your paper, this repo, and the doc at https://indelpost.readthedocs.io/en/latest/api.html and I still do not know how to construct a complex indel from simple variants. For example, from a VCF file generated by Mutect/Strelka/UVC as input, how do I get a new VCF file with small variants fused into complex indels.
The text was updated successfully, but these errors were encountered: