Mercurial > repos > nilesh > rseqc_dev2
view samtoolshelper.py @ 2:74e66b89b98b
Deleted selected files
author | nilesh |
---|---|
date | Mon, 29 Jul 2013 12:24:22 -0400 |
parents | dc3b3b88fbab |
children |
line wrap: on
line source
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:])