The 1000 genomes data is available with consequence annotations here: http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20130502/supporting/functional_annotation/. I then simply parse the variants for those that have a missense or synonymous consequence annotation.
1000 genomes coding SNPs not in exonic regions?
Hi, I am attempting to identify the exon that each of my synonymous or missense SNPs in the 1000 genomes data belongs to. I am using the GENCODE GTF files found here: https://www.gencodegenes.org/human/ and extracting all exons.
I then use bedtools to identify which exon each of my SNPs fall in. It appears that many of my SNPs' co-ordinates are not within any exon. What I would like to know is if and how synonymous or missense SNPs can fall in intronic regions?
1 answer
Log in to answer this question.
More posts like this
-
Identifying introns
written by bioinformatics.queries
-
Does bedtools intersect output overlaps with isoforms?
written by Kash
-
How to extract upstream and downstream genes of a list of locations from gtf file?
written by aboozar.soorni •
-
How to make an intronic GTF for alignment
written by piyushjo
-
extracting the start and end co-ordinates of a gene from ucsc table browser
written by cs18s003 •
-
Discrepencies between ANNOVAR and Gencode annotations of 1000 genomes
written by spiral01
-
Denisovan and Altai Neanderthal .bed files?
written by spiral01
-
Denisovan and Altai Neanderthal .bed files?
written by spiral01
-
Going from Snps to exons to protein domains
written by AB
-
Exon falls within a gap in the assembly
written by PoGibas
Can you confirm that the reference genomes are the same, so hg19 vs hg19 or hg38 vs hg38?
Hi, thanks for the response. Yes I can confirm that the ref genomes are the same, hg38.
Where are you getting your 1K genome SNPs from?
I am getting the data with GENCODE annotations here: http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20130502/supporting/functional_annotation/
Could these be artifacts from a liftOver operation, perhaps?
Why are you comparing to the GTF? There are tools designed to do exactly what you need.
I need to obtain the exon that each SNP lies in, as well as the start and end co-ordinates (because my ultimate goal is to identify the length of the specific exon that each SNP lies in). The available GENCODE annotation of 1000 genomes variants provides the exon number within the gene, but not the exon id or start and end coordinates?
Simply get the gencode annotation for hg19, extract exons, and use
bedtools intersectwhere-ais the SNPs and-bis the exon.gtf. Use option-wbto return the entire interval of the matching exon. From there you can cut or awk out what you need.