annotate lib/Report.pl @ 3:17ce4f3bffa2 default tip

Uploaded
author jesse-erdmann
date Tue, 24 Jan 2012 18:33:41 -0500
parents 1437a2df99c0
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
1 require 'config.pl';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
2
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
3
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
4
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
5 $sth = $dbh->prepare("drop table if exists lib_m_$proj");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
6 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
7 $sth = $dbh->prepare("create table lib_m_$proj select distinct library, sum(number) as total from illumina4dec2_$proj group by library");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
8 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
9 $sth = $dbh->prepare("update lib_m_$proj set library = left(library,LENGTH(library)-2) where library like '%-L'");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
10 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
11 $sth = $dbh->prepare("update lib_m_$proj set library = left(library,LENGTH(library)-2) where library like '%-R'");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
12 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
13 $sth = $dbh->prepare("drop table if exists lib_mappable_$proj");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
14 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
15 $sth = $dbh->prepare("create table lib_mappable_$proj select library, sum(total) as total from lib_m_$proj group by library");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
16 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
17
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
18 print "completed mapping!!!";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
19
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
20 ###
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
21 #Generate a report describing the mapping.
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
22 ###
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
23
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
24 open OUT, "> results/project_summary.txt";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
25 print OUT "Summary of project $proj\n\nA.Report of project based counts \n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
26
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
27 $sth = $dbh->prepare("select count(*) from illumina_raw_$proj");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
28 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
29
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
30 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
31 print OUT $row[0]." Total number of Seqeunces\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
32 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
33
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
34 $sth = $dbh->prepare("select count(*) from illumina_decoded_$proj;");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
35 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
36
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
37 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
38 print OUT $row[0]." Total number of Sequences that map to a barcode\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
39 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
40
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
41 $sth = $dbh->prepare("select count(*) from illumina_without_IRDR_$proj");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
42 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
43
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
44 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
45 print OUT $row[0]." Total number of Sequences with a barcode that also have an IRDR\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
46 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
47
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
48 $sth = $dbh->prepare("select count(*) from illumina_without_IRDR_$proj where type = 'idR';");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
49 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
50
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
51 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
52 print OUT $row[0]." Total number of Ideal Right linker sequences\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
53 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
54
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
55 $sth = $dbh->prepare("select count(*) from illumina_without_IRDR_$proj where type = 'idL';");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
56 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
57
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
58 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
59 print OUT $row[0]." Total number of Ideal left linkersequences\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
60 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
61
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
62 $sth = $dbh->prepare("select count(*) from illumina_without_IRDR_$proj where type = 'good';");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
63 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
64
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
65 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
66 print OUT $row[0]." Total number of good sequences\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
67 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
68
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
69 $sth = $dbh->prepare("select count(*) from illumina_without_IRDR_$proj where type = 'bam';");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
70 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
71
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
72 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
73 print OUT $row[0]." Total number of Bamh1 sequence Artifact (removed)\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
74 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
75
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
76 $sth = $dbh->prepare("select count(*) from illumina_without_IRDR_$proj where type = 'noTA';");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
77 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
78
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
79 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
80 print OUT $row[0]." Total number of sequences removed because of lack of TA sequence)\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
81 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
82
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
83 $sth = $dbh->prepare("select count(*) from illumina_without_IRDR_$proj where type = 'min';");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
84 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
85
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
86 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
87 print OUT $row[0]." Total number of sequences removed because of sequence length < 24)\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
88 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
89
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
90 $sth = $dbh->prepare("select count(*) from illumina4dec2_$proj;");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
91 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
92
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
93 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
94 print OUT $row[0]." Total number of unique seqeunces pryor to mapping\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
95 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
96
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
97 $sth = $dbh->prepare("select sum(number) from bowtie_$proj;");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
98 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
99 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
100 print OUT $row[0]." Total SUM of seqeunces that mapped\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
101 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
102
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
103 $sth = $dbh->prepare("select count(*) from bowtie_$proj;");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
104 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
105
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
106 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
107 print OUT $row[0]." Number of unique sequences that mapped\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
108 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
109
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
110 $sth = $dbh->prepare("select count(*) from bowtie_lib_$proj;");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
111 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
112
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
113 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
114 print OUT $row[0]." Number of distinct locations mapped\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
115 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
116
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
117 $sth = $dbh->prepare("select count(*) from illumina_hist_$proj;");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
118 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
119
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
120 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
121 print OUT $row[0]." Number of distinct regions mapped\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
122 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
123 close OUT;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
124
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
125 #report part B
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
126 open OUT, ">> results/project_summary.txt";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
127 print OUT"\n\nB. Library counts associated with project $proj\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
128
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
129 %sequence='';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
130 %barcode_count='';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
131 %IRDR_count='';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
132 %IRDR_good='';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
133 %IRDR_unique='';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
134 %map_count='';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
135 %map_total='';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
136 %nr_count='';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
137 @library='';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
138
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
139 $sth = $dbh->prepare("select * from barcode_$proj");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
140 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
141
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
142 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
143 $sequence{$row[1]} = $row[0];
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
144 push(@library,$row[1]);
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
145 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
146
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
147 $sth = $dbh->prepare("select library, count(id) from illumina_decoded_$proj group by library");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
148 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
149
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
150 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
151 $barcode_count{$row[0]} = $row[1];
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
152 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
153
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
154 $sth = $dbh->prepare("select library, count(id) from illumina_without_IRDR_$proj group by library");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
155 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
156
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
157 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
158 $IRDR_count{$row[0]} = $row[1];
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
159 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
160
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
161 $sth = $dbh->prepare("select library, sum(number) from illumina4dec2_$proj group by library;");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
162 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
163
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
164 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
165 $IRDR_good{$row[0]} = $row[1];
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
166 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
167
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
168 $sth = $dbh->prepare("select library, count(number) from illumina4dec2_$proj group by library;");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
169 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
170
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
171 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
172 $IRDR_unique{$row[0]} = $row[1];
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
173 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
174
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
175 $sth = $dbh->prepare("select library, sum(number) from bowtie_$proj group by library");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
176 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
177
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
178 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
179 $map_count{$row[0]} = $row[1];
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
180 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
181
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
182 $sth = $dbh->prepare("select library, count(id) from bowtie_$proj group by library");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
183 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
184
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
185 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
186 $map_total{$row[0]} = $row[1];
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
187 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
188
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
189 $sth = $dbh->prepare("select library, count(chromo) from bowtie_lib_$proj group by library");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
190 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
191
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
192 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
193 $map_nr{$row[0]} = $row[1];
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
194 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
195
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
196 print OUT "Library\tBarcode Sequence\tBarcode\tIRDR\tMappable\tUnique mappable\tTotal map\tUnique map\tNR map";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
197 foreach $item (@library) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
198 print OUT "$item\t$sequence{$item}\t$barcode_count{$item}\t$IRDR_count{$item}\t$IRDR_good{$item}\t$IRDR_unique{$item}\t$map_count{$item}\t$map_total{$item}\t$map_nr{$item}\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
199 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
200 close OUT;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
201
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
202 open OUT, ">> results/project_summary.txt";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
203 print OUT "\n\nC. Regions associated with project $proj\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
204
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
205 %mappable='';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
206 %name='';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
207 %count01='';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
208 %count001='';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
209 %count0001='';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
210 %count0='';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
211 @library='';
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
212
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
213 $sth = $dbh->prepare("select * from lib_count_$proj");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
214 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
215
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
216 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
217 $name{$row[0]} = $row[1];
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
218 push(@library,$row[0]);
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
219 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
220
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
221 $sth = $dbh->prepare("select * from lib_mappable_$proj");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
222 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
223
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
224 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
225 $mappable{$row[0]} = $row[1];
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
226 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
227
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
228 $sth = $dbh->prepare("select A.library, count(chromo) from illumina_hist_$proj A, lib_mappable_$proj B where A.library = B.library and count/total > .01 group by A.library order by chromo,start+0");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
229 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
230
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
231 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
232 $count01{$row[0]} = $row[1];
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
233 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
234
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
235 $sth = $dbh->prepare("select A.library, count(chromo) from illumina_hist_$proj A, lib_mappable_$proj B where A.library = B.library and count/total > .001 group by A.library order by chromo,start+0");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
236 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
237
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
238 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
239 $count001{$row[0]} = $row[1];
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
240 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
241
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
242 $sth = $dbh->prepare("select A.library, count(chromo) from illumina_hist_$proj A, lib_mappable_$proj B where A.library = B.library and count/total > .0001 group by A.library order by chromo,start+0");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
243 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
244
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
245 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
246 $count0001{$row[0]} = $row[1];
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
247 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
248
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
249 $sth = $dbh->prepare("select A.library, count(chromo) from illumina_hist_$proj A, lib_mappable_$proj B where A.library = B.library and count/total > 0 group by A.library order by chromo,start+0");
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
250 $sth->execute;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
251
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
252 print OUT "Library\tTotal mappable\ttotal map\tLocations >0.01\tLocations >0.001\tLocations >0.0001\t All mapped Locations";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
253 while ((@row) = $sth->fetchrow_array) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
254 $count0{$row[0]} = $row[1];
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
255 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
256
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
257 foreach $item (@library) {
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
258 print OUT "$item\t$mappable{$item}\t$name{$item}\t$count01{$item}\t$count001{$item}\t$count0001{$item}\t$count0{$item}\n";
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
259 }
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
260 close OUT;
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
261
1437a2df99c0 Uploaded
jesse-erdmann
parents:
diff changeset
262