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