SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
snpEff error sinnamone Bioinformatics 3 05-14-2014 11:05 PM
snpeff warnings mmmm Bioinformatics 2 01-23-2014 04:57 AM
snpEff annotation omerfaruk Bioinformatics 1 01-06-2014 12:25 AM
VarSifter and snpEFF v3 JKTeer Bioinformatics 0 01-15-2013 04:25 AM
snpEff error fc35802 Bioinformatics 19 02-22-2012 12:31 AM

Reply
 
Thread Tools
Old 04-04-2014, 03:42 AM   #1
mmmm
Senior Member
 
Location: UK

Join Date: Jul 2013
Posts: 131
Default snpeff

I have determined the snps (after mapping to the reference) and got vcf4 file with snps position within the chromosome
then I used snpEff to look at the effect of these snps BUT which output file includes the position of the snp and its effect (as in the snpeff.txt- only the effect of snps is listed within the gene but it does not refere to the snp position)??
mmmm is offline   Reply With Quote
Old 04-04-2014, 05:06 AM   #2
blakeoft
Member
 
Location: Connecticut

Join Date: Oct 2013
Posts: 79
Default

You should get a new vcf file as output. It should have all of the info of the original vcf, plus some new information in the INFO field following this format

Quote:
EFF= Effect ( Effect_Impact | Functional_Class | Codon_Change | Amino_Acid_Change| Amino_Acid_Length | Gene_Name | Transcript_BioType | Gene_Coding | Transcript_ID | Exon_Rank | Genotype_Number [ | ERRORS | WARNINGS ] )
See http://snpeff.sourceforge.net/SnpEff_manual.html#output for more details.

EDIT: According to the same page, it looks like you need to redirect the output. Their example:

Quote:
# Go to the directory where snpEff is installed
cd snpEff

# Do this only if you don't already have the database installed.
java -jar snpEff.jar download -v GRCh37.66

# Annotate the file
java -Xmx4g -jar snpEff.jar eff -v GRCh37.66 file.vcf > file.eff.vcf

Last edited by blakeoft; 04-04-2014 at 05:12 AM. Reason: added more useful information
blakeoft is offline   Reply With Quote
Old 04-04-2014, 06:36 AM   #3
mmmm
Senior Member
 
Location: UK

Join Date: Jul 2013
Posts: 131
Default

yes- thanks
I have got this, so this snp is synonomas in coding region but which gene/ region??- and what other information to be extracted?

DP=18;VDB=1.751185e-01;AF1=1;AC1=2;DP4=0,0,7,11;MQ=60;FQ=-81;EFF=DOWNSTREAM(MODIFIER||4283||257|yaaA||CODING|NP_459010.1||1),DOWNSTREAM(MODIFIER||576||21|thrL||CODING|NP_459006.1||1),SYNONYMOUS_CODING(LOW|SILENT|gtT/gtC|V165|820|thrA||CODING|NP_459007.1|1|1),UPSTREAM(MODIFIER||1970||309|thrB||CODING|NP_459008.1||1),UPSTREAM(MODIFIER||2903||428|thrC||CODING|NP_459009.1||1)
mmmm is offline   Reply With Quote
Old 04-04-2014, 06:54 AM   #4
blakeoft
Member
 
Location: Connecticut

Join Date: Oct 2013
Posts: 79
Default

It looks like snpEff has reported five different effects that this SNP could have. The gene names are in the sixth field inside each parentheses. They are yaaA, thrL, thrA, thrB, and thrC. The specific location of the variant should be in second column of the record. What you have posted is only the seventh (INFO) column.

Look back at http://snpeff.sourceforge.net/SnpEff_manual.html#output for more information about what snpEff adds to the INFO column.
blakeoft 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 08:43 PM.


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