We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent 702a388 commit bae6a99Copy full SHA for bae6a99
macrel/data/scripts/count.ngl
@@ -1,9 +1,10 @@
1
ngless "1.0"
2
3
input = samfile(ARGV[1],
4
- remove_flags=[256],
5
name="MACREL")
6
+filtered = select(input, keep_if=[{unique}])
7
+
8
counts = count(input,
9
features=["seqname"],
10
include_minus1=False,
0 commit comments