Mercurial > repos > bebatut > combine_metaphlan2_humann2
annotate combine_metaphlan2_humann2.py @ 1:e25efca0a49c draft
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
author | bebatut |
---|---|
date | Mon, 14 Sep 2020 12:19:49 +0000 |
parents | 31394a0c0242 |
children | fdfb35745104 |
rev | line source |
---|---|
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
1 #!/usr/bin/env python |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
2 # -*- coding: utf-8 -*- |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
3 |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
4 import argparse |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
5 |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
6 |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
7 def extract_clade_abundance(metaphlan2_fp): |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
8 clade_abund = {} |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
9 with open(metaphlan2_fp, 'r') as metaphlan2_f: |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
10 for line in metaphlan2_f.readlines(): |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
11 if line.find('g__') == -1: |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
12 continue |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
13 |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
14 split_line = line[:-1].split('\t') |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
15 taxo = split_line[0] |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
16 abundance = split_line[1] |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
17 |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
18 genus = taxo[(taxo.find('g__')+3):] |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
19 if genus.find('|') != -1: |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
20 genus = genus[:(genus.find('|'))] |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
21 clade_abund.setdefault(genus, {'abundance': 0, 'species': {}}) |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
22 if taxo.find('t__') != -1: |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
23 continue |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
24 elif taxo.find('s__') != -1: |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
25 species = taxo[(taxo.find('s__')+3):] |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
26 clade_abund[genus]['species'].setdefault( |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
27 species, |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
28 abundance) |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
29 else: |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
30 clade_abund[genus]['abundance'] = abundance |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
31 return clade_abund |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
32 |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
33 |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
34 def compute_overall_abundance(humann2_fp): |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
35 overall_abundance = 0 |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
36 with open(humann2_fp, 'r') as humann2_f: |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
37 for line in humann2_f.readlines(): |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
38 if line.find('|') != -1 or line.startswith('#'): |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
39 continue |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
40 split_line = line[:-1].split('\t') |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
41 overall_abundance += float(split_line[1]) |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
42 return overall_abundance |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
43 |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
44 |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
45 def format_characteristic_name(name): |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
46 formatted_n = name |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
47 formatted_n = formatted_n.replace('/', ' ') |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
48 formatted_n = formatted_n.replace('-', ' ') |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
49 formatted_n = formatted_n.replace("'", '') |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
50 if formatted_n.find('(') != -1 and formatted_n.find(')') != -1: |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
51 open_bracket = formatted_n.find('(') |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
52 close_bracket = formatted_n.find(')')+1 |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
53 formatted_n = formatted_n[:open_bracket] + formatted_n[close_bracket:] |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
54 return formatted_n |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
55 |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
56 |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
57 def combine_metaphlan2_humann2(args): |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
58 clade_abund = extract_clade_abundance(args.metaphlan2_fp) |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
59 overall_abund = compute_overall_abundance(args.humann2_fp) |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
60 |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
61 with open(args.output_fp, 'w') as output_f: |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
62 s = 'genus\tgenus_abundance\tspecies\tspecies_abundance\t' |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
63 s = '%s\t%s_id\t%s_name\t%s_abundance\n' % (s, args.type, args.type, args.type) |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
64 output_f.write(s) |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
65 with open(args.humann2_fp, 'r') as humann2_f: |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
66 for line in humann2_f.readlines(): |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
67 if line.find('|') == -1: |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
68 continue |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
69 |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
70 split_line = line[:-1].split('\t') |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
71 abundance = 100*float(split_line[1])/overall_abund |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
72 annotation = split_line[0].split('|') |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
73 charact = annotation[0].split(':') |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
74 charact_id = charact[0] |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
75 char_name = '' |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
76 if len(charact) > 1: |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
77 char_name = format_characteristic_name(charact[-1]) |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
78 taxo = annotation[1].split('.') |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
79 |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
80 if taxo[0] == 'unclassified': |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
81 continue |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
82 genus = taxo[0][3:] |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
83 species = taxo[1][3:] |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
84 |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
85 if genus not in clade_abund: |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
86 print("no %s found in %s" % (genus, args.metaphlan2_fp)) |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
87 continue |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
88 if species not in clade_abund[genus]['species']: |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
89 print("no %s found in %s for % s" % (species, args.metaphlan2_fp, genus)) |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
90 continue |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
91 |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
92 s = "%s\t%s\t" % (genus, clade_abund[genus]['abundance']) |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
93 s += "%s\t%s\t" % (species, clade_abund[genus]['species'][species]) |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
94 s += "%s\t%s\t%s\n" % (charact_id, char_name, abundance) |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
95 output_f.write(s) |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
96 |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
97 |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
98 if __name__ == '__main__': |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
99 parser = argparse.ArgumentParser() |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
100 parser.add_argument('--humann2_fp', required=True) |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
101 parser.add_argument('--metaphlan2_fp', required=True) |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
102 parser.add_argument('--output_fp', required=True) |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
103 parser.add_argument( |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
104 '--type', |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
105 required=True, |
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
106 choices=['gene_families', 'pathways']) |
0
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
107 args = parser.parse_args() |
31394a0c0242
planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit e6bee6545960c2a1ae3ca3031ec74d7c26d0b0ce-dirty
bebatut
parents:
diff
changeset
|
108 |
1
e25efca0a49c
"planemo upload for repository https://github.com/asaim/galaxytools/tree/master/tools/combine_metaphlan2_humann2 commit 4a40a73a6eb981310a3c1c54dfa98b5460a843f9"
bebatut
parents:
0
diff
changeset
|
109 combine_metaphlan2_humann2(args) |