SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
Base qualities for Illumina Sequencing fongchun Bioinformatics 6 11-29-2011 01:04 AM
Base qualities during qseq2fastq and export2sam conversion hrajasim Illumina/Solexa 1 05-26-2011 12:20 AM
SAMTOOLS AND GATK input file (base qualities) hrajasim Illumina/Solexa 0 05-13-2011 06:58 AM
Converting ABI colorspace qualities into base qualities szilva Bioinformatics 7 04-02-2010 12:38 AM
Bowtie, Color-space reads, and confusing base qualities at variable sites keebs42 Bioinformatics 6 03-13-2010 05:51 PM

Reply
 
Thread Tools
Old 04-15-2012, 10:03 PM   #1
fuad193
Member
 
Location: N/A

Join Date: Feb 2012
Posts: 17
Question All Base Qualities Are Low

Hi again

sorry for posting again about the same topic. I will try to explain my problem in more details:

I have recently received whole genomic reads with Phred+33 quality from CompleteGenomics (FileFormat v1.5), but the problem is that all bases in all the reads files are range from 33 to 61 ASCII which mean that Phred is range from 0 to 28 only.

In another old exome seq data which comes from our partner's lab the ASCII ranges from 33 to 126 and I already analysed it with no problems.

1) Why the reads quality from CompleteGenomics looks different (i.e range only from 33 to 61 ASCII) ?

2) Do I need to rescale the quality before run it at BWA, Bowtie2, or whatever aligner ?

Thanks
fuad193 is offline   Reply With Quote
Old 04-16-2012, 12:32 AM   #2
tonybolger
Senior Member
 
Location: berlin

Join Date: Feb 2010
Posts: 156
Default

Quote:
Originally Posted by fuad193 View Post
In another old exome seq data which comes from our partner's lab the ASCII ranges from 33 to 126 and I already analysed it with no problems.
I would question this data first, to be honest. I don't know of any technology which gives reads with Q scores in the order of 90 - AFAIK that would be 1 error per billion bases. Assembled data can have Q scores that high, but not raw data.

Quote:
Originally Posted by fuad193 View Post
1) Why the reads quality from CompleteGenomics looks different (i.e range only from 33 to 61 ASCII) ?
It's likely correct - scores of Q20 - Q30 (error rates of around 1:100 to 1:1000) are normal for raw reads, if a little below the maximum you could hope for.

Quote:
Originally Posted by fuad193 View Post
2) Do I need to rescale the quality before run it at BWA, Bowtie2, or whatever aligner ?
No, the Q scores you have are reasonable and easily sufficient for alignment.
tonybolger 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 06:58 PM.


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