SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > 454 Pyrosequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
EBARDenovo - A new RNA-seq do novo assembler for paired-end Illumina data htchu.taiwan RNA Sequencing 2 06-10-2013 01:13 AM
EBARDenovo - A new RNA-seq do novo assembler for paired-end Illumina data htchu.taiwan Illumina/Solexa 9 04-16-2013 12:08 AM
parallel de novo assembler tmy1018 Bioinformatics 3 10-22-2012 09:31 AM
Input parameters GS de novo Assembler grassgirl 454 Pyrosequencing 0 06-03-2011 12:40 PM
De Novo Short Read Assembler? doxologist De novo discovery 18 05-21-2010 06:55 AM

Reply
 
Thread Tools
Old 09-23-2011, 11:27 AM   #1
grassgirl
Member
 
Location: Oregon, USA

Join Date: Mar 2011
Posts: 51
Default Groom data prior to using GS de novo Assembler?

Hello forum members,

I'm running GS de novo Assembler on 454 junior data and I'm wondering if folks groom the data based on qual scores before running the assembler. I can pass in the sff file, which contains the qual data, so I assume the assembler does something with the qual data. The qual scores on isotigs from a recent run were very good and I have also read from others that using the sff file produced a better assembly than simply passing the read file without the quals. I tried to get an answer from the manual (as far as qual cutoff for reads), but found none.

Any ideas?
grassgirl is offline   Reply With Quote
Old 09-26-2011, 04:12 AM   #2
flxlex
Moderator
 
Location: Oslo, Norway

Join Date: Nov 2008
Posts: 415
Default

I never 'groom' as newbler (gsAssembler) is doing this for me :-) Using the quality data from the sff file, as you wrote...
flxlex is offline   Reply With Quote
Old 09-26-2011, 10:02 AM   #3
grassgirl
Member
 
Location: Oregon, USA

Join Date: Mar 2011
Posts: 51
Default

Hi flxlex,

Thanks for your reply. Do you know what criteria newbler uses (e.g., a lower cutoff of qual 10 or something)?
grassgirl is offline   Reply With Quote
Old 09-27-2011, 02:02 AM   #4
flxlex
Moderator
 
Location: Oslo, Norway

Join Date: Nov 2008
Posts: 415
Default

Sorry, no clue...
flxlex is offline   Reply With Quote
Old 09-27-2011, 10:55 AM   #5
grassgirl
Member
 
Location: Oregon, USA

Join Date: Mar 2011
Posts: 51
Thumbs up

Well, I will try and find an answer and if I do I'll post it here.

Thanks.
grassgirl is offline   Reply With Quote
Old 09-29-2011, 03:47 PM   #6
grassgirl
Member
 
Location: Oregon, USA

Join Date: Mar 2011
Posts: 51
Default

I contacted Roche and this is what the rep told me:

"The quality scores are not used to eliminate reads from the GS Assembler - the only input reads that are eliminated are those that are 'Too Short' (less than 50 bases)

During the assembly, it will generate consensus basecalls of the contigs by using quality and flow signal information for each nucleotide flow included in the contigs' multiple alignments.

If you were to just give it the fasta file and no quality scores it would make the assumption that the quality is good.

You can change the stringency of the quality filters during signal processing if you wanted to eliminate reads."

The second paragraph is from the manual. Still a bit of a black box, but I'm guessing that a particular base is called from the highest qual-scoring read(s) and the consensus drawn from that, then the qual score for the isotig/contig is given.

I guess the thing to do is filter your isotigs/contigs based on qual scores prior to further processing. Have you done this before Blasting flxlex?
grassgirl 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 03:38 PM.


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