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
|
BEGIN {
while (getline < TrLengthFile)
{
T[$1]=$2+0;
}
}
{
if (NR%2==1 && $3 in T)
{
s=$4;
if (s==0)
next;
i=$9;
l=T[$3];
n=substr($12,6);
#print s,i,l,n > "alignTr.txt";
iN[i]+=1/n;
t1=int((s-1)/l*100)+1;t2=int((s+i-2)/l*100)+1;
for (ii=t1;ii<=t2;ii++)
C[ii]+=1/n;
C5[t1]+=1/n;
C3[t2]+=1/n;
if (n==1)
{
for (ii=t1;ii<=t2;ii++)
Cu[ii]+=1;
C5u[t1]++;
C3u[t2]++;
iNu[i]++;
}
}
}
END {
for (ii in iN)
print ii,iN[ii] > "insertHist.txt";
for (ii in iNu)
print ii,iNu[ii] > "insertHistUnique.txt";
for (ii=1;ii<=101;ii++)
print ii,C[ii]+0,Cu[ii]+0,C5[ii]+0,C3[ii]+0,C5u[ii]+0,C3u[ii]+0 > "coverage.txt"
}
|
