# HG changeset patch # User iuc # Date 1614023982 0 # Node ID 6ca28e6144fab8bfcf3c9bb9902b981b20f547a7 # Parent fbf007713188c3041c54d58acc935a081264c235 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/medaka commit ed5a3aadbecc0decf9a797447f3ac7700683ea9a" diff -r fbf007713188 -r 6ca28e6144fa convert_VCF_info_fields.py --- a/convert_VCF_info_fields.py Wed Feb 10 08:28:25 2021 +0000 +++ b/convert_VCF_info_fields.py Mon Feb 22 19:59:42 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=