View Single Post
Old 09-13-2010, 06:55 AM   #1
edge
Senior Member
 
Location: China

Join Date: Sep 2009
Posts: 199
Default 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.
edge is offline   Reply With Quote