File: phred.py

package info (click to toggle)
freebayes 1.3.9-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 6,984 kB
  • sloc: cpp: 125,778; ansic: 4,581; sh: 1,084; python: 672; asm: 271; javascript: 94; lisp: 85; makefile: 37; perl: 27
file content (17 lines) | stat: -rw-r--r-- 299 bytes parent folder | download | duplicates (5) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
import math

M_LN10 = math.log(10)
M_LOG10E = math.log10(math.e)

def phred2ln(qual):
    return M_LN10 * qual * -.1

def ln2phred(prob):
    return -10 * M_LOG10E * prob

def phred2float(qual):
    return math.pow(10, qual * -.1)

def float2phred(prob):
    return min(-10 * math.log10(prob), 99)