How to get the gene name/region from a list of positions
1
0
Entering edit mode
9 months ago

I have a list of SNP positions:

761155 
4402247
472711
4027247
2867693

I have a gff reference annotation file and a reference fasta file.

How do I simply find the gene name and locus (e.g. "ctpJ", "Rv3743c", respectively; or "Rv3760", "Rv3760" if there is no gene name) where these SNPs occur?

SNP gene • 253 views
ADD COMMENT
1
Entering edit mode

Position data should include both chromosome and position within chromosome. Once you have that, you can use biomaRt to get gene data. Please search the site - that topic has been addressed multiple times.

ADD REPLY
0
Entering edit mode
9 months ago

Seems like you're looking for the Variant Effect Predictor. The link I've provided is for Mycobacterium_tuberculosis, which is what I presume you want based on the gene name provided. The tool works for any genome in Ensembl though.

ADD COMMENT

Login before adding your answer.

Traffic: 1638 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6