<div dir="ltr"><div>I received a weird @vep email, and now i m unsure of whether this was published (cannot see it in the archive), so I send it again. Apologies if this is not the case.</div><div><br></div><div>=</div><div><br></div>Hello there,<div><br></div><div>I m writing with further VEP whishes! The context is that I m matching the  variants given to VEP with the content of a given database(s) that report variant effects (e.g. biomarkers of drug response). I m currently passing these knowledgebases to VEP as a vcf (with the corresponding .tbi) with the --custom flag, which works lovely.</div><div><br></div><div>The thing is that I m also interested to have a match when the nucleotide change of the input variant is not the same than the one from the database *but* it leads to the same aminoacid change; something like:</div><div><br></div><div>database:  chr1:xxxxxA>C   (geneX p.V28E)</div><div><br></div><div>input:  chr1xxxxxxA>C   (geneX p.V28E)  produces a 'perfect match'</div><div>input: chr1yyyyyyyG>T    (geneX p.V28E) produces an 'aminoacid match'</div><div><br></div><div>My plan is to keep using your --custom flag to retrieve the 'perfect' match and then to have my own script to retrieve the 'aminoacid match'. However, as you know, this is kind of a hassle (need to map the variants of the original databases with VEP to ensure that the protein statements are fully compatible, and I need to do it each time there is an update etc etc).</div><div><br></div><div>In other words, it would be more neat and easy if this 'aminoacid match' is provided by VEP given the --custom vcf file(s). By any chance, do you have some option (that I did not see) that magically does it, or any option that you can implement this in some incoming release?</div><div><br></div><div>(as for other requests, I think that this is a feature that would be useful for other users)</div><div><br></div><div>Thanks a lot in advance!</div><div>br</div></div>