SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > 454 Pyrosequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Show filtered reads in the IGV? MBender Bioinformatics 0 01-25-2012 12:04 PM
Number of Reads Rfriedman Bioinformatics 0 01-11-2012 01:43 PM
quality filtered illumina PE reads Wallysb01 Bioinformatics 1 07-21-2011 10:04 AM
number of reads m_elena_bioinfo Bioinformatics 2 07-20-2010 08:47 AM
getting reads not filtered? gerald2545 454 Pyrosequencing 12 05-18-2009 12:00 AM

Reply
 
Thread Tools
Old 09-14-2011, 04:31 AM   #1
Dynamac
Member
 
Location: Cincinnati, OH

Join Date: Aug 2011
Posts: 46
Default Inordinate Number of Reads Being Filtered

I am losing about half of my data from the TrimBack Valley and Signal Intensity Filters (numTrimmedTooShortQuality). Has anyone had this problem, and if so, what are your solutions and/or recommendations?
Any comments at all would be appreciated, as I am still a newbie with this system, and support from Roche has been spotty at best.

Thanks!

Barry
Dynamac is offline   Reply With Quote
Old 10-11-2011, 09:10 AM   #2
vamosia
Member
 
Location: New York

Join Date: Mar 2009
Posts: 15
Default

Which pipeline are you using? The ends of the reads may be of lower quality, which may caused reads to be filtered out.

If your still interested in the those rejected reads anyway, you could try the standard shotgun pipeline (which supposedly is less stringent compare to the amplicon pipeline)
vamosia is offline   Reply With Quote
Old 10-11-2011, 09:32 AM   #3
Dynamac
Member
 
Location: Cincinnati, OH

Join Date: Aug 2011
Posts: 46
Default

Yes, we tried the shotgun pipeline and that roughly doubled our number of bases called.

Thanks!
Dynamac is offline   Reply With Quote
Old 10-11-2011, 09:59 AM   #4
vamosia
Member
 
Location: New York

Join Date: Mar 2009
Posts: 15
Default

Great, keep in mind that the amplicon pipeline does several additional correction to help reduce homopolymer error (amplicon library are generally very bright, hence more chance of producing homopolymer error).

While the using the shotgun pipeline to process amplicon samples gives more data, it may or may contain more homopolyer error. Just something to watch for.
vamosia is offline   Reply With Quote
Old 10-24-2011, 02:04 PM   #5
ajthomas
Senior Member
 
Location: Utah

Join Date: Mar 2010
Posts: 166
Default

In my experience, the extra data you get by processing amplicon data with the shotgun pipeline doesn't help. The extra data doesn't ususally equate to more information because the lower quality of the data offsets the reduced quantity. As for homopolymers, I do find more errors in homopolymer regions. Of course, your experience may differ from mine. It's probably worth it to try it both ways and see if that extra data helps or not.

There are a few TCBs relating to amplicon sequencing that you should read if you haven't already. There's one (sorry, I don't remember which one off the top of my head) that talks about the what makes the amplicon filtering pipeline different. You should certainly read that if you haven't already.
ajthomas is offline   Reply With Quote
Old 10-25-2011, 04:03 AM   #6
pmiguel
Senior Member
 
Location: Purdue University, West Lafayette, Indiana

Join Date: Aug 2008
Posts: 2,317
Default

But the amplicon pipeline does not do quality trim-back. It tosses reads that don't meet its quality criteria, but the reads it keeps may still have low quality 3' sequence in them.

--
Phillip
pmiguel is offline   Reply With Quote
Old 10-25-2011, 03:09 PM   #7
RCJK
Senior Member
 
Location: Australia

Join Date: May 2009
Posts: 155
Default

One of the recommended filtering strategies is to use the amplicon pipeline, but change the setting so that it will trim the reads instead of discarding them. This is in the latest (July 2011) guide to amplicon experimental design.
RCJK is offline   Reply With Quote
Reply

Tags
lost data, signal intensity, valley filter

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 11:08 PM.


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