# HG changeset patch # User iuc # Date 1614023913 0 # Node ID 08e0d74aac22e7b98c73a9d2414599e5219f8f7e # Parent ea1833858055d02ad6c30e629e14060037296142 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/medaka commit ed5a3aadbecc0decf9a797447f3ac7700683ea9a" 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=