Mercurial > repos > artbio > mircounts
comparison yac.py @ 15:ffcd42f85b61 draft default tip
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/mircounts commit 5eb8570dce4e22fb2759cc16c8e1ce9d304508fe
author | artbio |
---|---|
date | Sat, 10 Feb 2024 17:15:04 +0000 |
parents | 6b8adacd4750 |
children |
comparison
equal
deleted
inserted
replaced
14:c163574c246f | 15:ffcd42f85b61 |
---|---|
79 else: | 79 else: |
80 block = "@HWI-{0}\n{1}\n+\n{2}\n".format(id, read, qscore) | 80 block = "@HWI-{0}\n{1}\n+\n{2}\n".format(id, read, qscore) |
81 output.write(block) | 81 output.write(block) |
82 | 82 |
83 def handle_io(self): | 83 def handle_io(self): |
84 '''Open input file, pass read sequence and read qscore to clipping function. | 84 ''' |
85 Pass clipped read and qscore to output function.''' | 85 Open input file, pass read sequence and read qscore |
86 to clipping function. | |
87 Pass clipped read and qscore to output function. | |
88 ''' | |
86 id = 0 | 89 id = 0 |
87 output = open(self.outputfile, "a") | 90 output = open(self.outputfile, "a") |
88 with open(self.inputfile, "r") as input: | 91 with open(self.inputfile, "r") as input: |
89 block_gen = islice(input, 1, None, 2) | 92 block_gen = islice(input, 1, None, 2) |
90 for i, line in enumerate(block_gen): | 93 for i, line in enumerate(block_gen): |