diff samtoolshelper.py @ 2:4d701c67592a

Uploaded
author nilesh
date Mon, 08 Jul 2013 22:03:50 -0400
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/samtoolshelper.py	Mon Jul 08 22:03:50 2013 -0400
@@ -0,0 +1,20 @@
+import sys
+import subprocess as sp
+import os
+
+# Creates the sorted and indexed bam/bai files that are requried for both bam2wig and RSEQC_count
+def samtools_sorted(bam):
+	sortedbam = bam + ".sorted"
+	indexedbam = ".".join([sortedbam,"bam.bai"])
+	sp.call(['samtools', 'sort', '-m 1000000000', bam, sortedbam])
+	sortedbam = sortedbam + '.bam'
+	sp.call(['samtools', 'index', sortedbam, indexedbam])
+	return sortedbam
+
+def main(args):
+	args[2] = samtools_sorted(args[2])
+	sp.call(args)
+
+
+if __name__ == "__main__":
+	main(sys.argv[1:])
\ No newline at end of file