-
Notifications
You must be signed in to change notification settings - Fork 2
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
return no-effect when ter site is inserted around ter site and fix extension determine process #115
Conversation
Codecov ReportAttention: Patch coverage is
Additional details and impacted files@@ Coverage Diff @@
## master #115 +/- ##
==========================================
+ Coverage 46.93% 48.13% +1.19%
==========================================
Files 16 16
Lines 2124 2273 +149
Branches 71 72 +1
==========================================
+ Hits 997 1094 +97
- Misses 1056 1107 +51
- Partials 71 72 +1 ☔ View full report in Codecov by Sentry. |
(if new-aa-pos | ||
(coord/protein-coordinate new-aa-pos) | ||
(coord/unknown-coordinate))))))) | ||
:else |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
just change indent and move alt-prot-seq* let binding to line 581
(let [alt-prot-seq* (format-alt-prot-seq seq-info)] | ||
(cond | ||
(and (not= ppos 1) | ||
(ter-site-same-pos? ref-prot-seq alt-prot-seq*)) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
change c-ter-adjusted-alt-prot-seq
to alt-prot-seq*
t (cond | ||
ref-include-from-ter-start-and-over-ter-end :frame-shift | ||
(and ref-include-from-ter-start-and-over-ter-end | ||
(not first-diff-aa-is-ter-site)) :frame-shift |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
treat as extension when sequence is same until before ter site.
src/varity/vcf_to_hgvs/protein.clj
Outdated
(let [ref-ter-pos (count ref-prot-seq) | ||
alt-ter-pos (inc (count (first (string/split alt-prot-seq #"\*"))))] | ||
(and (string/includes? ref-prot-seq "*") | ||
(string/includes? alt-prot-seq "*") | ||
(= ref-ter-pos alt-ter-pos)))) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
There is a possibility of a NullPointerException occurring.
(let [ref-ter-pos (count ref-prot-seq) | |
alt-ter-pos (inc (count (first (string/split alt-prot-seq #"\*"))))] | |
(and (string/includes? ref-prot-seq "*") | |
(string/includes? alt-prot-seq "*") | |
(= ref-ter-pos alt-ter-pos)))) | |
(and (string/includes? ref-prot-seq "*") | |
(string/includes? alt-prot-seq "*") | |
(let [ref-ter-pos (count ref-prot-seq) | |
alt-ter-pos (inc (count (first (string/split alt-prot-seq #"\*"))))] | |
(= ref-ter-pos alt-ter-pos)))) |
083360e
to
aaf40df
Compare
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Thank you for promptly resolving the issue! LGTM 🙆♂️
I found
:clj-hgvs.coordinate/protein-coordinate
error when ter site is inserted around ter site and some extension varints are determined as frameshift.So I fixed protein-alteration type determine process.