This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Quality of whole genome sequencing and whole exome sequencing

I have some whole genome sequenced data with coverage ~6x and some whole exome sequenced data with coverage ~60x. The sequencing platforms are the same for the above data. However, I found there are more low phred scores in the whole genome sequenced data than that in the whole exome sequenced data. Could anyone tell me whether this is cased by the difference in coverages? Say deep sequencing leads to more high quality reads.

next-gen
• 3,589 views • link
updated by Ram Ram • written by devinliao0918 • devinliao0918

1 answer

Are you really talking about base-call phred scores or do you mean scores associated with variant calls? The former are unaffected by depth, the latter highly affected by it.

• 0 views • link
written by Devon Ryan Devon Ryan

I am talking about the base-call phred scores. All the scores are extracted from pileup files generated by "samtools pileup". I didn't call variants.

Also, the base-quality scores are sort of truncated at 41 because there is none phred score greater than 41 in the WGS data. The wiki page for FASTQ format says the following

For raw reads, the range of scores will depend on the technology and the base caller used, but will typically be up to 41 for recent Illumina chemistry. Since the maximum observed quality score was previously only 40, various scripts and tools break when they encounter data with quality values larger than 40. For processed reads, scores may be even higher. For example, quality values of 45 are observed in reads from Illumina's Long Read Sequencing Service (previously Moleculo).

• 0 views • link
updated by Ram Ram • written by devinliao0918 • devinliao0918

Phred scores are unlikely to be affected by WGS vs. exome seq.

• 0 views • link
written by Devon Ryan Devon Ryan

I agree with you. I guess the difference may be caused by different technologies, though I know both the WES and WGS are done on Illumina HiSeq 2000.

• 0 views • link
updated by Ram Ram • written by devinliao0918 • devinliao0918

I agree, this is more likely due to the library prep or difference in the source DNA.

• 0 views • link
updated by Ram Ram • written by User 59 User 59

Log in to answer this question.

More posts like this