Skip to content

Instantly share code, notes, and snippets.

@dreidpath

dreidpath/phred.R

Created September 30, 2018 12:45
Show Gist options
  • Save dreidpath/b0b66aad020a3c9f3ff5739836d39868 to your computer and use it in GitHub Desktop.
Save dreidpath/b0b66aad020a3c9f3ff5739836d39868 to your computer and use it in GitHub Desktop.
Download ZIP
The R functions for the caluclation of the phred quality score, " a measure of the quality of the identification of the nucleobases generated by automated DNA sequencing" https://en.wikipedia.org/wiki/Phred_quality_score
Raw
phred.R
phred <- function(q){
10^(-q/10)
}
inv_phred <- function(p){
10 * log10(p)
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment