SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   General (http://seqanswers.com/forums/forumdisplay.php?f=16)
-   -   Threshold quality score to determine the quality read of ILLUMINA reads problem (http://seqanswers.com/forums/showthread.php?t=6821)

edge 09-13-2010 06:55 AM

Threshold quality score to determine the quality read of ILLUMINA reads problem
 
Hi,

Does anybody have idea regarding the general threshold average quality score to determine whether the ILLUMINA raw reads is good or bad quality?

Below is one of my ILLUMINA raw reads in fastq format:

@HWI-EAS001
AACATTCAACGCTGNCGGTGAGTTGGAATTCTCGGGTGCCAAGGAACTCC
+HWI-EAS001
^`Y^aa__^_]\`_B\U][RV`W`^`][``__Z^P[UUZZUUa^Z[^^Z[

The above reads got 50 nucleotide in length.
Is it got any program or script able to calculate the average quality score of above read?
My purpose is hope to calculate the average quality score of each read.
Based on the average quality score of each read, I plan to filter out those "low quality reads" that below threshold of average quality score.
Thanks a lot for any advice.

svl 09-13-2010 02:22 PM

Maybe galaxy has what you're looking for:
http://main.g2.bx.psu.edu/

Or fastqc:
http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc/

Or have a look at the wiki:
http://seqanswers.com/wiki/Special:B...uality_Control


All times are GMT -8. The time now is 11:45 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.