-
Notifications
You must be signed in to change notification settings - Fork 40
/
Copy pathexample2.sh
125 lines (80 loc) · 6.22 KB
/
example2.sh
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
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
screen -RD HAIB_GATA2_K562_SE
DATA_ROOT=/srv/scratch/leepc12/data; WORK_DIR=/srv/scratch/leepc12/run/chipseq_test; SUFFIX=GATA2/HAIB_GATA2_K562_SE
mkdir -p ${WORK_DIR}/$SUFFIX; cd ${WORK_DIR}/$SUFFIX;
bds /users/leepc12/code/bds_atac/chipseq/chipseq.bds \
-fastq1 ${DATA_ROOT}/GATA2/SE/HAIB_GATA2_K562_SE/HAIB_GATA2_K562_SE_Rep1.fastq.gz \
-fastq2 ${DATA_ROOT}/GATA2/SE/HAIB_GATA2_K562_SE/HAIB_GATA2_K562_SE_Rep2.fastq.gz \
-ctl_fastq1 ${DATA_ROOT}/GATA2/SE/HAIB_GATA2_K562_SE/HAIB_Input_K562_SE_Rep1.fastq.gz \
-ctl_fastq2 ${DATA_ROOT}/GATA2/SE/HAIB_GATA2_K562_SE/HAIB_Input_K562_SE_Rep2.fastq.gz \
-url_base http://mitra.stanford.edu/kundaje/leepc12/chipseq_test/${SUFFIX}/out \
-species hg19 -bwa_idx ${DATA_ROOT}/CTCF/bwa_index/male.hg19.fa.gz -chrsz ${DATA_ROOT}/CTCF/male.hg19.chrom.sizes
screen -RD Snyder_CTCF_GM12878_PE
DATA_ROOT=/srv/scratch/leepc12/data; WORK_DIR=/srv/scratch/leepc12/run/chipseq_test; SUFFIX=CTCF/Snyder_CTCF_GM12878_PE
mkdir -p ${WORK_DIR}/$SUFFIX; cd ${WORK_DIR}/$SUFFIX;
bds /users/leepc12/code/bds_atac/chipseq/chipseq.bds \
-fastq1_1 ${DATA_ROOT}/CTCF/PE/Snyder_CTCF_GM12878_PE_Rep1_1.fastq.gz \
-fastq1_2 ${DATA_ROOT}/CTCF/PE/Snyder_CTCF_GM12878_PE_Rep1_2.fastq.gz \
-fastq2_1 ${DATA_ROOT}/CTCF/PE/Snyder_CTCF_GM12878_PE_Rep2_1.fastq.gz \
-fastq2_2 ${DATA_ROOT}/CTCF/PE/Snyder_CTCF_GM12878_PE_Rep2_2.fastq.gz \
-ctl_fastq1_1 ${DATA_ROOT}/CTCF/PE/Snyder_Input_GM12878_PE_1.fastq.gz \
-ctl_fastq1_2 ${DATA_ROOT}/CTCF/PE/Snyder_Input_GM12878_PE_2.fastq.gz \
-url_base http://mitra.stanford.edu/kundaje/leepc12/chipseq_test/${SUFFIX}/out \
-species hg19 -bwa_idx ${DATA_ROOT}/CTCF/bwa_index/male.hg19.fa.gz -chrsz ${DATA_ROOT}/CTCF/male.hg19.chrom.sizes
screen -RD UCD_GATA2_K562_SE
DATA_ROOT=/srv/scratch/leepc12/data; WORK_DIR=/srv/scratch/leepc12/run/chipseq_test; SUFFIX=GATA2/UCD_GATA2_K562_SE
mkdir -p ${WORK_DIR}/$SUFFIX; cd ${WORK_DIR}/$SUFFIX;
bds /users/leepc12/code/bds_atac/chipseq/chipseq.bds \
-species hg19 -bwa_idx ${DATA_ROOT}/CTCF/bwa_index/male.hg19.fa.gz -chrsz ${DATA_ROOT}/CTCF/male.hg19.chrom.sizes \
-url_base http://mitra.stanford.edu/kundaje/leepc12/chipseq_test/${SUFFIX}/out \
-fastq1 ${DATA_ROOT}/GATA2/SE/HAIB_GATA2_K562_SE/HAIB_GATA2_K562_SE_Rep1.fastq.gz \
-fastq2 ${DATA_ROOT}/GATA2/SE/HAIB_GATA2_K562_SE/HAIB_GATA2_K562_SE_Rep2.fastq.gz \
-ctl_fastq1 ${DATA_ROOT}/GATA2/SE/HAIB_GATA2_K562_SE/HAIB_Input_K562_SE_Rep1.fastq.gz \
-ctl_fastq2 ${DATA_ROOT}/GATA2/SE/HAIB_GATA2_K562_SE/HAIB_Input_K562_SE_Rep2.fastq.gz
screen -RD UChicago_GATA2_K562_SE
DATA_ROOT=/srv/scratch/leepc12/data; WORK_DIR=/srv/scratch/leepc12/run/chipseq_test; SUFFIX=GATA2/UChicago_GATA2_K562_SE
mkdir -p ${WORK_DIR}/$SUFFIX; cd ${WORK_DIR}/$SUFFIX;
bds /users/leepc12/code/bds_atac/chipseq/chipseq.bds \
-species hg19 -bwa_idx ${DATA_ROOT}/CTCF/bwa_index/male.hg19.fa.gz -chrsz ${DATA_ROOT}/CTCF/male.hg19.chrom.sizes \
-url_base http://mitra.stanford.edu/kundaje/leepc12/chipseq_test/${SUFFIX}/out \
-fastq1 ${DATA_ROOT}/GATA2/SE/UChicago_GATA2_K562_SE/UChicago_GATA2_K562_SE_Rep1.fastq.gz \
-fastq2 ${DATA_ROOT}/GATA2/SE/UChicago_GATA2_K562_SE/UChicago_GATA2_K562_SE_Rep2.fastq.gz \
-ctl_fastq1 ${DATA_ROOT}/GATA2/SE/UChicago_GATA2_K562_SE/UChicago_Input_K562_SE_Rep1.fastq.gz \
-ctl_fastq2 ${DATA_ROOT}/GATA2/SE/UChicago_GATA2_K562_SE/UChicago_Input_K562_SE_Rep2.fastq.gz
screen -RD Snyder_CTCF_GM12878_SE
DATA_ROOT=/srv/scratch/leepc12/data; WORK_DIR=/srv/scratch/leepc12/run/chipseq_test; SUFFIX=CTCF/Snyder_CTCF_GM12878_SE
mkdir -p ${WORK_DIR}/$SUFFIX; cd ${WORK_DIR}/$SUFFIX;
bds /users/leepc12/code/bds_atac/chipseq/chipseq.bds \
-species hg19 -bwa_idx ${DATA_ROOT}/CTCF/bwa_index/male.hg19.fa.gz -chrsz ${DATA_ROOT}/CTCF/male.hg19.chrom.sizes \
-url_base http://mitra.stanford.edu/kundaje/leepc12/chipseq_test/${SUFFIX}/out \
-fastq1 ${DATA_ROOT}/CTCF/SE/Snyder_CTCF_GM12878_SE_Rep1.fastq.gz \
-fastq2 ${DATA_ROOT}/CTCF/SE/Snyder_CTCF_GM12878_SE_Rep2.fastq.gz \
-ctl_fastq1 ${DATA_ROOT}/CTCF/SE/Snyder_Input_GM12878_SE.fastq.gz
screen -RD HudsonAlpha
DATA_ROOT=/srv/scratch/leepc12/data; WORK_DIR=/srv/scratch/leepc12/run/chipseq_test; SUFFIX=YY1/HudsonAlpha
mkdir -p ${WORK_DIR}/$SUFFIX; cd ${WORK_DIR}/$SUFFIX;
bds /users/leepc12/code/bds_atac/chipseq/chipseq.bds \
-species hg19 -bwa_idx ${DATA_ROOT}/CTCF/bwa_index/male.hg19.fa.gz -chrsz ${DATA_ROOT}/CTCF/male.hg19.chrom.sizes \
-url_base http://mitra.stanford.edu/kundaje/leepc12/chipseq_test/${SUFFIX}/out \
-fastq1 ${DATA_ROOT}/YY1/HudsonAlpha/ENCFF000OHH.fastq.gz \
-fastq2 ${DATA_ROOT}/YY1/HudsonAlpha/ENCFF000OHO.fastq.gz \
-ctl_fastq1 ${DATA_ROOT}/YY1/HudsonAlpha/ENCFF000ODP.fastq.gz
screen -RD Sydh
DATA_ROOT=/srv/scratch/leepc12/data; WORK_DIR=/srv/scratch/leepc12/run/chipseq_test; SUFFIX=YY1/Sydh
mkdir -p ${WORK_DIR}/$SUFFIX; cd ${WORK_DIR}/$SUFFIX;
bds /users/leepc12/code/bds_atac/chipseq/chipseq.bds \
-url_base http://mitra.stanford.edu/kundaje/leepc12/chipseq_test/${SUFFIX}/out \
-species hg19 -bwa_idx ${DATA_ROOT}/CTCF/bwa_index/male.hg19.fa.gz -chrsz ${DATA_ROOT}/CTCF/male.hg19.chrom.sizes \
-fastq1 ${DATA_ROOT}/YY1/Sydh/ENCFF000WGS.fastq.gz \
-fastq2 ${DATA_ROOT}/YY1/Sydh/ENCFF000WGT.fastq.gz \
-ctl_fastq1 ${DATA_ROOT}/YY1/Sydh/ENCFF000VWV.fastq.gz
WORK_DIR_OLD=/srv/scratch/leepc12/run/TF_chipseq_pipeline_test/$SUFFIX
mkdir -p $WORK_DIR/$SUFFIX/out/align/rep1; cp -pr $WORK_DIR_OLD/out/align_rep1/*.bam $WORK_DIR/$SUFFIX/out/align/rep1/
mkdir -p $WORK_DIR/$SUFFIX/out/align/rep1; cp -pr $WORK_DIR_OLD/out/align_rep1/*.flagstat.qc $WORK_DIR/$SUFFIX/out/align/rep1/
mkdir -p $WORK_DIR/$SUFFIX/out/align/rep2; cp -pr $WORK_DIR_OLD/out/align_rep2/*.bam $WORK_DIR/$SUFFIX/out/align/rep2/
mkdir -p $WORK_DIR/$SUFFIX/out/align/rep2; cp -pr $WORK_DIR_OLD/out/align_rep2/*.flagstat.qc $WORK_DIR/$SUFFIX/out/align/rep2/
mkdir -p $WORK_DIR/$SUFFIX/out/align/ctl1; cp -pr $WORK_DIR_OLD/out/align_ctl_rep1/*.bam $WORK_DIR/$SUFFIX/out/align/ctl1/
mkdir -p $WORK_DIR/$SUFFIX/out/align/ctl1; cp -pr $WORK_DIR_OLD/out/align_ctl_rep1/*.flagstat.qc $WORK_DIR/$SUFFIX/out/align/ctl1/
mkdir -p $WORK_DIR/$SUFFIX/out/align/ctl2; cp -pr $WORK_DIR_OLD/out/align_ctl_rep2/*.bam $WORK_DIR/$SUFFIX/out/align/ctl2/
mkdir -p $WORK_DIR/$SUFFIX/out/align/ctl2; cp -pr $WORK_DIR_OLD/out/align_ctl_rep2/*.flagstat.qc $WORK_DIR/$SUFFIX/out/align/ctl2/
find . -name '*.nmsrt.bam' -delete
find . -name '*.nodup.bam' -delete