Mercurial > repos > artbio > mircounts

comparison yac.py @ 15:ffcd42f85b61 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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):