Repository 'medaka_variant_pipeline'
hg clone https://toolshed.g2.bx.psu.edu/repos/iuc/medaka_variant_pipeline

Changeset 7:08e0d74aac22 (2021-02-22)
Previous changeset 6:ea1833858055 (2021-02-10) Next changeset 8:1040a307be0c (2021-02-23)
Commit message:
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/medaka commit ed5a3aadbecc0decf9a797447f3ac7700683ea9a"
modified:
convert_VCF_info_fields.py
b
diff -r ea1833858055 -r 08e0d74aac22 convert_VCF_info_fields.py
--- a/convert_VCF_info_fields.py Wed Feb 10 08:29:07 2021 +0000
+++ b/convert_VCF_info_fields.py Mon Feb 22 19:58:33 2021 +0000
[
@@ -15,7 +15,11 @@
 
 
 def pval_to_phredqual(pval):
-    return round(-10 * log10(pval))
+    try:
+        ret = round(-10 * log10(pval))
+    except ValueError:
+        ret = 2147483647  # transform pval of 0.0 to max signed 32 bit int
+    return ret
 
 
 def parseInfoField(info):
@@ -33,7 +37,7 @@
     to_skip = set(['SC', 'SR'])
     for i, line in enumerate(in_vcf):
         if i == 1:
-            out_vcf.write("##convert_VCF_info_fields=0.1\n")
+            out_vcf.write("##convert_VCF_info_fields=0.2\n")
         if line[0:2] == "##":
             if line[0:11] == "##INFO=<ID=":
                 id_ = line[11:].split(',')[0]