Snps annotation using snpEff or R

I have the following data in VCF format I want to annotate the Snps I have in this file What is the code needed in snpEff or the code by R

#CHROM  POS ID  
    chr7    78599583    rs987435
    chr15   33395779    rs345783
    chr1    189807684   rs955894
    chr20   33907909    rs6088791
    chr12   75664046    rs11180435
    chr1    218890658   rs17571465
    chr4    127630276   rs17011450
    chr6    90919465    rs6919430
    chr16   5748791     rs2342723

snp gene r

• 4,241 views • link

updated 9.6 years ago by zx8754

• written 9.6 years ago by mms140130

snpEff is very well documented, please read and try some example: http://snpeff.sourceforge.net/SnpEff_manual.html

• 0 views • link

written 9.6 years ago by zx8754

I answered question like this here A: What is the best method to map snps called using GATK to corresponding genes?

• 0 views • link

written 9.6 years ago by Medhat

0 answers

No answers yet.

More posts like this

How can I get this sorted by ascending order? (when the number and character is together)
written 7.1 years ago by woojoy14 •

samtools view NA12878_dedup.bam | cut -f 3 | sort -u > chr.list8 > chr1 chr10 chr11 chr12 chr13 chr14 chr15 chr16 chr17 chr18 chr19 chr2 …
Removing a chromosome from a fasta file
written 8.3 years ago by t86dan •

Hello, I have a fasta file for human Genome GRCh37 Reference Assembly, but for some reason chromosome 20 is repeated and I would like to …
difference between 2 snp's data annotations
written 8.8 years ago by mms140130

Hi, I have annotation from the http://www.scandb.org/ , The SCAN website , and from data I got from my advisor the Scan annotation as follows …
MAFTOOLS HELP - Mutational Signature
written 8.9 years ago by David_emir

Hi All, I am trying to run mutational signatures on my MAF file, using maftools (https://github.com/PoisonAlien/maftools). while trying to run this i am constantly My …
finding the gene contain the snps ID
written 9.6 years ago by mms140130

I have the following data set about the snps ID CHROM POS ID chr7 78599583 rs987435 chr15 33395779 rs345783 chr1 189807684 rs955894 chr20 33907909 rs6088791 …
Finding the Genes corresponding to given SNPs using python
written 9.6 years ago by mms140130

Hi, my Question is how to find the genes related to given SNPs the part of Snps data I have as follows: > rs987435 C …
coverageBed with -g option ?
written 10.9 years ago by bioguy24

I am running a coverageBed and getting the below error. ``` cmccabe@DTV-A5211QLM:~/Desktop/NGS$ coverageBed -d -sorted -a /home/cmccabe/Desktop/NGS/bed/bedtools/xgen_targets_sorted.bed -b /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam > /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_output.txt ERROR: Database file /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam …
problem with methylKit of the annotation
written 11.6 years ago by hua.peng1314

I am dealing with RRBS data follow [the manual][1]. I met some problem when doing annotation.The problem showed below. ``` > read.transcript.features("mm9.refseq.genes.bed",remove.unsual=TRUE,up.flank=1500,down.flank=500,unique.prom=TRUE) > diffAnn <- …
CONTRA problem. It does not stop.
written 12.1 years ago by mafonso •

Hi, I am trying to run an algorithm for CNV detection and was testing CONTRA. The test file works but using my own files the …
CONTRA for CNV detection. troubleshooting
written 12.2 years ago by Kizuna

Hi Biostars Users, I would like to know if you have encountered this issue while working with [CONTRA tool][1] for detecting CNVs in targeted NGS. …