annotate fasta_tabular_converter.py @ 4:87c99e4af616 draft default tip

planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
author drosofff
date Wed, 09 Nov 2016 11:24:50 -0500
parents 403f0769fc1c
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
4
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
1 #!/usr/bin/env python
0
70f2654ad7e2 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
2 #
4
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
3 import argparse
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
4 import logging
0
70f2654ad7e2 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
5 import sys
70f2654ad7e2 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
6 from collections import defaultdict
70f2654ad7e2 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
7
4
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
8
2
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
9 def Parser():
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
10 the_parser = argparse.ArgumentParser()
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
11 the_parser.add_argument(
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
12 '--input', action="store", type=str, help="input file")
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
13 the_parser.add_argument(
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
14 '--output', action="store", type=str, help="output converted file")
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
15 the_parser.add_argument(
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
16 '--type', action="store", type=str, help="type of convertion")
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
17 args = the_parser.parse_args()
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
18 return args
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
19
4
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
20
2
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
21 def readfasta_writetabular(fasta, tabular, mode="oneline"):
4
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
22 for line in fasta:
2
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
23 if line[0] == ">":
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
24 try:
4
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
25 seqdic["".join(stringlist)] += 1 # to dump the sequence of the previous item - try because of first missing stringlist variable
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
26 except NameError:
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
27 pass
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
28 stringlist = []
2
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
29 else:
4
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
30 try:
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
31 stringlist.append(line[:-1])
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
32 except UnboundLocalError: # if file went through filter and contains only empty lines
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
33 logging.info("first line is empty.")
3
403f0769fc1c planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 030207144f0811822dbdda9a10e036ff8e794d7c
drosofff
parents: 2
diff changeset
34 try:
4
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
35 seqdic["".join(stringlist)] += 1 # for the last sequence
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
36 except NameError:
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
37 logging.info("input file has not fasta sequences.")
2
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
38 for seq in sorted(seqdic, key=seqdic.get, reverse=True):
4
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
39 tabular.write("%s\t%s\n" % (seq, seqdic[seq]))
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
40
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
41
0
70f2654ad7e2 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
42 def readtabular_writefasta(tabular, fasta):
4
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
43 counter = 0
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
44 for line in tabular:
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
45 fields = line.split()
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
46 for i in range(int(fields[1])):
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
47 counter += 1
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
48 fasta.write(">%s\n%s\n" % (counter, fields[0]))
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
49
0
70f2654ad7e2 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
50
4
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
51 def readtabular_writefastaweighted(tabular, fasta):
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
52 counter = 0
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
53 for line in tabular:
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
54 counter += 1
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
55 fields = line[:-1].split()
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
56 fasta.write(">%s_%s\n%s\n" % (counter, fields[1], fields[0]))
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
57
0
70f2654ad7e2 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
58
4
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
59 def readfastaweighted_writefastaweighted(fastaweigthed_input, fastaweigthed_reparsed):
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
60 number_reads = 0
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
61 for line in fastaweigthed_input:
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
62 if line[0] == ">":
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
63 weigth = int(line[1:-1].split("_")[-1])
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
64 number_reads += weigth
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
65 else:
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
66 seqdic[line[:-1]] += weigth
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
67 n = 0
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
68 for seq in sorted(seqdic, key=seqdic.get, reverse=True):
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
69 n += 1
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
70 fastaweigthed_reparsed.write(">%s_%s\n%s\n" % (n, seqdic[seq], seq))
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
71 log.info("%s reads collapsed" % number_reads)
0
70f2654ad7e2 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
72
4
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
73
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
74 def readfastaweighted_writefasta(fastaweigthed, fasta):
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
75 counter = 0
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
76 for line in fastaweigthed:
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
77 if line[0] == ">":
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
78 weigth = int(line[1:-1].split("_")[-1])
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
79 else:
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
80 seq = line[:-1]
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
81 for i in range(weigth):
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
82 counter += 1
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
83 fasta.write(">%s\n%s\n" % (counter, seq))
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
84
0
70f2654ad7e2 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
85
2
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
86 def main(input, output, type):
4
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
87 with open(input, "r") as input:
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
88 with open(output, "w") as output:
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
89 if type == "fasta2tabular":
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
90 readfasta_writetabular(input, output)
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
91 elif type == "tabular2fasta":
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
92 readtabular_writefasta(input, output)
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
93 elif type == "tabular2fastaweight":
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
94 readtabular_writefastaweighted(input, output)
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
95 elif type == "fastaweight2fastaweight":
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
96 readfastaweighted_writefastaweighted(input, output)
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
97 elif type == "fastaweight2fasta":
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
98 readfastaweighted_writefasta(input, output)
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
99
0
70f2654ad7e2 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
100
2
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
101 if __name__ == "__main__":
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
102 seqdic = defaultdict(int)
87e27aa012d7 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit 6a93f2809e2939f9d847c3238bfbff8ead367d9f
drosofff
parents: 0
diff changeset
103 args = Parser()
4
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
104 log = logging.getLogger(__name__)
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
105 logging.basicConfig(stream=sys.stdout, level=logging.INFO)
87c99e4af616 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_fasta_tabular_converter commit b6de14061c479f0418cd89e26d6f5ac26e565a07
drosofff
parents: 3
diff changeset
106 main(args.input, args.output, args.type)