File: phred.py

package info (click to toggle)
freebayes 1.3.7-1~exp
  • links: PTS, VCS
  • area: main
  • in suites: experimental
  • size: 6,324 kB
  • sloc: cpp: 125,947; ansic: 4,539; sh: 1,084; python: 604; asm: 271; lisp: 75; perl: 27; makefile: 20
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)