Mercurial > repos > nilesh > rseqc

changeset 2:4d701c67592a

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author nilesh
date Mon, 08 Jul 2013 22:03:50 -0400
parents d9ac0d9cb48e
children 203bda4cf122
files samtoolshelper.py
diffstat 1 files changed, 20 insertions(+), 0 deletions(-) [+]
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