Mercurial > repos > iuc > ivar_filtervariants
changeset 10:0135f1d487a0 draft
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/ivar/ commit 8ce6fd9aee543d9e62db33a9c95f79d8dc4e6dea
| author | iuc |
|---|---|
| date | Wed, 13 Jul 2022 15:19:34 +0000 |
| parents | 47a27906d750 |
| children | 72ddd32392cd |
| files | prepare_amplicon_info.py sanitize_bed.py write_amplicon_info_file.py |
| diffstat | 3 files changed, 17 insertions(+), 7 deletions(-) [+] |
line wrap: on
line diff
--- a/prepare_amplicon_info.py Fri Aug 20 20:32:33 2021 +0000 +++ b/prepare_amplicon_info.py Wed Jul 13 15:19:34 2022 +0000 @@ -11,7 +11,10 @@ primer_starts = {} with open(sys.argv[1]) as i: for line in i: - f = line.strip().split('\t') + line = line.strip() + if not line: + continue + f = line.split('\t') try: if f[5] == '+': primer_starts[f[3]] = int(f[1]) @@ -32,8 +35,11 @@ with open(sys.argv[2]) as i: ret_lines = [] for line in i: + line = line.strip() + if not line: + continue first = last = None - for pname in line.strip().split('\t'): + for pname in line.split('\t'): try: primer_start = primer_starts[pname] except KeyError:
--- a/sanitize_bed.py Fri Aug 20 20:32:33 2021 +0000 +++ b/sanitize_bed.py Wed Jul 13 15:19:34 2022 +0000 @@ -9,10 +9,11 @@ sanitized_data = [] try: for record in bed_data: - fields = record.split('\t') - sanitized_data.append( - '\t'.join(fields[:4] + ['60'] + fields[5:]) - ) + if record.strip(): + fields = record.split('\t') + sanitized_data.append( + '\t'.join(fields[:4] + ['60'] + fields[5:]) + ) except IndexError: pass # leave column number issue to getmasked else:
--- a/write_amplicon_info_file.py Fri Aug 20 20:32:33 2021 +0000 +++ b/write_amplicon_info_file.py Wed Jul 13 15:19:34 2022 +0000 @@ -10,7 +10,10 @@ def write_amplicon_info_file(bed_file, amplicon_info_file): amplicon_sets = {} for line in bed_file: - fields = line.strip().split('\t') + line = line.strip() + if not line: + continue + fields = line.split('\t') start = int(fields[1]) name = fields[3] re_match = AMPLICON_PAT.match(name)