work on filters

1 files changed, 17 insertions, 0 deletions

diff --git a/src/remove_mers.py b/src/remove_mers.py
new file mode 100644
index 0000000..b18992a
--- /dev/null
+++ b/src/remove_mers.py
@@ -0,0 +1,17 @@
+#!/usr/bin/env python2.7
+import sys
+
+if __name__ == "__main__":
+
+	if len(sys.argv) is 1:
+		exit("Filter mers, input is stdin, output is stdout, mers are argv")
+	
+	mers_to_delete = set()
+
+	for mer in sys.argv[1:]:
+		mers_to_delete.add(mer)
+
+	for line in sys.stdin:
+		if line.split()[0] not in mers_to_delete:
+			sys.stdout.write(line)
+