SEQanswers

Go Back   SEQanswers > General



Similar Threads
Thread Thread Starter Forum Replies Last Post
VEP variant effect predictor DNAjunk Bioinformatics 0 11-12-2014 07:57 AM
SNP deleterious or tolerated? (VEP tool) clarissaboschi Bioinformatics 4 10-24-2014 05:50 PM
Some help with Variant Effect Predictor (VEP) moredd Bioinformatics 0 04-30-2014 12:21 PM
annovar cmccabe Bioinformatics 0 06-25-2013 02:46 PM
the result of VEP bewlib Bioinformatics 0 10-20-2012 01:12 AM

Reply
 
Thread Tools
Old 03-13-2015, 07:20 AM   #1
akvarel
Member
 
Location: Berlin

Join Date: Nov 2011
Posts: 11
Default ANNOVAR vs VEP

I was comparing the results from ANNOVAR and VEP and what I found is that


gene-based vs transcript based
ANNOVAR (ENSEMBL ids) seems to output information only gene-based, not transcript based, therefore out of 81 SNPs I have 81(+-) genes as those that fall in intronic/exonic/ncRNA has one gene annotation and others neighboring genes.
VEP outputs all ensembld ids of transcripts that are affected by the SNP.

1. Is there any possibility to force ANNOVAR to take all transcripts affected by SNP into consideration?

annovar line: annotate_variation.pl -build hg19 SNPs_annovar_input.avinput humandb/ -dbtype ensGene
vep: variant_effect_predictor.pl -i SNPs_vep_input.txt --cache --force_overwrite --symbol


intergenic-variantion handling
​If a variation is intergentic, ANNOVAR outputs the neighboring genes. VEP tells nothing about the neighboring genes.

2. ​Is there any possibility to force VEP to output the neighboring genes for an internecine variations?


strand information in VEP
Strand information should be specified for VEP input and is crucial for a prediction about how the variant affects AS. ANNULAR does not need this specification.
I do not have a strand information for my SNPs (they are taken from a paper and authors provided no information about the strand, only MAF, normal allele and risk allele).

3. Is there a possibility to assign each SNP (using a MAF and risk,normal allele information, position) a strand using R and some database (better command line approach) ?



Thanks

Last edited by akvarel; 03-13-2015 at 07:23 AM.
akvarel is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 07:02 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO