SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
454 quality filtering GSCHALLA 454 Pyrosequencing 0 01-30-2012 02:57 PM
Should I do any quality filtering in such a case? GloriaFu Bioinformatics 4 06-10-2011 11:51 AM
454 amplicon quality filtering JackieBadger Bioinformatics 10 03-16-2011 01:28 PM
Quality filtering Farhat SOLiD 0 05-13-2010 11:42 PM
Read quality filtering for long, PE runs kmcarr Illumina/Solexa 0 07-21-2009 07:07 AM

Reply
 
Thread Tools
Old 05-18-2010, 04:12 AM   #1
Farhat
Member
 
Location: Pune, India

Join Date: Apr 2008
Posts: 21
Default Filtering on quality

What kind of trimming and/or quality filter do you use before mapping reads? Are there any commonly accepted filtering criteria?

I am trying to map some data from an older sequencer and my filtering (Q>15, no more than 6 bases outside quality range allowed) seems to filter out about 85% of the reads.
__________________
Farhat Habib
Farhat is offline   Reply With Quote
Old 05-18-2010, 10:37 PM   #2
Jose Blanca
Member
 
Location: Valencia, Spain

Join Date: Aug 2009
Posts: 70
Default

I use a combination of several tools lucy, exonerate, blast, Univec, quality trimming. If you want to take a look I wrote a document about our read cleaning process.
Jose Blanca is offline   Reply With Quote
Old 05-19-2010, 12:32 AM   #3
strob
Member
 
Location: Belgium

Join Date: Nov 2008
Posts: 79
Default

fastx toolkit (incorporated in Galaxy)
strob is offline   Reply With Quote
Old 05-19-2010, 03:12 AM   #4
Farhat
Member
 
Location: Pune, India

Join Date: Apr 2008
Posts: 21
Default

Thanks. I am aware of software that does this, what I would like to know is the criteria used in filtering, (e.g., Q>15, maximum six bases outside the quality range, etc.)
__________________
Farhat Habib
Farhat is offline   Reply With Quote
Old 05-19-2010, 06:24 AM   #5
Jose Blanca
Member
 
Location: Valencia, Spain

Join Date: Aug 2009
Posts: 70
Default

Our default pipelines depend on the kind of sequences (mainly on the sequencing platform). Each pipepile has several steps (you can look for them in the ngs_backbone page). For these steps our defaults are:

For the unknown vectors Blast against the UniVec database. We keep the hits with more than 15 bases and a similarity above 96%. We remove this regions from the reads.

For the adaptors longer than 15 bases, exonerate alignment. Again we keep the hits with more than 15 bases and a similarity above 96%.

We look for the adaptors shorted than 15 bases using regular expressions.

Quality for short sequences. We apply a quality threshold of 20 and a window of 1. We keep the regions with more than 3 good bases.

Quality for the long sequences. We use lucy with the default parameters. If we have the cloning site we use it.

For the sequences without quality we use trimpoly to look for N and remove the regions with too many Ns.

Sometimes, we remove 3 or 4 nucleotides from the edges.
Jose Blanca 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 02:09 PM.


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