Skip to content
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

Fix vcf-to-hgvs protein to handle initiation or stop codon variant more precisely #97

Open
nokara26 opened this issue Apr 12, 2024 · 0 comments

Comments

@nokara26
Copy link
Contributor

I encountered a failure while converting with vcf-to-hgvs for the variant below.

These are stop codon deletion variant.

We couldn't handle precisely these variants.
So we need to add conceivable cases which affect initiation or stop codon and implement more precisely.

error case 1

{:chr "chr3" :pos 41239336 :ref "CCTGTAAATCATCCTT" :alt "C"}

message

Cannot invoke "Object.getClass()" because "x" is null

...
varity.vcf-to-hgvs.protein/alt-exon-ranges (protein.clj:60)
...
varity.vcf-to-hgvs.protein/alt-exon-ranges (protein.clj:49)
varity.vcf-to-hgvs.protein/alt-exon-ranges (protein.clj:35)
varity.vcf-to-hgvs.protein/apply-offset (protein.clj:183)
varity.vcf-to-hgvs.protein/apply-offset (protein.clj:185)
varity.vcf-to-hgvs.protein/apply-offset (protein.clj:189)
varity.vcf-to-hgvs.protein/apply-offset (protein.clj:180)

error case 2

{:chr "chr3" :pos 41239336 :ref "CCTGTAAATC" :alt "C"}

message

Range [779, 713) out of bounds for length 783

...
varity.vcf-to-hgvs.protein/->protein-variant (protein.clj:306)
varity.vcf-to-hgvs.protein/->protein-variant (protein.clj:278)
varity.vcf-to-hgvs.protein/mutation (protein.clj:510)
varity.vcf-to-hgvs.protein/mutation (protein.clj:507)
varity.vcf-to-hgvs.protein/->hgvs (protein.clj:532)
varity.vcf-to-hgvs.protein/->hgvs (protein.clj:528)
...
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant