| Previous changeset 6:d915879a88bc (2021-02-10) Next changeset 8:2c72b8c09a70 (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 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=<ID=": id_ = line[11:].split(',')[0] |