Skip to content

Commit 702a388

Browse files
Update count.ngl
1 parent a8d30b3 commit 702a388

File tree

1 file changed

+5
-1
lines changed

1 file changed

+5
-1
lines changed

macrel/data/scripts/count.ngl

Lines changed: 5 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,9 +1,13 @@
11
ngless "1.0"
22

3-
input = samfile(ARGV[1], name='MACREL')
3+
input = samfile(ARGV[1],
4+
remove_flags=[256],
5+
name="MACREL")
6+
47
counts = count(input,
58
features=["seqname"],
69
include_minus1=False,
710
normalization={raw})
811

12+
counts = as_integer(round(counts))
913
write(counts, ofile=ARGV[2])

0 commit comments

Comments
 (0)