Mercurial > repos > artbio > repenrich2

diff RepEnrich2.py @ 5:08e50af788f7 draft
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/main/tools/repenrich2 commit b1761de76fd068b86a06d88e70c1ba1d8644e7b5
author: artbio
date: Sun, 21 Apr 2024 21:52:40 +0000
parents: 4905a332a094
children: 388a47ca4199
--- a/RepEnrich2.py	Sat Apr 20 23:23:40 2024 +0000
+++ b/RepEnrich2.py	Sun Apr 21 21:52:40 2024 +0000
@@ -115,6 +115,12 @@
     sumofrepeatreads += int(line[4])
 print(f"Identified {sumofrepeatreads} unique reads that mapped to repeats.")
 
+# print unique mapper counts
+with open("unique_mapper_counts.tsv", 'w') as fout:
+    fout.write("#element\tcount\n")
+    for count in sorted(counts):
+        fout.write(f"{count}\t{counts[count]}\n")
+
 # multimapper parsing
 if not paired_end:
     args_list = [(metagenome, fastqfile_1) for