VEP output SIFT_score unclear Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern) Announcing the arrival of Valued Associate #679: Cesar Manara Unicorn Meta Zoo #1: Why another podcast?Why are there missing calls in a VCF file from exome sequencing?Selecting sites from VCF which have an alt AD > 10Keep Format and Individual fields when annotating VCF with VEPupdate dbSNP IDWhy Ti/Tv ratio?Meaning of the FORMAT fields of the VCF file coming from GIAB projectAnnotation with Prokka or RAST.Efficiently aligning a lot of reads on the same small reference sequenceFastqc- Per Base Sequence QualityHow to correctly call a VCF file using damaged DNA? (IonTorrent & FFPE)
What are the motives behind Cersei's orders given to Bronn?
Are variable time comparisons always a security risk in cryptography code?
Stars Make Stars
Does accepting a pardon have any bearing on trying that person for the same crime in a sovereign jurisdiction?
Do you forfeit tax refunds/credits if you aren't required to and don't file by April 15?
Bonus calculation: Am I making a mountain out of a molehill?
When -s is used with third person singular. What's its use in this context?
How was the dust limit of 546 satoshis was chosen? Why not 550 satoshis?
Is there a concise way to say "all of the X, one of each"?
What's the purpose of writing one's academic bio in 3rd person?
What would be the ideal power source for a cybernetic eye?
Area of a 2D convex hull
Is there a documented rationale why the House Ways and Means chairman can demand tax info?
What does '1 unit of lemon juice' mean in a grandma's drink recipe?
How much radiation do nuclear physics experiments expose researchers to nowadays?
How can whole tone melodies sound more interesting?
How can players work together to take actions that are otherwise impossible?
Did Xerox really develop the first LAN?
What is this single-engine low-wing propeller plane?
How to motivate offshore teams and trust them to deliver?
If a contract sometimes uses the wrong name, is it still valid?
Java 8 stream max() function argument type Comparator vs Comparable
Why was the term "discrete" used in discrete logarithm?
How discoverable are IPv6 addresses and AAAA names by potential attackers?
VEP output SIFT_score unclear
Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern)
Announcing the arrival of Valued Associate #679: Cesar Manara
Unicorn Meta Zoo #1: Why another podcast?Why are there missing calls in a VCF file from exome sequencing?Selecting sites from VCF which have an alt AD > 10Keep Format and Individual fields when annotating VCF with VEPupdate dbSNP IDWhy Ti/Tv ratio?Meaning of the FORMAT fields of the VCF file coming from GIAB projectAnnotation with Prokka or RAST.Efficiently aligning a lot of reads on the same small reference sequenceFastqc- Per Base Sequence QualityHow to correctly call a VCF file using damaged DNA? (IonTorrent & FFPE)
$begingroup$
We have been experimenting with VEP (Variant Effect Predictor). One of the meta data attributes that we are interested in is the SIFT score, indeed when we apply the dbNSFP plug we get a column containing the scores (named SIFT_score). However, I don't understand why there are sometimes dots or multiple values in the fields. For example, the gene ENSG00000196924 below has 5 transcripts:
The SIFT_score column contains several values, not 1 per transcript/rs-number...
Here is another example that confuses me (I added the SIFT_pred column this time):
There are two mutations, the lower one can be expressed in 4 transcripts and I thus understand that there can be 4 SIFT scores, but why are all for in given in every row? Is the first one the SIFT_score for the first transcript?
One last example, again 4 transcripts, but now 2 of the scores are dots, what does that mean?
I have been looking for quite some time now how to interpret this data, any help is appreciated.
ngs variant-calling vep variant-effect-predictor
asked Apr 10 at 8:17
FreekFreek
2176
$endgroup$
add a comment |
$begingroup$
We have been experimenting with VEP (Variant Effect Predictor). One of the meta data attributes that we are interested in is the SIFT score, indeed when we apply the dbNSFP plug we get a column containing the scores (named SIFT_score). However, I don't understand why there are sometimes dots or multiple values in the fields. For example, the gene ENSG00000196924 below has 5 transcripts:
The SIFT_score column contains several values, not 1 per transcript/rs-number...
Here is another example that confuses me (I added the SIFT_pred column this time):
There are two mutations, the lower one can be expressed in 4 transcripts and I thus understand that there can be 4 SIFT scores, but why are all for in given in every row? Is the first one the SIFT_score for the first transcript?
One last example, again 4 transcripts, but now 2 of the scores are dots, what does that mean?
I have been looking for quite some time now how to interpret this data, any help is appreciated.
ngs variant-calling vep variant-effect-predictor
asked Apr 10 at 8:17
FreekFreek
2176
$endgroup$
add a comment |
$begingroup$
We have been experimenting with VEP (Variant Effect Predictor). One of the meta data attributes that we are interested in is the SIFT score, indeed when we apply the dbNSFP plug we get a column containing the scores (named SIFT_score). However, I don't understand why there are sometimes dots or multiple values in the fields. For example, the gene ENSG00000196924 below has 5 transcripts:
The SIFT_score column contains several values, not 1 per transcript/rs-number...
Here is another example that confuses me (I added the SIFT_pred column this time):
There are two mutations, the lower one can be expressed in 4 transcripts and I thus understand that there can be 4 SIFT scores, but why are all for in given in every row? Is the first one the SIFT_score for the first transcript?
One last example, again 4 transcripts, but now 2 of the scores are dots, what does that mean?
I have been looking for quite some time now how to interpret this data, any help is appreciated.
ngs variant-calling vep variant-effect-predictor
asked Apr 10 at 8:17
FreekFreek
2176
$endgroup$
We have been experimenting with VEP (Variant Effect Predictor). One of the meta data attributes that we are interested in is the SIFT score, indeed when we apply the dbNSFP plug we get a column containing the scores (named SIFT_score). However, I don't understand why there are sometimes dots or multiple values in the fields. For example, the gene ENSG00000196924 below has 5 transcripts:
The SIFT_score column contains several values, not 1 per transcript/rs-number...
Here is another example that confuses me (I added the SIFT_pred column this time):
There are two mutations, the lower one can be expressed in 4 transcripts and I thus understand that there can be 4 SIFT scores, but why are all for in given in every row? Is the first one the SIFT_score for the first transcript?
One last example, again 4 transcripts, but now 2 of the scores are dots, what does that mean?
I have been looking for quite some time now how to interpret this data, any help is appreciated.
ngs variant-calling vep variant-effect-predictor
ngs variant-calling vep variant-effect-predictor
asked Apr 10 at 8:17
FreekFreek
2176
asked Apr 10 at 8:17
FreekFreek
2176
asked Apr 10 at 8:17
FreekFreek
2176
asked Apr 10 at 8:17
FreekFreek
2176
asked Apr 10 at 8:17
FreekFreek
2176
2176
add a comment |
add a comment |
2 Answers
2
active
oldest
votes
$begingroup$
The dbNSFP plugin from VEP accesses tables of data for each variant from dbNSFP and pulls out the values. dbNSFP provide their SIFT scores in that format: a score for every transcript affected by the variant, all on one line. The lookup is just for the variant, not for the variant/transcript combo, so they provide scores for all variant/transcript combos. You can also get a column that gives you a list of the transcripts or proteins (Ensembl_transcriptid or Ensembl_proteinid) in order so you know which score goes with which transcript.
A better way to get SIFT scores with VEP is to get them directly from VEP, rather than using dbNSFP. This will get the SIFT score that goes with the transcript on the line with the relevant transcript.
answered Apr 10 at 10:05
Emily_EnsemblEmily_Ensembl
1,06918
$endgroup$
$begingroup$
I am guessing the dots are there for cases where dbNSFP doesn't have a value for the relevant transcript, right?
$endgroup$
– terdon♦
Apr 10 at 12:18
$begingroup$
Yes, that's it. Could be that the variant isn't missense in that transcript.
$endgroup$
– Emily_Ensembl
Apr 10 at 13:00
add a comment |
$begingroup$
The first gene you mention, ENSG00000196924, actually has 6 transcripts (link to the VarSome.com page of variant rs371839875), not 5. It's just that one of them is non-coding:
So the Sift scores you see are indeed one per transcript, it's just that there are 6 because dbNSFP also includes a score for the non-coding transcript of the gene.
The dots are just there as placeholders, they mean there was no value associated with that transcript. Many tools will show some sort of symbol instead of an empty field both for clarity and for practical technical reasons.
Visiting the variant's page on VarSome gives you a clearer picture since we collapse the identical scores and also include the converted rankscore provided by dbNSFP so you can have a single number for your variant:
Disclaimer: I work for the company behind VarSome, but it's a free tool. You need to pay to annotate VCF files (unlike the 100% free VEP), but it's free to use as a lookup tool for single variants.
answered Apr 10 at 12:06
terdon♦terdon
4,7752830
$endgroup$
add a comment |
Your Answer
StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "676"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);
else
createEditor();
);
function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);
);
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fbioinformatics.stackexchange.com%2fquestions%2f7430%2fvep-output-sift-score-unclear%23new-answer', 'question_page');
);
Post as a guest
Required, but never shown
2 Answers
2
active
oldest
votes
2 Answers
2
active
oldest
votes
active
oldest
votes
active
oldest
votes
$begingroup$
The dbNSFP plugin from VEP accesses tables of data for each variant from dbNSFP and pulls out the values. dbNSFP provide their SIFT scores in that format: a score for every transcript affected by the variant, all on one line. The lookup is just for the variant, not for the variant/transcript combo, so they provide scores for all variant/transcript combos. You can also get a column that gives you a list of the transcripts or proteins (Ensembl_transcriptid or Ensembl_proteinid) in order so you know which score goes with which transcript.
A better way to get SIFT scores with VEP is to get them directly from VEP, rather than using dbNSFP. This will get the SIFT score that goes with the transcript on the line with the relevant transcript.
answered Apr 10 at 10:05
Emily_EnsemblEmily_Ensembl
1,06918
$endgroup$
$begingroup$
I am guessing the dots are there for cases where dbNSFP doesn't have a value for the relevant transcript, right?
$endgroup$
– terdon♦
Apr 10 at 12:18
$begingroup$
Yes, that's it. Could be that the variant isn't missense in that transcript.
$endgroup$
– Emily_Ensembl
Apr 10 at 13:00
add a comment |
$begingroup$
The dbNSFP plugin from VEP accesses tables of data for each variant from dbNSFP and pulls out the values. dbNSFP provide their SIFT scores in that format: a score for every transcript affected by the variant, all on one line. The lookup is just for the variant, not for the variant/transcript combo, so they provide scores for all variant/transcript combos. You can also get a column that gives you a list of the transcripts or proteins (Ensembl_transcriptid or Ensembl_proteinid) in order so you know which score goes with which transcript.
A better way to get SIFT scores with VEP is to get them directly from VEP, rather than using dbNSFP. This will get the SIFT score that goes with the transcript on the line with the relevant transcript.
answered Apr 10 at 10:05
Emily_EnsemblEmily_Ensembl
1,06918
$endgroup$
$begingroup$
I am guessing the dots are there for cases where dbNSFP doesn't have a value for the relevant transcript, right?
$endgroup$
– terdon♦
Apr 10 at 12:18
$begingroup$
Yes, that's it. Could be that the variant isn't missense in that transcript.
$endgroup$
– Emily_Ensembl
Apr 10 at 13:00
add a comment |
$begingroup$
The dbNSFP plugin from VEP accesses tables of data for each variant from dbNSFP and pulls out the values. dbNSFP provide their SIFT scores in that format: a score for every transcript affected by the variant, all on one line. The lookup is just for the variant, not for the variant/transcript combo, so they provide scores for all variant/transcript combos. You can also get a column that gives you a list of the transcripts or proteins (Ensembl_transcriptid or Ensembl_proteinid) in order so you know which score goes with which transcript.
A better way to get SIFT scores with VEP is to get them directly from VEP, rather than using dbNSFP. This will get the SIFT score that goes with the transcript on the line with the relevant transcript.
answered Apr 10 at 10:05
Emily_EnsemblEmily_Ensembl
1,06918
$endgroup$
The dbNSFP plugin from VEP accesses tables of data for each variant from dbNSFP and pulls out the values. dbNSFP provide their SIFT scores in that format: a score for every transcript affected by the variant, all on one line. The lookup is just for the variant, not for the variant/transcript combo, so they provide scores for all variant/transcript combos. You can also get a column that gives you a list of the transcripts or proteins (Ensembl_transcriptid or Ensembl_proteinid) in order so you know which score goes with which transcript.
A better way to get SIFT scores with VEP is to get them directly from VEP, rather than using dbNSFP. This will get the SIFT score that goes with the transcript on the line with the relevant transcript.
answered Apr 10 at 10:05
Emily_EnsemblEmily_Ensembl
1,06918
edited Apr 10 at 10:45
answered Apr 10 at 10:05
Emily_EnsemblEmily_Ensembl
1,06918
answered Apr 10 at 10:05
Emily_EnsemblEmily_Ensembl
1,06918
answered Apr 10 at 10:05
Emily_EnsemblEmily_Ensembl
1,06918
1,06918
$begingroup$
I am guessing the dots are there for cases where dbNSFP doesn't have a value for the relevant transcript, right?
$endgroup$
– terdon♦
Apr 10 at 12:18
$begingroup$
Yes, that's it. Could be that the variant isn't missense in that transcript.
$endgroup$
– Emily_Ensembl
Apr 10 at 13:00
add a comment |
$begingroup$
I am guessing the dots are there for cases where dbNSFP doesn't have a value for the relevant transcript, right?
$endgroup$
– terdon♦
Apr 10 at 12:18
$begingroup$
Yes, that's it. Could be that the variant isn't missense in that transcript.
$endgroup$
– Emily_Ensembl
Apr 10 at 13:00
$begingroup$
I am guessing the dots are there for cases where dbNSFP doesn't have a value for the relevant transcript, right?
$endgroup$
– terdon♦
Apr 10 at 12:18
$begingroup$
I am guessing the dots are there for cases where dbNSFP doesn't have a value for the relevant transcript, right?
$endgroup$
– terdon♦
Apr 10 at 12:18
$begingroup$
Yes, that's it. Could be that the variant isn't missense in that transcript.
$endgroup$
– Emily_Ensembl
Apr 10 at 13:00
$begingroup$
Yes, that's it. Could be that the variant isn't missense in that transcript.
$endgroup$
– Emily_Ensembl
Apr 10 at 13:00
add a comment |
$begingroup$
The first gene you mention, ENSG00000196924, actually has 6 transcripts (link to the VarSome.com page of variant rs371839875), not 5. It's just that one of them is non-coding:
So the Sift scores you see are indeed one per transcript, it's just that there are 6 because dbNSFP also includes a score for the non-coding transcript of the gene.
The dots are just there as placeholders, they mean there was no value associated with that transcript. Many tools will show some sort of symbol instead of an empty field both for clarity and for practical technical reasons.
Visiting the variant's page on VarSome gives you a clearer picture since we collapse the identical scores and also include the converted rankscore provided by dbNSFP so you can have a single number for your variant:
Disclaimer: I work for the company behind VarSome, but it's a free tool. You need to pay to annotate VCF files (unlike the 100% free VEP), but it's free to use as a lookup tool for single variants.
answered Apr 10 at 12:06
terdon♦terdon
4,7752830
$endgroup$
add a comment |
$begingroup$
The first gene you mention, ENSG00000196924, actually has 6 transcripts (link to the VarSome.com page of variant rs371839875), not 5. It's just that one of them is non-coding:
So the Sift scores you see are indeed one per transcript, it's just that there are 6 because dbNSFP also includes a score for the non-coding transcript of the gene.
The dots are just there as placeholders, they mean there was no value associated with that transcript. Many tools will show some sort of symbol instead of an empty field both for clarity and for practical technical reasons.
Visiting the variant's page on VarSome gives you a clearer picture since we collapse the identical scores and also include the converted rankscore provided by dbNSFP so you can have a single number for your variant:
Disclaimer: I work for the company behind VarSome, but it's a free tool. You need to pay to annotate VCF files (unlike the 100% free VEP), but it's free to use as a lookup tool for single variants.
answered Apr 10 at 12:06
terdon♦terdon
4,7752830
$endgroup$
add a comment |
$begingroup$
The first gene you mention, ENSG00000196924, actually has 6 transcripts (link to the VarSome.com page of variant rs371839875), not 5. It's just that one of them is non-coding:
So the Sift scores you see are indeed one per transcript, it's just that there are 6 because dbNSFP also includes a score for the non-coding transcript of the gene.
The dots are just there as placeholders, they mean there was no value associated with that transcript. Many tools will show some sort of symbol instead of an empty field both for clarity and for practical technical reasons.
Visiting the variant's page on VarSome gives you a clearer picture since we collapse the identical scores and also include the converted rankscore provided by dbNSFP so you can have a single number for your variant:
Disclaimer: I work for the company behind VarSome, but it's a free tool. You need to pay to annotate VCF files (unlike the 100% free VEP), but it's free to use as a lookup tool for single variants.
answered Apr 10 at 12:06
terdon♦terdon
4,7752830
$endgroup$
The first gene you mention, ENSG00000196924, actually has 6 transcripts (link to the VarSome.com page of variant rs371839875), not 5. It's just that one of them is non-coding:
So the Sift scores you see are indeed one per transcript, it's just that there are 6 because dbNSFP also includes a score for the non-coding transcript of the gene.
The dots are just there as placeholders, they mean there was no value associated with that transcript. Many tools will show some sort of symbol instead of an empty field both for clarity and for practical technical reasons.
Visiting the variant's page on VarSome gives you a clearer picture since we collapse the identical scores and also include the converted rankscore provided by dbNSFP so you can have a single number for your variant:
Disclaimer: I work for the company behind VarSome, but it's a free tool. You need to pay to annotate VCF files (unlike the 100% free VEP), but it's free to use as a lookup tool for single variants.
answered Apr 10 at 12:06
terdon♦terdon
4,7752830
edited Apr 10 at 12:13
answered Apr 10 at 12:06
terdon♦terdon
4,7752830
answered Apr 10 at 12:06
terdon♦terdon
4,7752830
answered Apr 10 at 12:06
terdon♦terdon
4,7752830
4,7752830
add a comment |
add a comment |
Thanks for contributing an answer to Bioinformatics Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
Use MathJax to format equations. MathJax reference.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fbioinformatics.stackexchange.com%2fquestions%2f7430%2fvep-output-sift-score-unclear%23new-answer', 'question_page');
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
