Mercurial > repos > iuc > ivar_trim
changeset 11:5671e1d3d5ee 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:18:33 +0000 |
parents | 9f978da6528a |
children | e3bb03cf207a |
files | ivar_trim.xml prepare_amplicon_info.py sanitize_bed.py write_amplicon_info_file.py |
diffstat | 4 files changed, 18 insertions(+), 8 deletions(-) [+] |
line wrap: on
line diff
--- a/ivar_trim.xml Sun Mar 20 11:43:34 2022 +0000 +++ b/ivar_trim.xml Wed Jul 13 15:18:33 2022 +0000 @@ -1,4 +1,4 @@ -<tool id="ivar_trim" name="ivar trim" version="@VERSION@+galaxy3"> +<tool id="ivar_trim" name="ivar trim" version="@VERSION@+galaxy4" profile="21.01"> <description>Trim reads in aligned BAM</description> <macros> <import>macros.xml</import>
--- a/prepare_amplicon_info.py Sun Mar 20 11:43:34 2022 +0000 +++ b/prepare_amplicon_info.py Wed Jul 13 15:18:33 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 Sun Mar 20 11:43:34 2022 +0000 +++ b/sanitize_bed.py Wed Jul 13 15:18:33 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 Sun Mar 20 11:43:34 2022 +0000 +++ b/write_amplicon_info_file.py Wed Jul 13 15:18:33 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)