Mercurial > repos > earlhaminst > t_coffee
comparison filter_by_fasta_ids.py @ 1:b3833e5b50d4 draft
planemo upload for repository https://github.com/TGAC/earlham-galaxytools/tree/master/tools/t_coffee commit 011cabb2a2b3237bbbc4850ed26972816702a2ba-dirty
author | earlhaminst |
---|---|
date | Mon, 19 Dec 2016 17:47:31 -0500 |
parents | 794a6e864a96 |
children | 78dd29aa7fc1 |
comparison
equal
deleted
inserted
replaced
0:794a6e864a96 | 1:b3833e5b50d4 |
---|---|
80 | 80 |
81 used_sequences = set() | 81 used_sequences = set() |
82 work_summary = {'wanted': 0, 'found': 0, 'duplicates': 0} | 82 work_summary = {'wanted': 0, 'found': 0, 'duplicates': 0} |
83 targets = [] | 83 targets = [] |
84 | 84 |
85 f_target = open(sys.argv[1]) | 85 with open(sys.argv[1]) as f_target: |
86 for line in f_target.readlines(): | 86 for line in f_target.readlines(): |
87 targets.append(">%s" % line.strip().upper()) | 87 targets.append(">%s" % line.strip().upper()) |
88 f_target.close() | |
89 | 88 |
90 work_summary['wanted'] = len(targets) | 89 work_summary['wanted'] = len(targets) |
91 homd_db = FASTAReader(sys.argv[2]) | 90 homd_db = FASTAReader(sys.argv[2]) |
92 | 91 |
93 # output = open(sys.argv[3], "w") | 92 # output = open(sys.argv[3], "w") |