# HG changeset patch # User iuc # Date 1614023961 0 # Node ID 57beb8ea55b4214ae8ad43ac087813f5e3aae297 # Parent d915879a88bc8ec8cd0b7891fa9a4d70010cda51 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/medaka commit ed5a3aadbecc0decf9a797447f3ac7700683ea9a" diff -r d915879a88bc -r 57beb8ea55b4 convert_VCF_info_fields.py --- a/convert_VCF_info_fields.py Wed Feb 10 08:28:46 2021 +0000 +++ b/convert_VCF_info_fields.py Mon Feb 22 19:59:21 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=