<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 12 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:"Myriad Pro";
panose-1:2 11 5 3 3 4 3 2 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.MsoAcetate, li.MsoAcetate, div.MsoAcetate
{mso-style-priority:99;
mso-style-link:"Balloon Text Char";
margin:0in;
margin-bottom:.0001pt;
font-size:8.0pt;
font-family:"Tahoma","sans-serif";}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
span.BalloonTextChar
{mso-style-name:"Balloon Text Char";
mso-style-priority:99;
mso-style-link:"Balloon Text";
font-family:"Tahoma","sans-serif";}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoNormal>Hi,<o:p></o:p></p><p class=MsoNormal>I’m having an issue with the --pick filtering option using the variant_effect_predictor.pl script (version 75). For a certain variant (see below), it is unclear why this filter is “picking” the specified consequence over seemingly higher ranked consequences.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Here is the variant in .vcf format:<o:p></o:p></p><p class=MsoNormal>19 33793007 . TCGC T . PASS<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>The unfiltered/unpicked consequence is as follows:<o:p></o:p></p><p class=MsoNormal>CSQ=-|ENSG00000245848|ENST00000498907|Transcript|inframe_deletion|461-463|311-313|104-105|GD/D|gGCGac/gac|COSM18270&TMP_ESP_19_33793008_33793010||||||||-1||YES|CEBPA|HGNC||||protein_coding|Low_complexity_(Seg):Seg&PROSITE_profiles:PS50315&PIRSF_domain:PIRSF005879|CCDS54243.1|ENST00000498907.2:c.311_313delGCG|ENSP00000427514.1:p.Gly104del|||||,-|ENSG00000267130|ENST00000593041|Transcript|upstream_gene_variant||||||COSM18270&TMP_ESP_19_33793008_33793010|||||||2932|1||YES|CTD-2540B15.9|Clone_based_vega_gene||||lincRNA|||||||||,-|ENSG00000267727|ENST00000587312|Transcript|downstream_gene_variant||||||COSM18270&TMP_ESP_19_33793008_33793010|||||||162|1||YES|CTD-2540B15.7|Clone_based_vega_gene||||antisense|||||||||,-|ENSG00000178863|ENST00000320232|Transcript|upstream_gene_variant||||||COSM18270&TMP_ESP_19_33793008_33793010|||||||966|1||YES|CEBPA-AS1|HGNC||||pseudogene|||||||||,-|ENSG00000267580|ENST00000589932|Transcript|downstream_gene_variant||||||COSM18270&TMP_ESP_19_33793008_33793010|||||||934|1||YES|CTD-2540B15.11|Clone_based_vega_gene||||antisense|||||||||,-|ENSG00000267296|ENST00000592982|Transcript|upstream_gene_variant||||||COSM18270&TMP_ESP_19_33793008_33793010|||||||753|1||YES|CEBPA-AS1|Clone_based_vega_gene||||antisense|||||||||,-|ENSG00000230259|ENST00000425420|Transcript|intron_variant&nc_transcript_variant&feature_truncation||||||COSM18270&TMP_ESP_19_33793008_33793010||||||||-1||YES|AC008738.1|Clone_based_ensembl_gene||||pseudogene|||ENST00000425420.2:n.246+175_246+177delGCG||||||,-||ENSR00000345628|RegulatoryFeature|regulatory_region_variant||||||COSM18270&TMP_ESP_19_33793008_33793010|||||||||||||||||||||||||<o:p></o:p></p><p class=MsoNormal><br>Although hard to see, the first consequence is inframe_deletion of a protein_coding transcript. This seems to be the consequence that should be highest ranked. However, when I use --pick, I get the following consequence:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>CSQ=-|ENSG00000230259|ENST00000425420|Transcript|intron_variant&nc_transcript_variant&feature_truncation||||||COSM18270&TMP_ESP_19_33793008_33793010||||||||-1||YES|AC008738.1|Clone_based_ensembl_gene||||pseudogene|||ENST00000425420.2:n.246+175_246+177delGCG||||||<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>This is an intron_variant&nc_transcript_variant&feature_truncation of a pseudogene. I believe this should be of lower rank based on biotype and consequence type. <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Is there a reason why the --pick filter is choosing this consequence over other (seemingly more relevant) consequences? Is it because there are multiple consequences (intron_variant + nc_transcript_variant + feature_truncation) in the same consequence?<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Note, when I run this same variant through the VEP Web interface, I get the same consequence when selecting “Show one selected consequence per variant”.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Any help would be greatly appreciated!<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><span style='font-family:"Myriad Pro","sans-serif"'>Andrew R. Carson, Ph.D.</span> <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p></div></body></html>