Mercurial > repos > nilesh > rseqc
changeset 2:4d701c67592a
Uploaded
author | nilesh |
---|---|
date | Mon, 08 Jul 2013 22:03:50 -0400 (2013-07-09) |
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