Just for completeness, I'll point out snpEff (http://snpeff.sourceforge.net/, and referenced in the Biostars link above) which achieves similar goals. I found snpEff to be a bit easier to get working for non-human data, for what it's worth.
Determine Variant Location From Vcf
hi all, i have a rather large VCF files with about 2.8 million variations. i want to know if each of the variation is within an exon, intron or promoter region?
i have downloaded the corresponding files from RefSeq and I wanted to know if i could use BED Intersect to accomplish this?
can anyone help with the syntax especially i need the output file in a form that i can parse easily and which will tell me if the variant was exonic or intronic or in the promoter region.
thanks in advance
3 answers
You can use the ensembl variant effect predictor
http://www.ensembl.org/Homo_sapiens/UserData/UploadVariations
http://www.ensembl.org/tools.html
To work this out, it will take vcf both for the web interface or the perl script (for a large number of variants the perl script is required) and compare to any given species ensembl database
you can use snpEff wich a felexible software you can use any genome you want since you have a gff file, vcf file and a reference genome. I hope that it will help you. http://snpeff.sourceforge.net/
• 0 views
•
link
written
by
Mchimich •
To annotate this kind of stuff in your VCF upload it to SeattleSeq, or you can download and install a tool like Annovar to accomplish the same thing.These tools will annotate the coding and noncoding regions, and then from there you can go after details like promotor regions, etc.
You can also read more about these here.
• 0 views
•
link
written
by
Alex Paciorkowski
Yeah, kind of assumed the OP was referring to human data of course. Maybe the OP can clarify what organism the data is from?
• 0 views
•
link
written
by
Alex Paciorkowski
Log in to answer this question.
More posts like this
-
Intron variant identifier
written by Sabeen •
-
Extract reads mapped within chosen intronic region from BAM file (no junctions)
written by gspirito •
-
Extract only RNA-seq reads in a BAM file that extend into an intron
written by foxdie •
-
Extract codon information from VCF file
written by k.rajain1212 •
-
Help regarding variant filtration and variant density calculation
written by mirza
-
somatic mutations from WGS of CGI
written by kanwarjag
-
How To Understand Some Annotated Terms Made By Annovar
written by zengtony743
-
A Question About Annotation To Exome Seq
written by camelbbs
-
Promoter Or Enhancer Regions Bed Format
written by win
-
Determine Where Is My Oligo Probe In The Gene? In Exonic Region Or Interonic Region
written by Omid
Unfortunately I won't be able to flesh out a more complete answer, but have you taken a look at using a tool like 'annovar' to annotate your VCF file?
http://www.openbioinformatics.org/annovar/