simple seq to line script

jgurtowski · Jun 27, 2014 · e29ec82 · e29ec82
1 parent 5bae9e2
commit e29ec82
Show file tree

Hide file tree

Showing 2 changed files with 18 additions and 0 deletions.
diff --git a/sequenceToLine.py b/sequenceToLine.py
@@ -0,0 +1,17 @@
+#!/usr/bin/env
+
+import sys
+
+from seqio import iteratorFromExtension
+from nucio import fileIterator
+
+if not len(sys.argv) == 2:
+    sys.exit("sequencToLine.py in.{fa.fq}\n")
+
+it = iteratorFromExtension(sys.argv[1])
+for record in fileIterator(sys.argv[1], it):
+    if hasattr(record, "desc"):
+        print "\t".join([record.name, record.seq, record.desc, record.qual])
+    else:
+        print "\t".join([record.name, record.seq])
+
diff --git a/zmwProductivityHeatmap.py b/zmwProductivityHeatmap.py
@@ -44,6 +44,7 @@
 plt.suptitle(cellname)
 
 plt.savefig(pp, format="pdf")
+plt.legend()
 
 plt.figure()