SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Reply
 
Thread Tools
Old 06-18-2014, 02:26 AM   #1
mmmm
Senior Member
 
Location: UK

Join Date: Jul 2013
Posts: 131
Default SNPs_analysis

Dear all, I have 2 draft genomes of bacteria assembled by velvet and I used Mauve to extract snps (but I need your advice to:
- know if these snps are synonomas or nonsynonomas
-in which genes they are located?

is there a softwrae that can do this?

(I know snpeff- can be used but I do not have a reference strain)

your advice is very appreciated
mmmm is offline   Reply With Quote
Old 06-19-2014, 05:03 AM   #2
TiborNagy
Senior Member
 
Location: Budapest

Join Date: Mar 2010
Posts: 329
Default

First you need to create an annotation (RAST) and after that you can use snpeff.
TiborNagy is offline   Reply With Quote
Old 06-27-2014, 05:47 AM   #3
mmmm
Senior Member
 
Location: UK

Join Date: Jul 2013
Posts: 131
Default

thanks- have created the annotation using RAST on the draft genome- now how to us the annotation file in snpeff (can I download the genbank file to snpeff database)???
mmmm is offline   Reply With Quote
Old 06-27-2014, 12:02 PM   #4
BurlEarl
Member
 
Location: Salt lake city, UT, US

Join Date: Jun 2012
Posts: 19
Default

This is a pain. To add a custom genome to SnpEff you need to add lines to the snpEff.config file. There is good documentation here:

http://snpeff.sourceforge.net/SnpEff...html#databases

Under building databases. It is a little picky.

i added lines to snpEff.config like this:

#Friend Murine Leukemia Virus
X02794.genome : FMuLV

cd PATH/TO/snpEff/data
mkdir X02794
cp PATH/TO/GENBANKFILE/genes.gb X02794/genes.gb

java -jar snpEff.jar build -genbank -v X02794

I think that is what i did...Anyway, documentation is pretty good. I remember having a problem with my generated .gb file. Turns out there needed to be a couple of spaces after ORIGIN and the sequence.

good luck
BurlEarl 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:21 AM.


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