Mercurial > repos > iuc > medaka_consensus

comparison convert_VCF_info_fields.py @ 11:4e65663b2567 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/medaka commit 0faf0ade3f13d7c78d93869823ea9fdf25c21b13"
author iuc
date Thu, 22 Apr 2021 20:21:17 +0000
parents a7543e0dccc0
children a5b6aed7e6d4
comparison
equal deleted inserted replaced
10:b8d86e7780e8 11:4e65663b2567
9 9
10 import sys 10 import sys
11 from collections import OrderedDict 11 from collections import OrderedDict
12 from math import log10 12 from math import log10
13 13
14 from scipy.stats import fisher_exact 14 import scipy
15 import scipy.stats
15 16
16 17
17 def pval_to_phredqual(pval): 18 def pval_to_phredqual(pval):
18 try: 19 try:
19 ret = round(-10 * log10(pval)) 20 ret = round(-10 * log10(pval))
67 dpspf += sr_list[i] 68 dpspf += sr_list[i]
68 dpspr += sr_list[i + 1] 69 dpspr += sr_list[i + 1]
69 for j, i in enumerate(range(2, len(sr_list), 2)): 70 for j, i in enumerate(range(2, len(sr_list), 2)):
70 dp4 = (sr_list[ref_fwd], sr_list[ref_rev], sr_list[i], sr_list[i + 1]) 71 dp4 = (sr_list[ref_fwd], sr_list[ref_rev], sr_list[i], sr_list[i + 1])
71 dp2x2 = [[dp4[0], dp4[1]], [dp4[2], dp4[3]]] 72 dp2x2 = [[dp4[0], dp4[1]], [dp4[2], dp4[3]]]
72 _, p_val = fisher_exact(dp2x2) 73 _, p_val = scipy.stats.fisher_exact(dp2x2)
73 sb = pval_to_phredqual(p_val) 74 sb = pval_to_phredqual(p_val)
74 75
75 as_ = (sc_list[ref_fwd], sc_list[ref_rev], sc_list[i], sc_list[i + 1]) 76 as_ = (sc_list[ref_fwd], sc_list[ref_rev], sc_list[i], sc_list[i + 1])
76 77
77 info = [] 78 info = []