diff overlapping_reads.py @ 3:c0e72bd28453 draft

planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_signatures commit bc672cf9d5fe9c1f7eaf02abdc906cf2c1763668

--- a/overlapping_reads.py	Sat Oct 22 23:49:14 2022 +0000
+++ b/overlapping_reads.py	Sat Feb 25 09:54:18 2023 +0000
@@ -90,12 +90,10 @@
         target_range = self.target_range
         uppers = [seq for seq in uppers if (len(seq) in query_range or len(seq)
                   in target_range)]
-        print(uppers)
         uppers_expanded = []
         for seq in uppers:
             expand = [seq for i in range(self.readdic[seq])]
             uppers_expanded.extend(expand)
-        print(uppers_expanded)
         uppers = uppers_expanded
         lowers = [seq for seq in lowers if (len(seq) in query_range or len(seq)
                   in target_range)]
@@ -123,7 +121,6 @@
         stringresult = []
         header_template = '>%s|coord=%s|strand %s|size=%s|nreads=%s\n%s\n'
         total_pairs = 0
-        print('Chromosome\tNbre of pairs')
         for chrom in sorted(self.chromosomes):
             number_pairs = 0
             for pos in self.all_query_positions[chrom]:
@@ -149,9 +146,6 @@
                                      len(downread), self.readdic[downread],
                                      self.revcomp(downread)))
                 stringresult.extend(sorted(set(stringbuffer)))
-            print('%s\t%s' % (chrom, number_pairs))
-        print('Total nbre of pairs that can be simultaneously formed\t%s'
-              % total_pairs)
         F.write(''.join(stringresult))
 
     def revcomp(self, sequence):