SVA only support hg18 so no point using it for the moment if you aligned to hg19.
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
I'm an intern working with exome analysis and I am facing the same thing. I want to implement this annotation into my own java pipeline and although the solution may seem very easy, I am troubled in finding the right approach.
I have tried using the RefSeq annotation file of SeqCap EZ Exome v2 (matching UCSC genome browser with HG19), which holds information on cdsStart, -End en exon starts and endings. This file also holds an Ensembl gene reference for each RefSeq gene, which should make it easy to link with the cDNA fasta file of Ensembl and get exactly what I want...
... a few problems though:
1) RefSeq- and ensembl genes overlap and multiple of the same ensembl references may occur in the ensembl fasta file, making it hard to differentiate. This is most likely due to different isoforms.
2) Looking at a few cases, I noticed that some RefSeq genes show cdsStart and cdsEnd positions that can not be traced back to ensembl. In other words: when I read the ensembl reference from the RefSeq file and look them up in the ensembl file, I can find multiple isoforms, but none with the same cdsStart and/or cdsEnd. I already take into account that RefSeq and Ensembl differ 1 nuc. in cdsStart. Both files are based on HG19, so that can't be the problem either.
What would be the best approach on solving this puzzle? Should I just walk through the entire genome and annotate all the information to my SNPs as I go along? Any thoughts are welcome.
Thanks a bunch!
Comment
-
If you have hg19/GRCh37 positions for all your snps I would suggest using a tool like the ensembl variant effect predictor to get the consequences of your snps and then tracing to refseq ids using the ensembl xref system rather than doing it the other way around
Refseq models and Ensembl models should be mostly the same for the cds coordinates (though not in all models) but to get the models which are identical across both sets it best to look at the ccds models http://www.ncbi.nlm.nih.gov/projects...CcdsBrowse.cgi
Do remember that utr coordinate may be different across both sets
Comment
-
We're looking at using Ensembl's Variant Effect Predictor right now ... but with a non-model organism (or at least a model organism with some regions replaced by a "better" assembly). It seems like we'll need to set up a local database in order to provide the reference and gene models. Does anyone have experience setting up a db like this?
Comment
-
You could also try :
The requirement is that the input snps be in VCF format.
-Abhi
Comment
-
For help with setting up a custom ensembl database I would suggest emailing [email protected] to get help
Comment
-
Originally posted by jnfass View PostThanks, but it seems to only do SNPs ... we need to annotate indels, too.Mendelian Disorder: A blogshare of random useful information for general public consumption. [Blog]
Breakway: A Program to Identify Structural Variations in Genomic Data [Website] [Forum Post]
Projects: U87MG whole genome sequence [Website] [Paper]
Comment
-
@Michael.James.Clark ... Annovar seems to require annotation databases from the UCSC Genome Browser. That doesn't exist for at least part of the genome I'm working with.
@nexgengirl ... Polyphen-2 seems to be restricted to human
@laura ... thanks for the suggestion; I might try contacting ensembl.
But in the meantime someone responded to my Biostar post and suggested snpEff (snpeff.sourceforge.net) and it seems to fit the bill. If anyone has had good or bad experience with it, I'd appreciate hearing about it.
Comment
-
Originally posted by jnfass View Post@Michael.James.Clark ... Annovar seems to require annotation databases from the UCSC Genome Browser. That doesn't exist for at least part of the genome I'm working with.
@nexgengirl ... Polyphen-2 seems to be restricted to human
@laura ... thanks for the suggestion; I might try contacting ensembl.
But in the meantime someone responded to my Biostar post and suggested snpEff (snpeff.sourceforge.net) and it seems to fit the bill. If anyone has had good or bad experience with it, I'd appreciate hearing about it.Mendelian Disorder: A blogshare of random useful information for general public consumption. [Blog]
Breakway: A Program to Identify Structural Variations in Genomic Data [Website] [Forum Post]
Projects: U87MG whole genome sequence [Website] [Paper]
Comment
-
Yeah, I'm not saying it is necessarily the easiest solution for you, but it is incorrect to state that it requires annotation databases from UCSC.Mendelian Disorder: A blogshare of random useful information for general public consumption. [Blog]
Breakway: A Program to Identify Structural Variations in Genomic Data [Website] [Forum Post]
Projects: U87MG whole genome sequence [Website] [Paper]
Comment
Latest Articles
Collapse
-
by seqadmin
The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...-
Channel: Articles
02-24-2025, 06:31 AM -
-
by seqadmin
Like all molecular biology applications, next-generation sequencing (NGS) workflows require diligent quality control (QC) measures to ensure accurate and reproducible results. Proper QC begins at nucleic acid extraction and continues all the way through to data analysis. This article outlines the key QC steps in an NGS workflow, along with the commonly used tools and techniques.
Nucleic Acid Quality Control
Preparing for NGS starts with isolating the...-
Channel: Articles
02-10-2025, 01:58 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 03-03-2025, 01:15 PM
|
0 responses
46 views
0 likes
|
Last Post
by seqadmin
03-03-2025, 01:15 PM
|
||
Started by seqadmin, 02-28-2025, 12:58 PM
|
0 responses
167 views
0 likes
|
Last Post
by seqadmin
02-28-2025, 12:58 PM
|
||
Started by seqadmin, 02-24-2025, 02:48 PM
|
0 responses
525 views
0 likes
|
Last Post
by seqadmin
02-24-2025, 02:48 PM
|
||
Started by seqadmin, 02-21-2025, 02:46 PM
|
0 responses
256 views
0 likes
|
Last Post
by seqadmin
02-21-2025, 02:46 PM
|
Comment