54 lines
1.7 KiB
Python
54 lines
1.7 KiB
Python
|
import farm_commands
|
||
|
|
import os.path
|
||
|
|
import sys
|
||
|
|
import getopt
|
||
|
|
|
||
|
|
defaultCommands = ['CountReads', 'Pileup']
|
||
|
|
|
||
|
|
def usage():
|
||
|
|
print "Optional arguments:"
|
||
|
|
print " -f QUEUE Farm jobs to QUEUE on LSF"
|
||
|
|
print " -c cmd1,cmd2 Walkers to execute, otherwise", ' '.join(defaultCommands)
|
||
|
|
|
||
|
|
if __name__ == "__main__":
|
||
|
|
opts = None
|
||
|
|
try:
|
||
|
|
opts, args = getopt.getopt(sys.argv[1:], "f:c:i", ["farm", "commands", "ignoreExistingFiles"])
|
||
|
|
except getopt.GetoptError:
|
||
|
|
print sys.argv
|
||
|
|
usage()
|
||
|
|
sys.exit(2)
|
||
|
|
|
||
|
|
farm_sub = False
|
||
|
|
commandsList = defaultCommands
|
||
|
|
ignoreExistingFiles = False
|
||
|
|
|
||
|
|
for opt, arg in opts:
|
||
|
|
if opt in ("-f", "--farm"):
|
||
|
|
farm_sub = arg
|
||
|
|
if opt in ("-c", "--commands"):
|
||
|
|
commandsList = arg.split(',')
|
||
|
|
if opt in ("-e", "--ignoreExistingFiles"):
|
||
|
|
ignoreExistingFiles = True
|
||
|
|
|
||
|
|
directory = args[1]
|
||
|
|
for line in open(args[0]):
|
||
|
|
lane = line.strip()
|
||
|
|
|
||
|
|
if not os.path.exists(lane):
|
||
|
|
print 'Input SAM/BAM file: "', lane, '" does not exist, skipping...'
|
||
|
|
continue
|
||
|
|
|
||
|
|
head, lane_filename = os.path.split(lane)
|
||
|
|
filebase = os.path.splitext(lane_filename)[0]
|
||
|
|
|
||
|
|
# convert the fasta
|
||
|
|
for analysis in commandsList:
|
||
|
|
output = os.path.join(directory, filebase + '.' + analysis + '.output')
|
||
|
|
if ignoreExistingFiles or not os.path.exists(output):
|
||
|
|
cmd = "java -jar ~/dev/GenomeAnalysisTK/trunk/dist/GenomeAnalysisTK.jar T=" + analysis + " I= " + lane + " R= /seq/references/Homo_sapiens_assembly18/v0/Homo_sapiens_assembly18.fasta "
|
||
|
|
print cmd
|
||
|
|
farm_commands.cmd(cmd, farm_sub, output)
|
||
|
|
|
||
|
|
|