SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
mpileup SNP detection SWP Bioinformatics 2 05-06-2011 06:33 AM
Coverage required for SNP detection? Kasycas De novo discovery 3 04-12-2011 10:38 PM
SNP detection with gsMapper without reference dschika Bioinformatics 10 04-03-2011 08:00 PM
coverage and SNP detection rudi283 Bioinformatics 0 03-29-2011 07:22 AM
SNP detection doxologist Genomic Resequencing 14 09-06-2010 06:49 AM

Reply
 
Thread Tools
Old 07-23-2010, 12:43 AM   #1
Oxygen81
Member
 
Location: Germany

Join Date: Aug 2009
Posts: 11
Default SNP Detection

Hello everybody,

I am searching for some software with which it is possible to align 454 contigs against a reference sequence and to detect SNPs. It is essential that the quality file is considered and that a table with Position, SNP and quality is written in the output. Has anyone an idea?

Thank you in advance
Janina
Oxygen81 is offline   Reply With Quote
Old 07-23-2010, 02:57 AM   #2
nickloman
Senior Member
 
Location: Birmingham, UK

Join Date: Jul 2009
Posts: 356
Default

Try gsMapper (part of Newbler) supplied by Roche.
nickloman is offline   Reply With Quote
Old 07-23-2010, 06:33 AM   #3
krobison
Senior Member
 
Location: Boston area

Join Date: Nov 2007
Posts: 747
Default

The software wiki lists a number of packages specifically in this space (454 SNP Discovery); you might also look at the other Snp Discovery packages as many are platform agnostic.
krobison is offline   Reply With Quote
Old 07-23-2010, 08:05 AM   #4
SoftGenetics
Registered Vendor
 
Location: pa

Join Date: Apr 2009
Posts: 32
Default

Quote:
Originally Posted by Oxygen81 View Post
Hello everybody,

I am searching for some software with which it is possible to align 454 contigs against a reference sequence and to detect SNPs. It is essential that the quality file is considered and that a table with Position, SNP and quality is written in the output. Has anyone an idea?

Thank you in advance
Janina
NextGENe will do a great and easy job on this, also has a quality ranking, like Phred, for all found variants including indels
SoftGenetics is offline   Reply With Quote
Old 07-25-2010, 11:06 PM   #5
geschickten
Member
 
Location: India

Join Date: Jul 2009
Posts: 31
Default

Hi Janina, do you have gsMapper or any other similar s/w? if not write to me; I can probably help you...
geschickten is offline   Reply With Quote
Old 07-26-2010, 09:23 AM   #6
Oxygen81
Member
 
Location: Germany

Join Date: Aug 2009
Posts: 11
Default

Hello,

thank you all for your fast answers! I already tried CLC genomic workbench and Genious which are great but they don't display the quality of the detected SNPs. Maybe I have to make clear that I don't mean the quality of the SNP but of the sequenced base (for 454: 0-64). I asked our sequencing facility about the use of the gsMapper but they told me that it is not possible to import assembled contigs and no quality files, too... Is that right?
I will try NextGENe if there is a trial version available. Thank you all!
Oxygen81 is offline   Reply With Quote
Old 08-12-2010, 06:13 PM   #7
Nomijill
Member
 
Location: Southwest Florida

Join Date: Sep 2009
Posts: 24
Default Display of confidence scores

You should be able to view the quality scores of all of your bases in the Genomics Workbench. It is a display option. If you want, I can show you how.
Nomijill is offline   Reply With Quote
Old 08-12-2010, 06:17 PM   #8
link1
Member
 
Location: US

Join Date: May 2010
Posts: 21
Default

agreed gsmapper
link1 is offline   Reply With Quote
Old 08-13-2010, 12:43 AM   #9
Dethecor
Member
 
Location: Germany

Join Date: May 2010
Posts: 24
Default Why gsmapper?

Hi all,

just a question that popped into my mind. (I normally use output from Illumina machines)

Why do you need to use this gsmapper tool instead of a standard SNP-calling pipeline like:

sequencer -> read-file -> [bwa|bowtie|...(aligner of your choice)] -> output.sam-file

and then just run 'samtools pileup' on your output file and get the SNP's with quality and all from there?

I was just wondering if there is something special about 454-reads, that renders this approach unfeasible?

Cheers,
Paul
__________________

"You are only young once, but you can stay immature indefinitely."
Dethecor is offline   Reply With Quote
Old 10-07-2010, 01:12 PM   #10
sklages
Senior Member
 
Location: Berlin, DE

Join Date: May 2008
Posts: 628
Default

Quote:
Originally Posted by link1 View Post
agreed gsmapper
gsMapper will align 454 reads and will refuse mapping sequences longer than a certain threshold, like contigs.
sklages is offline   Reply With Quote
Old 10-07-2010, 01:17 PM   #11
sklages
Senior Member
 
Location: Berlin, DE

Join Date: May 2008
Posts: 628
Default

Quote:
Originally Posted by Dethecor View Post
Hi all,

just a question that popped into my mind. (I normally use output from Illumina machines)

Why do you need to use this gsmapper tool instead of a standard SNP-calling pipeline like:

sequencer -> read-file -> [bwa|bowtie|...(aligner of your choice)] -> output.sam-file

and then just run 'samtools pileup' on your output file and get the SNP's with quality and all from there?

I was just wondering if there is something special about 454-reads, that renders this approach unfeasible?

Cheers,
Paul
454 reads are longer and are not fixed-length; the short read aligners are not the approbiate tools for that kind of sequencing data.

And, he was asking to map contigs, not reads. Although, I would map the reads and look for variations .. :-)
sklages is offline   Reply With Quote
Old 10-08-2010, 05:40 AM   #12
SoftGenetics
Registered Vendor
 
Location: pa

Join Date: Apr 2009
Posts: 32
Default

You may wish to test the aligner in NextGENe which also provides a high level of visualization and annotation.
SoftGenetics is offline   Reply With Quote
Old 10-08-2010, 11:50 AM   #13
DNASTAR
Registered Vendor
 
Location: Madison, WI

Join Date: Aug 2010
Posts: 48
Default

Hi Janina,

Lasergene's SeqMan app can assemble 454 data against a reference, and detect SNPs. The SNP reports will display the imported quality scores, position, SNP base, as well as a lot of other useful info. You can learn more or get a free trial here.

Thanks,
Anne
DNASTAR 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 05:54 AM.


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