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") |