-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathgenerate.py
executable file
·70 lines (69 loc) · 1.97 KB
/
generate.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
#!/usr/bin/python
import sys
import math
import os
import numpy as np
import scipy
import json
from collections import OrderedDict
if __name__=="__main__":
config = OrderedDict([
("Rscriptbin", "/n/app/R/3.3.3/bin"),
("pythonbin", "/n/app/python/2.7.12/bin"),
("trimmomaticbin", "/n/app/trimmomatic/0.36/bin"),
("trimmomaticjarfile", "trimmomatic-0.36.jar"),
("bowtie2bin", "/n/app/bowtie2/2.2.9/bin"),
("samtoolsbin", "/n/app/samtools/1.3.1/bin"),
("adapterpath", "/home/qz64"),
("picardbin", "/n/app/picard/2.8.0/bin"),
("picardjarfile", "picard-2.8.0.jar"),
("macs2bin", "/n/app/macs2/2.1.1.20160309/bin"),
("kseqbin", "/home/qz64"),
("memebin", "/home/qz64/meme/bin"),
("bedopsbin", "/home/qz64/bin"),
("bedtoolsbin", "/n/app/bedtools/2.26.0/bin"),
("makecutmatrixbin", "/home/qz64/.local/bin"),
("bt2idx", "/n/groups/shared_databases/bowtie2_indexes"),
("genome_sequence", "/home/qz64/chrom.hg19/hg19.fa"),
("extratoolsbin", "/home/qz64"),
("extrasettings", "/home/qz64"),
("input/output", OrderedDict([
("fastq_directory", "/n/scratch2/qz64/Nan_18_aug23/Nan_run_19"),
("workdir", "/n/scratch2/qz64/workdir"),
("fastq_sequence_length", 40),
("organism_build", "hg19"),
])),
("motif_finding", OrderedDict([
("num_bp_from_summit", 150),
("num_peaks", 5000),
("total_peaks", 15000),
("motif_scanning_pval", 0.0005),
("num_motifs", 20),
])),
("cluster", OrderedDict([
("email", "bernardzhu@gmail.com"),
("step_alignment", OrderedDict([
("queue", "short"),
("memory", 32000),
("time_limit", "0-12:00"),
])),
("step_process_bam", OrderedDict([
("queue", "short"),
("memory", 32000),
("time_limit", "0-12:00"),
])),
("step_motif_find", OrderedDict([
("queue", "short"),
("memory", 32000),
("time_limit", "0-12:00"),
])),
("step_footprinting", OrderedDict([
("queue", "short"),
("memory", 32000),
("time_limit", "0-12:00"),
])),
])),
])
fw = open("config.json", "w")
json.dump(config, fw, indent=4)
fw.close()