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
|
#!/bin/bash
cd "$(dirname "${BASH_SOURCE[0]}")"
echo -n '.'
weblogo --format png_print --size large \
-i -5 -l 1 -u 20 \
--title "The DNA-binding helix-turn-helix motif of the CAP family" \
< cap_hth.fa \
> cap_hth.png || exit
echo -n '.'
weblogo --format png_print --size large \
--title "58 CAP Binding Sites" \
-i -10 \
< cap_dna.fa \
> cap_dna.png || exit
echo -n '.'
weblogo --format png_print --size large \
--title "19 LexA Binding Sites" \
-i -9 \
< lexA.fa \
> lexA.png || exit
echo -n '.'
weblogo --format png_print --size large \
--title "-10 region of E. coli promoters" \
-i -21 --lower 0 -u 7 \
< ecoli10.fa \
> ecoli10.png || exit
echo -n '.'
weblogo --format png_print --size large \
-l 63 -u 83 \
< globins.fa \
> globins.png || exit
#echo -n '.'
#weblogo --format PNG --size large \
# -l 31 -u 150 \
# < globins.fa \
# > more_globins.png || exit
echo -n '.'
weblogo --format png_print --size large \
--title "Helix-Turn-Helix Motifs" \
-i -11 -l 1 -u 17 --yaxis 3.5 \
< hth.fa \
> hth.png || exit
echo -n '.'
weblogo --format png_print --size large \
--title "exon|intron" \
-i -11 -l -6 -u 8 \
< exon-intron.fa \
> exon-intron.png || exit
echo -n '.'
weblogo --format png_print --size large \
--title "intron | exon" \
-i -21 -l -20 -u 3 \
< intron-exon.fa \
> intron-exon.png || exit
echo
|
