We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent 36cf937 commit a86d4a7Copy full SHA for a86d4a7
macrel/data/scripts/count.ngl
@@ -3,7 +3,7 @@ ngless "1.0"
3
input = samfile(ARGV[1], name='MACREL')
4
5
# Filter to keep only primary alignments (equivalent to -F 256)
6
-filtered = filter(input, keep_if="not secondary")
+filtered = select(input, keep_if=[!is_secondary])
7
8
# Count features from filtered alignments
9
counts = count(filtered,
0 commit comments