Mercurial > repos > galaxyp > metagene_annotator
annotate convert_mga.py @ 1:17c7ab82bfbc draft default tip
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit f66aa9a65aaf45ac4888ca26f5af661f73b46fdc
author | galaxyp |
---|---|
date | Mon, 18 Mar 2024 12:34:58 +0000 |
parents | b04960a7abf5 |
children |
rev | line source |
---|---|
0
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
1 #!/usr/bin/env python |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
2 |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
3 from __future__ import print_function |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
4 |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
5 import argparse |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
6 import math |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
7 import re |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
8 import sys |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
9 |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
10 |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
11 def __main__(): |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
12 parser = argparse.ArgumentParser( |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
13 description='Convert mga output to bed and tsv') |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
14 parser.add_argument( |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
15 'input_mga', |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
16 help="mga output to convert, '-' for stdin") |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
17 parser.add_argument( |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
18 '-t', '--tsv', default=None, |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
19 help='Path to output mga.tsv') |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
20 parser.add_argument( |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
21 '-b', '--bed', default=None, |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
22 help='Path to output mga.bed') |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
23 parser.add_argument('-v', '--verbose', action='store_true', help='Verbose') |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
24 args = parser.parse_args() |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
25 |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
26 input_rdr = open(args.input_mga, 'r')\ |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
27 if args.input_mga != '-' else sys.stdin |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
28 |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
29 bed_wtr = open(args.bed, 'w') if args.bed is not None else None |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
30 tsv_wtr = open(args.tsv, 'w') if args.bed is not None else None |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
31 if tsv_wtr: |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
32 tsv_wtr.write('#%s\n' % '\t'.join([ |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
33 'seq_id', 'seq_model', 'seq_gc', 'seq_rbs', |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
34 'gene ID', 'start pos', 'end pos', 'strand', 'frame', |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
35 'complete/partial', 'gene score', 'used model', |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
36 'rbs start', 'rbs end', 'rbs score'])) |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
37 |
1
17c7ab82bfbc
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit f66aa9a65aaf45ac4888ca26f5af661f73b46fdc
galaxyp
parents:
0
diff
changeset
|
38 gc_rbs_pat = '# gc = (-?[0-9]*[.]?[0-9]+), rbs = (-?[0-9]*[.]?[0-9]+)' |
0
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
39 seq_count = 0 |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
40 gene_count = 0 |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
41 for i, line in enumerate(input_rdr): |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
42 # 1317/1 |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
43 # gc = 0.272955, rbs = -1 |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
44 # self: - |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
45 if line.startswith('# gc'): |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
46 try: |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
47 m = re.match(gc_rbs_pat, line.strip()) |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
48 seq_gc, seq_rbs = m.groups() |
1
17c7ab82bfbc
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit f66aa9a65aaf45ac4888ca26f5af661f73b46fdc
galaxyp
parents:
0
diff
changeset
|
49 except ValueError: |
0
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
50 seq_gc = seq_rbs = '' |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
51 elif line.startswith('# self:'): |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
52 seq_type = re.sub('# self:', '', line.rstrip()) |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
53 elif line.startswith('# '): |
1
17c7ab82bfbc
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit f66aa9a65aaf45ac4888ca26f5af661f73b46fdc
galaxyp
parents:
0
diff
changeset
|
54 seq_name = re.sub(r'# (\S+).*$', '\\1', line.rstrip()) |
0
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
55 seq_count += 1 |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
56 else: |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
57 fields = line.split('\t') |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
58 if len(fields) == 11: |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
59 gene_count += 1 |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
60 start = int(fields[1]) - 1 |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
61 end = int(fields[2]) |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
62 if tsv_wtr: |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
63 tsv_wtr.write('%s\t%s\t%s\t%s\t%s' % ( |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
64 seq_name, |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
65 seq_type, |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
66 seq_gc, |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
67 seq_rbs, |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
68 line)) |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
69 if bed_wtr: |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
70 bed_wtr.write( |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
71 '%s\t%d\t%d\t%s\t%s\t%s\t%d\t%d\t%s\t%s\t%s\t%s\n' % ( |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
72 seq_name, |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
73 start, |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
74 end, |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
75 '%s:%s' % (seq_name, fields[0]), |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
76 int(math.ceil(float(fields[6]))), |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
77 fields[3], |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
78 start, |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
79 end, |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
80 0, |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
81 1, |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
82 abs(end - start), |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
83 0)) |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
84 |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
85 if args.verbose: |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
86 print("sequences: %d\tgenes: %d" |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
87 % (seq_count, gene_count), file=sys.stdout) |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
88 |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
89 |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
90 if __name__ == "__main__": |
b04960a7abf5
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/metagene_annotator commit 6d8b6e0fa2f1b47b337dbf21f5bc320586ccbd4c
galaxyp
parents:
diff
changeset
|
91 __main__() |