SNP annotated outside gene location (Variant Effect Predictor)
Entering edit mode
5.2 years ago
serpalma.v ▴ 70

I have the following SNP annotated with variant effect predictor to two genes, their locations in the genome are shown as well (UMD3.1):

  • rs469441731: chr1:83,591,886-83,591,886
  • ENSBTAG00000020106: chr1: 83,583,028-83,591,912 (reverse strand)
  • ENSBTAG00000019336: chr1: 83,595,844-83,609,419 (reverse strand)

I see that the location of rs469441731 lies within the boundaries of ENSBTAG00000020106 but not of ENSBTAG00000019336.

I wonder why is rs469441731 annotated to ENSBTAG00000019336

vep SNP variant effect predictor ensembl • 1.1k views
Entering edit mode

Thinking about it quickly, there is the (default) option Upstream/Downstream distance (bp): 5000, which checks the flanking regions of a variant for overlap with a transcripts. As your third range lies within 5kb, maybe that is the reason. Therefore, please provide the full output and the command line used.

Entering edit mode
5.2 years ago
Denise CS ★ 5.2k

You need to zoom out to view the boundaries of ENSBTAG00000019336. Variant rs469441731 is 5' UTR with regards to ENSBTAG00000020106, but downstream in regards to ENSBTAG00000019336.

This is the results of the VEP seen on the Location tab of the Ensembl Browser. Note the new data track called 'VEP result' and the variant as a grey box right in the middle of that view. The grey box maps right in the middle of the first exon of ENSBTAG00000020106 (i.e. AP2M1, which is a reverse stranded gene).

You need to zoom out to view a larger region (e.g. 6 kb either side of the SNP of interest), to see the same rs469441731 is downstream of ENSBTAG00000019336 (i.e. DVL3, another reverse stranded gene). This is the zoomed out view.


Login before adding your answer.

Traffic: 1029 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6