HCAP
HCAP (Homogneous Composition Among Positions) is a command line program written in Java to crop high-throughput short (Illumina) sequencing reads in order to observe homogeneous base composition among all read positions. Although this preprocessing step is generally not required (see e.g. here), HCAP was developed for dealing with unexpectedly biased composition that can be sometimes observed at the 5' and/or 3' ends (see e.g. here).
Installation
Clone this repository with the following command line:
git clone https://gitlab.pasteur.fr/GIPhy/HCAP.gitCompilation and execution
The source code of HCAP is inside the src directory. It requires Java 11 (or higher) and can be compiled and executed in two different ways.
Building an executable jar file
On computers with Oracle JDK (11 or higher) installed, a Java executable jar file can be created. In a command-line window, go to the src directory and type:
javac HCAP.java
echo Main-Class: HCAP > MANIFEST.MF
jar -cmvf MANIFEST.MF HCAP.jar HCAP.class
rm MANIFEST.MF HCAP.classThis will create the executable jar file HCAP.jar that can be run with the following command line model:
java -jar HCAP.jar [options]Building a native code binary
On computers with GraalVM installed, a native executable can be built. In a command-line window, go to the src directory, and type:
javac HCAP.java
native-image HCAP HCAP
rm HCAP.classThis will create the native executable HCAP that can be run with the following command line model:
./HCAP [options]Usage
Run HCAP without option to read the following documentation:
HCAP
Cropping FASTQ short reads to reach homogeneous composition among positions
USAGE:
HCAP -i <FASTQ> -o <basefile> [-z] [-l <int>] [-n] [-v] [-h]
HCAP -1 <FASTQ.R1> -2 <FASTQ.R2> -o <basefile> [-z] [-l <int>] [-n] [-v] [-h]
OPTIONS:
-i <infile> FASTQ-formatted input file; filename should end with .gz when gzipped
-i <infile> FASTQ-formatted R1 input file; filename should end with .gz when gzipped
-i <infile> FASTQ-formatted R2 input file; filename should end with .gz when gzipped
-o <string> outfile basename (mandatory)
-z gzipped output file(s) (default: not set)
-c z-score cutoff (defaut: 1.64)
-n replacing read prefix/suffix to crop with Ns (default: not set)
-l <integer> minimum length of output reads (default: 50)
-v verbose mode (default: not set)
-h prints this help and exitMethod
-
Given a high-throughput sequencing read of length ℓ, let np(r) be the number of the observed bases r = A, C, G, T at each sequenced position p = 1, ..., ℓ, and np = np(A) + np(C) + np(G) + np(T). As each of the four proportions f(A), f(C), f(G), f(T) should be almost identical among every sequenced position p, this expectation can be assessed using a one proportion z-test, i.e. |zp(r)| ≤ c, where zp(r) := [np(r) / np − f(r)] [f(r) (1 − f(r)) / np]−1/2 and c is a fixed cutoff (e.g. c = 1.64). In consequence, every large value zp(r) indicates that the base r is overrepresented at position p.
-
First, HCAP parses the s (= 8,192) first reads to determine (i) the (maximum) read length ℓ and (ii) an estimate of the proportions fp(r) = np(r) / np for each sequenced position p = 1, ..., ℓ. For each r = A, C, G, T, the expected proportion f(r) is then approximated by the median of the ℓ observed proportions fp(r).
-
To obtain an overall homogeneous composition among positions, HCAP crops the read prefixes and/or suffixes that are mainly constituted of bases inducing large z-scores values zp, i.e. every value zp(r) associated to the base r at position p. Therefore, the first b bases of a read should be cropped when zp > c and F5(b) := Σp=1...b zp > 0, whereas the last b bases should be cropped when zp > c and F3(b) := Σp=b...ℓ zp > 0. Of note, every negative value zp is multiplied by a small constant (= 1.5) to avoid a too aggressive cropping. In practice, HCAP computes b5 := argmaxb, zb > 0 F5(b) and b3 := argmaxb, zb > 0 F3(b), and crops the corresponding prefix and/or suffix when F5(b5) > 0 and/or F3(b3) > 0, respectively.
-
Second, for each the s first reads, when a cropping is performed (according to the above criteria), every corresponding value np(r) is decremented (i.e. for p = 1, ..., b5 and/or p = b3 ..., ℓ). To avoid a too aggressive cropping for the first s reads, the cutoff is defined as a decreasing function, varying from 10 to c depending on the nth parsed read. After processing the s first reads, the four expected proportions f(r) are updated, i.e. median of the ℓ observed proportions fp(r).
-
Finally, for each of the remaining reads, every corresponding value np(r) is first incremented (i.e. for p = 1, ..., ℓ); next, if the read is cropped (according to the above criteria using the specified cutoff), then every corresponding value np(r) is decremented (i.e. for p = 1, ..., b5 and/or p = b3 ..., ℓ). The four expected proportions f(r) are updated every s parsed reads.
-
Following this read cropping procedure, the ℓ + 1 proportions fp(r) and f(r) generally converge to the same value for each r = A, C, G, T, as expected.
Notes
-
Input files are FASTQ-formatted ones (input files compressed using gzip should ends with the prefix .gz). Single-ends FASTQ files should be specified using option
-i, whereas paired-ends FASTQ files should be specified using options-1and-2. The prefix of the output files should be specified using option-o. Reading gzipped FASTQ files does not seem to have a negative impact on the overall running times. However, writing gzipped FASTQ output files (option-z) requires approximately 5 times slower running times. -
It is worth noting that when the input reads are of different lengths, HCAP considers they were clipped only at the 3' ends (e.g. technical oligonucleotide sequence removed during demultiplexing). If the input reads were trimmed/cliiped at the 5' ends, HCAP results will be strongly biased (as the initial sequenced positions can not be determined).
-
The z-score cutoff is set to c = 1.64 by default. Increasing this value (e.g. 3) enables a less aggressive cropping.
-
Every cropped position can be replaced with 'N' using the option
-n. -
Every read of length lower than a specified threshold (50, by default; option
-l) is discarded. This criteria is applied before and after cropping. When dealing with paired-end data, if one of the read is too short, then the whole pair is discarded. -
The verbose mode (option -v) can be used to observe the trimming details.
Example
The following Bash command lines enable to download the two compressed FASTQ files associated to the run accession ERR6798852:
EBIURL="https://ftp.sra.ebi.ac.uk/vol1/fastq";
wget $EBIURL/ERR679/002/ERR6798852/ERR6798852_1.fastq.gz ;
wget $EBIURL/ERR679/002/ERR6798852/ERR6798852_2.fastq.gz ;The downloaded FASTQ files ERR6798852_1.fastq.gz (47 Mb) and ERR6798852_2.fastq.gz (52 Mb) correspond to the 2×250 bps whole genome sequencing of a Listeria monocytogenes strain.
Associated descriptive statistics were computed using fastq_info. The two obtained figures (see below) show a typical biased composition at the 5' ends (positions 1-16), as well as the absence of the base C at position 251.
 |
 |
To obtain homogeneous composition among all sequenced positions, read cropping can be performed using HCAP with the following command line:
HCAP -1 ERR6798852_1.fastq.gz -2 ERR6798852_2.fastq.gz -o ERR6798852The above command line displays some details (see below) that confirm important deviations from the expected proportions f(A), f(C), f(G), f(T) at the ends (e.g. positions 1-14 and 251 induce z-scores ≫ 1.64).
=== initial estimates =====================================================================================================================================
[ fA fC fG fT R1] [ fA fC fG fT R2]
median 0.303 0.197 0.197 0.304 0.304 0.196 0.198 0.303
pos fA fC fG fT zA zC zG zT fA fC fG fT zA zC zG zT
1 0.264 0.219 0.359 0.158 -7.72 4.99 37.11 -28.66 0.261 0.219 0.364 0.156 -8.48 5.38 37.85 -28.96
2 0.232 0.196 0.137 0.435 -13.93 -0.18 -13.53 25.77 0.225 0.196 0.135 0.445 -15.57 -0.02 -14.34 28.03
3 0.294 0.209 0.159 0.337 -1.66 2.87 -8.52 6.54 0.300 0.206 0.155 0.339 -0.70 2.46 -9.73 7.01
4 0.188 0.240 0.116 0.456 -22.73 9.82 -18.31 30.04 0.189 0.249 0.119 0.444 -22.59 12.12 -17.95 27.69
5 0.400 0.107 0.113 0.380 19.09 -20.44 -18.95 14.98 0.391 0.108 0.117 0.384 17.15 -20.10 -18.22 15.98
6 0.444 0.118 0.234 0.205 27.79 -18.02 8.44 -19.48 0.439 0.120 0.236 0.206 26.64 -17.37 8.64 -19.15
7 0.355 0.159 0.213 0.273 10.24 -8.52 3.71 -6.07 0.352 0.160 0.216 0.272 9.49 -8.04 4.20 -6.19
8 0.409 0.155 0.202 0.233 20.96 -9.49 1.29 -13.86 0.405 0.155 0.195 0.245 19.89 -9.27 -0.49 -11.48
9 0.191 0.318 0.227 0.264 -21.96 27.47 6.88 -7.75 0.191 0.324 0.217 0.268 -22.25 29.25 4.50 -6.89
10 0.283 0.218 0.157 0.342 -3.95 4.79 -9.02 7.60 0.278 0.224 0.155 0.343 -5.02 6.36 -9.68 7.93
11 0.290 0.153 0.293 0.265 -2.63 -9.96 21.87 -7.66 0.286 0.152 0.291 0.271 -3.56 -9.96 21.32 -6.31
12 0.277 0.199 0.187 0.337 -5.20 0.54 -2.10 6.54 0.278 0.205 0.190 0.327 -5.05 2.07 -1.74 4.78
13 0.239 0.216 0.254 0.291 -12.63 4.43 13.14 -2.56 0.234 0.211 0.244 0.311 -13.65 3.43 10.44 1.65
14 0.228 0.280 0.183 0.309 -14.75 19.02 -3.13 0.99 0.228 0.283 0.187 0.303 -15.00 19.87 -2.44 -0.03
15 0.318 0.204 0.223 0.255 3.05 1.60 5.99 -9.60 0.324 0.205 0.220 0.250 4.06 2.21 5.11 -10.40
16 0.329 0.190 0.197 0.284 5.16 -1.63 0.12 -3.86 0.339 0.197 0.190 0.275 6.92 0.23 -1.80 -5.56
17 0.325 0.190 0.177 0.308 4.42 -1.49 -4.49 0.75 0.328 0.182 0.178 0.312 4.85 -3.17 -4.46 1.75
18 0.324 0.182 0.187 0.307 4.23 -3.49 -2.10 0.61 0.328 0.177 0.186 0.309 4.68 -4.31 -2.55 1.24
19 0.312 0.185 0.193 0.310 1.80 -2.71 -0.74 1.19 0.309 0.195 0.188 0.307 1.03 -0.08 -2.13 0.88
20 0.314 0.191 0.190 0.304 2.28 -1.32 -1.43 0.10 0.321 0.185 0.189 0.305 3.31 -2.41 -1.91 0.42
21 0.314 0.190 0.192 0.305 2.09 -1.57 -1.04 0.18 0.305 0.192 0.195 0.308 0.31 -0.83 -0.60 0.93
22 0.299 0.196 0.203 0.302 -0.73 -0.18 1.46 -0.38 0.306 0.202 0.191 0.301 0.36 1.46 -1.49 -0.32
23 0.305 0.198 0.183 0.314 0.45 0.29 -3.05 1.93 0.300 0.194 0.199 0.307 -0.72 -0.33 0.20 0.83
24 0.309 0.183 0.197 0.311 1.17 -3.10 0.10 1.43 0.302 0.195 0.201 0.302 -0.34 -0.16 0.76 -0.18
25 0.297 0.198 0.194 0.310 -1.11 0.29 -0.52 1.31 0.302 0.193 0.202 0.303 -0.39 -0.58 1.03 -0.01
26 0.304 0.196 0.189 0.311 0.26 -0.27 -1.63 1.38 0.308 0.195 0.189 0.307 0.86 -0.10 -1.88 0.86
27 0.314 0.186 0.193 0.307 2.25 -2.57 -0.71 0.58 0.307 0.188 0.200 0.305 0.60 -1.75 0.64 0.35
28 0.307 0.194 0.192 0.307 0.76 -0.54 -1.13 0.68 0.305 0.199 0.194 0.302 0.31 0.65 -0.71 -0.25
29 0.304 0.195 0.194 0.308 0.11 -0.38 -0.68 0.80 0.310 0.195 0.190 0.306 1.22 -0.21 -1.82 0.54
30 0.315 0.190 0.199 0.295 2.47 -1.52 0.60 -1.67 0.311 0.193 0.195 0.301 1.49 -0.60 -0.63 -0.42
31 0.307 0.203 0.197 0.293 0.74 1.40 0.15 -2.08 0.310 0.198 0.194 0.299 1.29 0.43 -0.94 -0.85
32 0.303 0.204 0.194 0.299 -0.08 1.74 -0.63 -0.88 0.299 0.201 0.192 0.308 -1.01 1.23 -1.24 1.02
33 0.308 0.188 0.205 0.299 0.98 -2.10 2.01 -0.90 0.310 0.193 0.194 0.304 1.25 -0.72 -0.85 0.11
34 0.306 0.193 0.203 0.299 0.55 -0.99 1.51 -1.00 0.302 0.207 0.189 0.301 -0.27 2.63 -1.96 -0.30
35 0.312 0.198 0.181 0.309 1.85 0.18 -3.63 1.14 0.294 0.200 0.199 0.307 -1.93 1.01 0.31 0.78
36 0.305 0.191 0.197 0.308 0.33 -1.43 0.12 0.80 0.300 0.194 0.201 0.305 -0.70 -0.49 0.76 0.47
37 0.305 0.192 0.199 0.304 0.33 -1.02 0.51 0.10 0.305 0.198 0.194 0.303 0.33 0.59 -0.88 -0.08
38 0.301 0.201 0.194 0.304 -0.37 0.90 -0.46 -0.02 0.311 0.191 0.195 0.303 1.37 -0.99 -0.63 0.04
39 0.309 0.189 0.197 0.305 1.20 -1.79 0.15 0.22 0.305 0.192 0.194 0.309 0.21 -0.88 -0.85 1.29
40 0.312 0.201 0.189 0.298 1.75 1.01 -1.74 -1.12 0.314 0.188 0.201 0.298 2.01 -1.86 0.67 -1.00
41 0.303 0.196 0.200 0.301 0.07 -0.15 0.71 -0.54 0.313 0.190 0.192 0.305 1.75 -1.33 -1.24 0.47
42 0.302 0.190 0.203 0.306 -0.27 -1.66 1.54 0.37 0.303 0.189 0.203 0.305 -0.08 -1.58 1.28 0.33
43 0.295 0.198 0.204 0.303 -1.54 0.26 1.71 -0.16 0.311 0.192 0.199 0.297 1.46 -0.77 0.39 -1.14
44 0.314 0.194 0.189 0.303 2.28 -0.77 -1.68 -0.16 0.305 0.201 0.198 0.296 0.26 1.21 0.01 -1.31
45 0.302 0.187 0.199 0.312 -0.15 -2.21 0.51 1.62 0.310 0.191 0.196 0.303 1.32 -1.08 -0.38 -0.06
46 0.300 0.197 0.198 0.304 -0.51 0.07 0.37 0.13 0.304 0.198 0.197 0.301 -0.00 0.62 -0.22 -0.35
47 0.312 0.195 0.187 0.306 1.82 -0.49 -2.10 0.42 0.315 0.193 0.190 0.302 2.18 -0.55 -1.77 -0.18
48 0.313 0.188 0.189 0.309 2.06 -1.93 -1.71 1.09 0.314 0.196 0.183 0.307 1.97 0.15 -3.24 0.71
49 0.304 0.190 0.200 0.306 0.14 -1.54 0.85 0.46 0.302 0.196 0.202 0.301 -0.34 -0.02 0.89 -0.42
50 0.305 0.200 0.188 0.306 0.50 0.76 -1.88 0.46 0.311 0.192 0.192 0.306 1.34 -0.91 -1.27 0.54
51 0.305 0.194 0.204 0.296 0.42 -0.55 1.76 -1.46 0.309 0.188 0.190 0.313 1.10 -1.81 -1.75 1.98
52 0.301 0.196 0.201 0.302 -0.39 -0.18 1.10 -0.40 0.305 0.198 0.195 0.302 0.24 0.54 -0.49 -0.27
53 0.301 0.198 0.197 0.304 -0.33 0.22 0.05 0.10 0.311 0.199 0.181 0.310 1.36 0.69 -3.88 1.40
54 0.303 0.191 0.202 0.304 -0.05 -1.23 1.28 0.01 0.307 0.189 0.203 0.301 0.55 -1.40 1.13 -0.32
55 0.307 0.200 0.193 0.300 0.81 0.74 -0.85 -0.71 0.295 0.198 0.203 0.305 -1.74 0.43 1.21 0.33
56 0.301 0.202 0.198 0.298 -0.30 1.14 0.41 -1.04 0.301 0.204 0.187 0.309 -0.64 1.83 -2.35 1.09
57 0.312 0.189 0.195 0.304 1.77 -1.71 -0.43 0.09 0.305 0.190 0.200 0.305 0.32 -1.39 0.57 0.38
58 0.295 0.206 0.199 0.300 -1.48 2.00 0.60 -0.77 0.309 0.197 0.194 0.300 0.98 0.40 -0.85 -0.59
59 0.299 0.192 0.198 0.312 -0.81 -1.15 0.24 1.60 0.309 0.201 0.184 0.306 1.08 1.22 -3.17 0.61
60 0.302 0.193 0.201 0.304 -0.15 -0.94 1.10 0.01 0.301 0.196 0.203 0.300 -0.58 0.01 1.27 -0.52
61 0.301 0.189 0.201 0.309 -0.33 -1.70 0.96 0.97 0.309 0.192 0.201 0.298 0.96 -0.76 0.67 -0.89
62 0.298 0.198 0.197 0.307 -0.97 0.16 0.19 0.67 0.306 0.193 0.194 0.306 0.49 -0.52 -0.76 0.62
63 0.304 0.189 0.195 0.312 0.15 -1.73 -0.41 1.69 0.308 0.188 0.204 0.300 0.74 -1.71 1.49 -0.56
64 0.301 0.207 0.193 0.299 -0.30 2.26 -0.86 -0.91 0.294 0.201 0.193 0.312 -1.85 1.19 -0.99 1.68
65 0.305 0.192 0.195 0.308 0.42 -1.14 -0.27 0.80 0.302 0.195 0.198 0.306 -0.40 -0.25 0.09 0.54
66 0.303 0.193 0.205 0.299 -0.02 -0.83 1.89 -0.90 0.301 0.203 0.201 0.294 -0.54 1.72 0.85 -1.67
67 0.307 0.200 0.193 0.301 0.75 0.67 -0.82 -0.62 0.306 0.189 0.201 0.304 0.44 -1.50 0.84 0.14
68 0.301 0.199 0.196 0.304 -0.29 0.54 -0.22 0.02 0.311 0.196 0.190 0.304 1.32 0.03 -1.79 0.21
69 0.308 0.191 0.203 0.298 0.93 -1.33 1.57 -1.14 0.308 0.198 0.200 0.294 0.78 0.53 0.62 -1.77
70 0.299 0.206 0.187 0.308 -0.75 2.04 -2.07 0.77 0.303 0.196 0.198 0.303 -0.15 0.07 0.16 -0.04
71 0.310 0.196 0.191 0.304 1.28 -0.25 -1.32 0.08 0.305 0.191 0.193 0.311 0.25 -1.10 -0.98 1.55
72 0.302 0.196 0.194 0.307 -0.10 -0.13 -0.53 0.67 0.307 0.188 0.195 0.311 0.57 -1.75 -0.69 1.54
73 0.306 0.192 0.193 0.308 0.62 -1.07 -0.68 0.89 0.302 0.198 0.205 0.295 -0.29 0.42 1.61 -1.47
74 0.311 0.200 0.188 0.301 1.59 0.65 -1.90 -0.51 0.299 0.195 0.197 0.309 -0.88 -0.12 -0.11 1.08
75 0.300 0.195 0.203 0.302 -0.47 -0.51 1.47 -0.36 0.298 0.193 0.202 0.307 -1.14 -0.54 1.05 0.70
76 0.307 0.200 0.199 0.294 0.77 0.77 0.51 -1.88 0.300 0.197 0.196 0.308 -0.76 0.20 -0.36 0.90
77 0.309 0.200 0.189 0.301 1.26 0.77 -1.70 -0.46 0.304 0.196 0.189 0.311 0.13 -0.00 -2.02 1.63
78 0.302 0.193 0.203 0.301 -0.21 -0.79 1.55 -0.44 0.312 0.193 0.198 0.298 1.54 -0.61 0.00 -1.02
79 0.301 0.202 0.195 0.302 -0.44 1.24 -0.27 -0.40 0.305 0.199 0.196 0.300 0.26 0.69 -0.43 -0.49
80 0.302 0.199 0.198 0.301 -0.26 0.51 0.28 -0.43 0.311 0.199 0.192 0.297 1.48 0.84 -1.18 -1.18
81 0.313 0.198 0.197 0.292 1.87 0.34 0.16 -2.30 0.313 0.197 0.189 0.301 1.79 0.36 -1.89 -0.46
82 0.298 0.204 0.195 0.303 -1.00 1.57 -0.28 -0.12 0.299 0.193 0.209 0.299 -0.95 -0.63 2.53 -0.70
83 0.308 0.189 0.197 0.307 0.91 -1.85 0.02 0.68 0.308 0.192 0.195 0.305 0.90 -0.87 -0.57 0.34
84 0.312 0.187 0.196 0.305 1.68 -2.13 -0.00 0.17 0.311 0.190 0.197 0.303 1.30 -1.29 -0.22 0.01
85 0.299 0.198 0.197 0.306 -0.81 0.28 0.07 0.51 0.303 0.201 0.200 0.296 -0.08 1.21 0.51 -1.40
86 0.306 0.197 0.188 0.309 0.59 0.05 -1.95 1.05 0.301 0.205 0.195 0.298 -0.48 2.15 -0.57 -0.88
87 0.308 0.187 0.198 0.307 0.97 -2.17 0.41 0.55 0.309 0.193 0.202 0.296 0.99 -0.56 1.02 -1.39
88 0.304 0.195 0.197 0.304 0.15 -0.38 0.10 0.10 0.307 0.194 0.196 0.303 0.63 -0.34 -0.43 0.03
89 0.310 0.196 0.194 0.301 1.29 -0.26 -0.58 -0.57 0.309 0.193 0.194 0.304 0.96 -0.53 -0.85 0.23
90 0.305 0.200 0.201 0.293 0.49 0.79 0.93 -1.98 0.309 0.198 0.190 0.304 0.95 0.41 -1.65 0.12
91 0.305 0.202 0.196 0.296 0.45 1.17 -0.01 -1.45 0.300 0.192 0.207 0.301 -0.69 -0.87 2.09 -0.36
92 0.310 0.191 0.194 0.304 1.33 -1.20 -0.46 0.11 0.309 0.195 0.198 0.299 1.00 -0.19 0.01 -0.84
93 0.297 0.188 0.204 0.311 -1.07 -1.90 1.61 1.32 0.302 0.203 0.198 0.297 -0.36 1.55 0.08 -1.04
94 0.298 0.199 0.201 0.303 -0.95 0.40 0.96 -0.23 0.303 0.206 0.197 0.294 -0.09 2.27 -0.08 -1.81
95 0.298 0.196 0.190 0.317 -1.00 -0.25 -1.54 2.56 0.296 0.190 0.198 0.316 -1.55 -1.21 0.16 2.46
96 0.306 0.190 0.196 0.308 0.52 -1.49 -0.01 0.77 0.289 0.194 0.207 0.310 -2.76 -0.45 2.00 1.41
97 0.297 0.196 0.199 0.308 -1.08 -0.21 0.56 0.78 0.298 0.196 0.199 0.308 -1.04 0.00 0.20 0.87
98 0.303 0.204 0.192 0.301 -0.06 1.57 -0.93 -0.49 0.305 0.196 0.192 0.307 0.26 0.11 -1.24 0.72
99 0.309 0.194 0.197 0.300 1.11 -0.71 0.20 -0.67 0.301 0.191 0.202 0.306 -0.54 -1.03 0.89 0.65
100 0.303 0.204 0.195 0.298 0.03 1.50 -0.22 -1.14 0.310 0.196 0.202 0.292 1.28 0.01 0.87 -2.04
101 0.300 0.201 0.196 0.303 -0.54 0.91 -0.13 -0.14 0.300 0.198 0.196 0.306 -0.63 0.43 -0.32 0.54
102 0.308 0.195 0.206 0.291 0.89 -0.48 2.16 -2.34 0.302 0.193 0.206 0.300 -0.27 -0.68 1.73 -0.64
103 0.303 0.195 0.190 0.311 0.05 -0.35 -1.36 1.42 0.298 0.196 0.207 0.299 -1.08 0.08 2.04 -0.75
104 0.308 0.194 0.191 0.307 1.02 -0.72 -1.21 0.65 0.300 0.192 0.201 0.307 -0.67 -0.86 0.66 0.83
105 0.306 0.200 0.190 0.304 0.57 0.61 -1.36 0.08 0.300 0.187 0.202 0.311 -0.79 -1.84 0.98 1.53
106 0.304 0.200 0.197 0.299 0.22 0.58 0.15 -0.85 0.299 0.204 0.194 0.303 -0.87 1.80 -0.80 0.01
107 0.304 0.200 0.193 0.304 0.16 0.64 -0.84 0.01 0.294 0.192 0.201 0.313 -1.77 -0.81 0.63 1.92
108 0.305 0.190 0.199 0.306 0.44 -1.53 0.59 0.37 0.299 0.205 0.197 0.299 -1.00 2.02 -0.09 -0.67
109 0.300 0.197 0.198 0.305 -0.58 -0.03 0.36 0.30 0.308 0.191 0.205 0.296 0.75 -0.99 1.67 -1.34
110 0.307 0.205 0.183 0.305 0.73 1.89 -3.03 0.26 0.306 0.201 0.185 0.307 0.49 1.25 -2.75 0.82
111 0.306 0.188 0.200 0.306 0.67 -2.05 0.81 0.40 0.312 0.189 0.193 0.306 1.60 -1.50 -0.92 0.49
112 0.305 0.202 0.194 0.299 0.47 1.11 -0.58 -0.93 0.309 0.197 0.190 0.304 0.95 0.35 -1.63 0.16
113 0.298 0.191 0.197 0.313 -0.87 -1.22 0.16 1.79 0.312 0.181 0.199 0.308 1.52 -3.23 0.25 1.05
114 0.309 0.186 0.198 0.307 1.13 -2.41 0.39 0.63 0.309 0.191 0.193 0.307 0.93 -0.97 -0.90 0.69
115 0.301 0.201 0.203 0.296 -0.38 0.89 1.32 -1.53 0.291 0.201 0.201 0.307 -2.48 1.27 0.78 0.71
116 0.309 0.200 0.195 0.295 1.22 0.72 -0.25 -1.63 0.303 0.198 0.190 0.309 -0.07 0.51 -1.68 1.09
117 0.308 0.186 0.200 0.306 1.01 -2.45 0.84 0.38 0.300 0.191 0.207 0.302 -0.64 -1.04 2.07 -0.26
118 0.305 0.193 0.203 0.299 0.31 -0.90 1.51 -0.84 0.302 0.194 0.205 0.300 -0.33 -0.47 1.57 -0.63
119 0.295 0.200 0.192 0.313 -1.51 0.76 -1.05 1.76 0.299 0.194 0.196 0.310 -0.82 -0.29 -0.28 1.32
120 0.304 0.205 0.190 0.301 0.16 1.80 -1.43 -0.49 0.300 0.195 0.202 0.303 -0.73 -0.09 0.85 0.08
121 0.303 0.201 0.197 0.298 0.09 1.00 0.18 -1.11 0.306 0.198 0.203 0.293 0.43 0.59 1.18 -1.96
122 0.299 0.204 0.185 0.311 -0.70 1.57 -2.45 1.46 0.302 0.196 0.198 0.304 -0.29 0.02 0.18 0.12
123 0.304 0.198 0.202 0.296 0.15 0.27 1.13 -1.36 0.313 0.191 0.199 0.297 1.82 -0.93 0.21 -1.20
124 0.301 0.199 0.204 0.296 -0.41 0.41 1.65 -1.37 0.300 0.203 0.194 0.303 -0.76 1.69 -0.78 -0.02
125 0.303 0.192 0.196 0.309 0.03 -1.13 -0.03 0.97 0.320 0.187 0.192 0.301 3.11 -1.92 -1.27 -0.36
126 0.304 0.200 0.196 0.300 0.21 0.70 -0.15 -0.69 0.312 0.189 0.192 0.308 1.46 -1.55 -1.25 0.96
127 0.296 0.204 0.197 0.302 -1.26 1.62 0.19 -0.30 0.301 0.194 0.197 0.309 -0.55 -0.46 -0.21 1.12
128 0.312 0.186 0.203 0.299 1.71 -2.28 1.30 -0.86 0.306 0.190 0.193 0.311 0.36 -1.27 -0.93 1.55
129 0.303 0.195 0.198 0.305 0.00 -0.50 0.31 0.16 0.301 0.192 0.206 0.301 -0.53 -0.77 1.82 -0.38
130 0.299 0.193 0.192 0.316 -0.73 -0.80 -0.92 2.21 0.312 0.191 0.202 0.295 1.52 -0.95 0.97 -1.54
131 0.305 0.195 0.191 0.308 0.47 -0.37 -1.11 0.81 0.299 0.197 0.201 0.303 -0.87 0.18 0.72 0.09
132 0.298 0.201 0.202 0.298 -0.86 0.86 1.26 -0.98 0.306 0.201 0.197 0.297 0.43 1.06 -0.24 -1.14
133 0.299 0.202 0.202 0.296 -0.72 1.14 1.25 -1.35 0.302 0.195 0.198 0.305 -0.24 -0.07 0.01 0.29
134 0.298 0.201 0.188 0.313 -0.86 0.85 -1.92 1.78 0.306 0.203 0.190 0.301 0.50 1.61 -1.74 -0.38
135 0.309 0.193 0.198 0.300 1.13 -0.77 0.33 -0.74 0.308 0.188 0.197 0.307 0.76 -1.58 -0.18 0.76
136 0.305 0.201 0.197 0.297 0.45 0.88 0.02 -1.22 0.312 0.196 0.196 0.296 1.50 0.14 -0.43 -1.24
137 0.306 0.197 0.189 0.307 0.57 0.12 -1.56 0.68 0.308 0.201 0.186 0.305 0.76 1.20 -2.42 0.30
138 0.301 0.197 0.199 0.303 -0.43 0.12 0.58 -0.17 0.299 0.196 0.198 0.307 -0.93 0.11 0.19 0.67
139 0.296 0.201 0.201 0.302 -1.34 0.83 0.99 -0.24 0.302 0.201 0.192 0.305 -0.36 1.14 -1.25 0.46
140 0.306 0.186 0.198 0.311 0.50 -2.39 0.22 1.38 0.300 0.197 0.195 0.308 -0.70 0.33 -0.53 0.87
141 0.295 0.192 0.204 0.309 -1.56 -1.00 1.69 0.96 0.302 0.193 0.200 0.306 -0.40 -0.58 0.44 0.51
142 0.289 0.203 0.197 0.311 -2.68 1.27 0.20 1.41 0.296 0.198 0.201 0.305 -1.49 0.55 0.63 0.47
143 0.309 0.194 0.195 0.302 1.04 -0.64 -0.30 -0.23 0.301 0.192 0.196 0.311 -0.48 -0.68 -0.40 1.41
144 0.305 0.187 0.194 0.313 0.39 -2.03 -0.50 1.80 0.299 0.194 0.206 0.301 -0.87 -0.45 1.90 -0.39
145 0.304 0.199 0.191 0.306 0.16 0.55 -1.20 0.40 0.306 0.203 0.190 0.302 0.33 1.55 -1.71 -0.18
146 0.299 0.205 0.194 0.302 -0.74 1.70 -0.53 -0.27 0.302 0.198 0.192 0.308 -0.37 0.59 -1.24 0.93
147 0.304 0.190 0.194 0.312 0.22 -1.56 -0.50 1.56 0.308 0.194 0.197 0.301 0.80 -0.42 -0.19 -0.27
148 0.306 0.196 0.195 0.303 0.48 -0.10 -0.34 -0.10 0.296 0.203 0.200 0.301 -1.39 1.61 0.41 -0.35
149 0.301 0.198 0.198 0.304 -0.43 0.14 0.31 0.04 0.310 0.192 0.191 0.307 1.07 -0.73 -1.32 0.71
150 0.299 0.200 0.196 0.305 -0.81 0.77 -0.13 0.25 0.301 0.194 0.204 0.301 -0.48 -0.43 1.34 -0.31
151 0.300 0.200 0.195 0.305 -0.59 0.61 -0.26 0.29 0.311 0.197 0.198 0.294 1.29 0.30 0.03 -1.58
152 0.302 0.200 0.196 0.303 -0.22 0.69 -0.18 -0.21 0.304 0.195 0.197 0.304 0.08 -0.19 -0.16 0.23
153 0.313 0.193 0.196 0.298 1.89 -0.88 -0.09 -1.05 0.301 0.195 0.202 0.302 -0.56 -0.15 0.97 -0.15
154 0.305 0.190 0.193 0.312 0.44 -1.50 -0.77 1.52 0.304 0.198 0.199 0.299 0.02 0.56 0.35 -0.81
155 0.301 0.202 0.193 0.304 -0.32 1.06 -0.76 0.05 0.303 0.207 0.189 0.301 -0.15 2.33 -1.74 -0.36
156 0.309 0.192 0.198 0.301 1.05 -1.04 0.32 -0.43 0.306 0.193 0.200 0.301 0.40 -0.58 0.55 -0.37
157 0.297 0.201 0.200 0.302 -1.11 0.83 0.76 -0.26 0.304 0.193 0.200 0.303 0.11 -0.65 0.60 -0.07
158 0.306 0.199 0.192 0.303 0.64 0.38 -1.00 -0.10 0.296 0.200 0.193 0.311 -1.48 0.87 -0.97 1.57
159 0.312 0.188 0.199 0.301 1.60 -1.82 0.64 -0.58 0.299 0.200 0.202 0.298 -0.83 0.93 0.98 -0.83
160 0.296 0.202 0.200 0.302 -1.27 1.05 0.77 -0.30 0.307 0.200 0.198 0.296 0.50 0.82 0.10 -1.30
161 0.303 0.192 0.191 0.315 0.05 -1.13 -1.26 2.02 0.306 0.200 0.194 0.300 0.35 0.85 -0.70 -0.48
162 0.308 0.190 0.197 0.304 0.98 -1.39 0.21 0.04 0.299 0.194 0.201 0.306 -0.88 -0.30 0.62 0.61
163 0.296 0.201 0.199 0.304 -1.31 0.97 0.46 0.07 0.316 0.192 0.199 0.293 2.27 -0.81 0.33 -1.86
164 0.307 0.198 0.189 0.307 0.71 0.19 -1.61 0.52 0.313 0.189 0.191 0.306 1.74 -1.38 -1.38 0.65
165 0.301 0.200 0.199 0.300 -0.32 0.69 0.44 -0.65 0.313 0.200 0.196 0.292 1.66 0.85 -0.38 -2.07
166 0.307 0.202 0.198 0.293 0.71 1.12 0.36 -1.99 0.302 0.200 0.199 0.300 -0.38 0.86 0.29 -0.62
167 0.297 0.199 0.195 0.310 -1.16 0.45 -0.40 1.11 0.303 0.197 0.192 0.308 -0.18 0.34 -1.12 0.85
168 0.306 0.190 0.201 0.303 0.65 -1.49 0.99 -0.22 0.296 0.199 0.204 0.301 -1.45 0.74 1.40 -0.40
169 0.299 0.202 0.201 0.298 -0.78 1.16 0.97 -1.05 0.298 0.201 0.201 0.300 -1.05 1.20 0.63 -0.54
170 0.308 0.192 0.191 0.309 0.85 -1.06 -1.13 1.04 0.305 0.190 0.196 0.309 0.25 -1.11 -0.45 1.10
171 0.304 0.189 0.204 0.303 0.23 -1.71 1.49 -0.04 0.292 0.203 0.208 0.296 -2.16 1.63 2.29 -1.23
172 0.301 0.207 0.195 0.296 -0.34 2.17 -0.24 -1.33 0.312 0.196 0.200 0.292 1.47 0.11 0.51 -2.01
173 0.302 0.196 0.195 0.307 -0.12 -0.10 -0.42 0.57 0.313 0.192 0.189 0.306 1.77 -0.88 -1.74 0.49
174 0.309 0.192 0.191 0.307 1.17 -1.04 -1.06 0.64 0.300 0.193 0.199 0.308 -0.74 -0.65 0.35 1.00
175 0.305 0.202 0.187 0.306 0.43 1.05 -2.09 0.47 0.299 0.202 0.198 0.300 -0.83 1.41 0.10 -0.47
176 0.306 0.200 0.194 0.300 0.56 0.57 -0.43 -0.68 0.298 0.200 0.193 0.308 -0.97 0.97 -0.98 0.98
177 0.299 0.200 0.206 0.295 -0.63 0.56 2.08 -1.66 0.299 0.203 0.192 0.306 -0.88 1.56 -1.17 0.55
178 0.295 0.200 0.198 0.307 -1.45 0.63 0.41 0.55 0.299 0.207 0.201 0.293 -0.86 2.36 0.79 -1.86
179 0.301 0.204 0.195 0.300 -0.29 1.42 -0.40 -0.59 0.305 0.193 0.192 0.310 0.19 -0.57 -1.13 1.28
180 0.297 0.193 0.198 0.312 -1.16 -0.83 0.40 1.54 0.308 0.198 0.201 0.293 0.83 0.44 0.72 -1.84
181 0.300 0.199 0.199 0.302 -0.58 0.52 0.47 -0.28 0.305 0.192 0.199 0.304 0.14 -0.75 0.35 0.21
182 0.305 0.196 0.196 0.302 0.46 -0.12 -0.02 -0.34 0.309 0.193 0.192 0.307 0.91 -0.64 -1.17 0.65
183 0.296 0.189 0.205 0.310 -1.24 -1.68 1.69 1.23 0.311 0.185 0.199 0.305 1.31 -2.14 0.24 0.33
184 0.292 0.198 0.208 0.301 -1.96 0.28 2.50 -0.45 0.297 0.197 0.206 0.301 -1.26 0.22 1.65 -0.35
185 0.298 0.194 0.194 0.314 -0.86 -0.52 -0.60 1.83 0.302 0.193 0.197 0.308 -0.32 -0.50 -0.12 0.85
186 0.307 0.193 0.202 0.299 0.68 -0.83 1.12 -0.93 0.312 0.195 0.204 0.289 1.53 -0.05 1.23 -2.56
187 0.296 0.199 0.194 0.312 -1.27 0.45 -0.62 1.41 0.302 0.203 0.209 0.286 -0.39 1.54 2.43 -3.05
188 0.296 0.198 0.194 0.312 -1.20 0.13 -0.55 1.56 0.307 0.198 0.190 0.305 0.65 0.43 -1.64 0.40
189 0.307 0.200 0.198 0.295 0.78 0.65 0.33 -1.62 0.311 0.197 0.208 0.283 1.39 0.31 2.26 -3.62
190 0.304 0.200 0.202 0.294 0.17 0.62 1.12 -1.68 0.305 0.195 0.203 0.297 0.18 -0.05 1.03 -1.03
191 0.305 0.202 0.199 0.293 0.43 1.15 0.53 -1.89 0.307 0.195 0.190 0.308 0.57 -0.23 -1.54 0.96
192 0.309 0.193 0.201 0.298 1.03 -0.83 0.86 -1.05 0.309 0.191 0.205 0.295 0.85 -1.02 1.63 -1.39
193 0.306 0.196 0.198 0.300 0.60 -0.16 0.34 -0.75 0.307 0.198 0.193 0.302 0.55 0.59 -0.97 -0.22
194 0.302 0.205 0.187 0.305 -0.10 1.76 -1.97 0.28 0.308 0.191 0.197 0.304 0.83 -0.96 -0.15 0.13
195 0.299 0.199 0.197 0.305 -0.63 0.34 0.11 0.23 0.307 0.203 0.196 0.294 0.62 1.43 -0.32 -1.57
196 0.298 0.201 0.199 0.302 -0.80 0.89 0.48 -0.38 0.308 0.203 0.195 0.294 0.74 1.52 -0.47 -1.64
197 0.310 0.198 0.188 0.304 1.20 0.30 -1.85 0.14 0.301 0.193 0.195 0.311 -0.58 -0.57 -0.52 1.52
198 0.293 0.190 0.204 0.313 -1.77 -1.46 1.54 1.70 0.304 0.194 0.203 0.300 0.06 -0.44 1.04 -0.58
199 0.314 0.196 0.194 0.295 2.07 -0.11 -0.49 -1.55 0.301 0.199 0.196 0.303 -0.43 0.68 -0.28 0.09
200 0.306 0.202 0.193 0.299 0.53 1.14 -0.82 -0.80 0.295 0.198 0.195 0.311 -1.52 0.55 -0.53 1.50
201 0.305 0.197 0.192 0.306 0.30 0.01 -0.86 0.42 0.290 0.196 0.208 0.306 -2.51 0.03 2.15 0.62
202 0.292 0.198 0.197 0.313 -1.98 0.30 0.10 1.63 0.301 0.199 0.201 0.299 -0.54 0.78 0.67 -0.71
203 0.300 0.199 0.191 0.310 -0.50 0.46 -1.11 1.05 0.304 0.196 0.192 0.308 0.01 0.12 -1.14 0.88
204 0.299 0.201 0.200 0.300 -0.71 0.83 0.69 -0.61 0.297 0.192 0.194 0.316 -1.15 -0.68 -0.69 2.33
205 0.298 0.198 0.196 0.308 -0.83 0.16 -0.19 0.86 0.301 0.199 0.204 0.296 -0.53 0.60 1.36 -1.17
206 0.302 0.192 0.185 0.321 -0.14 -1.04 -2.34 3.06 0.304 0.192 0.195 0.309 0.08 -0.75 -0.63 1.10
207 0.303 0.202 0.197 0.297 0.02 1.12 0.18 -1.15 0.304 0.193 0.202 0.300 0.01 -0.46 0.97 -0.45
208 0.298 0.187 0.205 0.309 -0.82 -1.94 1.78 0.95 0.296 0.195 0.195 0.315 -1.36 -0.21 -0.64 2.10
209 0.295 0.207 0.194 0.303 -1.43 2.19 -0.46 -0.06 0.301 0.197 0.190 0.311 -0.52 0.37 -1.51 1.51
210 0.309 0.190 0.194 0.307 1.04 -1.35 -0.51 0.58 0.302 0.192 0.213 0.293 -0.28 -0.82 3.10 -1.70
211 0.304 0.196 0.194 0.306 0.15 -0.13 -0.55 0.44 0.306 0.196 0.193 0.306 0.35 0.00 -1.04 0.55
212 0.305 0.197 0.195 0.303 0.32 0.12 -0.39 -0.09 0.294 0.198 0.197 0.312 -1.82 0.48 -0.17 1.56
213 0.296 0.202 0.209 0.293 -1.25 1.01 2.63 -1.90 0.306 0.194 0.210 0.290 0.39 -0.40 2.51 -2.23
214 0.290 0.204 0.201 0.306 -2.31 1.37 0.92 0.32 0.301 0.194 0.210 0.295 -0.58 -0.26 2.50 -1.36
215 0.299 0.201 0.191 0.308 -0.64 0.91 -1.14 0.84 0.300 0.201 0.202 0.297 -0.67 1.00 0.94 -1.01
216 0.297 0.192 0.198 0.312 -0.98 -0.90 0.40 1.41 0.298 0.207 0.200 0.294 -0.98 2.43 0.48 -1.53
217 0.292 0.202 0.199 0.307 -1.87 1.09 0.43 0.55 0.307 0.198 0.194 0.301 0.63 0.45 -0.79 -0.34
218 0.283 0.201 0.204 0.313 -3.55 0.77 1.46 1.62 0.305 0.197 0.201 0.297 0.20 0.33 0.76 -1.14
219 0.298 0.193 0.202 0.307 -0.85 -0.81 1.08 0.61 0.300 0.196 0.209 0.295 -0.64 0.14 2.23 -1.41
220 0.290 0.206 0.200 0.305 -2.33 1.81 0.68 0.18 0.299 0.195 0.199 0.307 -0.88 -0.09 0.34 0.66
221 0.291 0.200 0.197 0.313 -2.12 0.54 0.06 1.59 0.297 0.193 0.205 0.305 -1.17 -0.55 1.48 0.36
222 0.299 0.192 0.200 0.309 -0.61 -0.95 0.67 0.86 0.309 0.189 0.204 0.299 0.84 -1.37 1.23 -0.72
223 0.298 0.197 0.196 0.308 -0.86 0.13 -0.08 0.82 0.300 0.198 0.203 0.300 -0.72 0.39 1.13 -0.60
224 0.306 0.201 0.192 0.302 0.47 0.78 -0.94 -0.34 0.302 0.205 0.188 0.305 -0.32 1.95 -1.95 0.33
225 0.299 0.193 0.200 0.309 -0.76 -0.85 0.71 0.87 0.295 0.200 0.207 0.297 -1.55 0.94 1.99 -0.98
226 0.302 0.199 0.190 0.309 -0.19 0.42 -1.27 0.93 0.300 0.202 0.201 0.298 -0.75 1.22 0.63 -0.84
227 0.311 0.191 0.190 0.308 1.35 -1.14 -1.35 0.80 0.299 0.203 0.192 0.306 -0.77 1.48 -1.16 0.50
228 0.315 0.184 0.201 0.300 2.11 -2.61 0.97 -0.69 0.296 0.192 0.211 0.300 -1.35 -0.68 2.82 -0.51
229 0.301 0.197 0.205 0.298 -0.42 -0.06 1.73 -1.02 0.308 0.205 0.189 0.298 0.79 1.91 -1.81 -0.87
230 0.295 0.199 0.200 0.307 -1.36 0.36 0.62 0.52 0.306 0.197 0.196 0.301 0.39 0.21 -0.28 -0.33
231 0.301 0.189 0.196 0.314 -0.34 -1.51 -0.19 1.81 0.297 0.198 0.200 0.306 -1.22 0.49 0.40 0.45
232 0.298 0.198 0.193 0.312 -0.93 0.13 -0.73 1.44 0.305 0.196 0.200 0.299 0.15 0.14 0.46 -0.67
233 0.303 0.205 0.191 0.300 0.08 1.69 -1.10 -0.59 0.303 0.194 0.191 0.312 -0.11 -0.39 -1.38 1.65
234 0.304 0.191 0.200 0.305 0.17 -1.12 0.76 0.14 0.312 0.191 0.204 0.293 1.48 -0.99 1.22 -1.68
235 0.305 0.201 0.194 0.299 0.35 0.93 -0.45 -0.76 0.303 0.196 0.195 0.306 -0.07 0.15 -0.62 0.48
236 0.294 0.200 0.194 0.312 -1.53 0.59 -0.55 1.49 0.301 0.202 0.187 0.310 -0.46 1.37 -2.17 1.16
237 0.293 0.190 0.204 0.313 -1.78 -1.37 1.57 1.61 0.302 0.203 0.203 0.293 -0.39 1.49 1.06 -1.81
238 0.303 0.196 0.196 0.304 0.08 -0.11 -0.07 0.08 0.307 0.203 0.199 0.291 0.50 1.51 0.33 -2.08
239 0.305 0.199 0.195 0.300 0.41 0.52 -0.31 -0.59 0.295 0.197 0.200 0.307 -1.49 0.33 0.52 0.75
240 0.306 0.190 0.196 0.308 0.60 -1.44 -0.03 0.68 0.302 0.205 0.202 0.291 -0.36 1.93 0.87 -2.05
241 0.302 0.196 0.188 0.314 -0.11 -0.27 -1.63 1.75 0.303 0.206 0.200 0.291 -0.18 2.14 0.42 -2.03
242 0.296 0.201 0.193 0.310 -1.16 0.84 -0.81 1.13 0.296 0.204 0.194 0.306 -1.34 1.62 -0.69 0.54
243 0.301 0.194 0.203 0.302 -0.40 -0.56 1.30 -0.24 0.301 0.200 0.203 0.295 -0.45 0.96 1.18 -1.41
244 0.304 0.207 0.189 0.300 0.23 2.00 -1.59 -0.59 0.298 0.199 0.194 0.309 -0.93 0.65 -0.72 1.00
245 0.299 0.201 0.205 0.295 -0.69 0.84 1.76 -1.56 0.305 0.194 0.188 0.313 0.29 -0.43 -1.98 1.79
246 0.305 0.188 0.204 0.304 0.29 -1.81 1.43 0.04 0.301 0.192 0.211 0.296 -0.53 -0.75 2.67 -1.14
247 0.303 0.203 0.200 0.294 -0.01 1.15 0.72 -1.61 0.305 0.195 0.207 0.292 0.26 -0.04 1.97 -1.93
248 0.307 0.186 0.209 0.298 0.65 -2.11 2.52 -1.01 0.296 0.204 0.201 0.299 -1.30 1.74 0.65 -0.75
249 0.316 0.184 0.203 0.297 2.25 -2.67 1.33 -1.09 0.311 0.187 0.209 0.293 1.33 -1.71 2.20 -1.76
250 0.329 0.156 0.205 0.310 4.31 -7.88 1.66 1.07 0.318 0.158 0.221 0.303 2.39 -7.38 4.56 0.02
251 0.411 0 0.249 0.339 16.15 . 9.08 5.31 0.392 0 0.274 0.334 13.09 . 13.12 4.68
=== final results =========================================================================================================================================
[ fA fC fG fT R1] [ fA fC fG fT R2]
median 0.304 0.196 0.197 0.303 0.303 0.197 0.198 0.303
pos fA fC fG fT zA zC zG zT nbases fA fC fG fT zA zC zG zT nbases
1 0.303 0.198 0.200 0.299 -0.33 1.15 1.62 -2.06 61099 0.303 0.199 0.200 0.299 -0.22 1.25 1.63 -2.28 60307
2 0.302 0.197 0.196 0.306 -1.11 0.40 -1.02 1.65 74687 0.301 0.197 0.196 0.305 -1.10 0.38 -0.95 1.59 73973
3 0.303 0.197 0.195 0.305 -0.49 0.53 -1.19 1.06 78104 0.303 0.197 0.196 0.304 -0.27 0.31 -1.15 0.99 78037
4 0.301 0.199 0.195 0.306 -1.52 1.64 -1.74 1.62 78258 0.300 0.199 0.195 0.305 -1.68 1.62 -1.56 1.63 78200
5 0.306 0.194 0.195 0.305 1.62 -1.83 -1.93 1.63 96179 0.305 0.194 0.195 0.305 1.64 -1.92 -1.87 1.64 96339
6 0.306 0.194 0.199 0.301 1.64 -2.24 1.64 -1.13 113389 0.305 0.194 0.199 0.301 1.63 -2.10 1.63 -1.23 113290
7 0.305 0.195 0.198 0.302 1.30 -1.41 0.95 -0.91 148439 0.304 0.195 0.199 0.302 1.12 -1.49 1.38 -1.02 147285
8 0.306 0.195 0.197 0.302 1.60 -0.78 0.22 -1.11 148920 0.305 0.196 0.198 0.301 1.62 -0.74 0.20 -1.16 147653
9 0.300 0.198 0.199 0.303 -2.80 1.62 1.63 -0.01 156454 0.300 0.198 0.199 0.303 -2.80 1.63 1.61 -0.00 155790
10 0.303 0.197 0.195 0.305 -0.59 1.17 -2.24 1.51 187364 0.303 0.198 0.195 0.304 -0.46 1.22 -2.31 1.41 186685
11 0.304 0.194 0.199 0.303 0.62 -2.29 1.64 -0.06 188588 0.304 0.195 0.199 0.302 0.57 -1.79 1.64 -0.45 188460
12 0.302 0.197 0.196 0.305 -1.60 0.78 -0.74 1.57 213519 0.302 0.197 0.197 0.304 -1.29 0.55 -0.74 1.46 212780
13 0.301 0.198 0.198 0.303 -2.90 1.62 1.59 0.13 214292 0.300 0.198 0.199 0.303 -2.82 1.34 1.63 0.25 213910
14 0.301 0.198 0.197 0.305 -3.26 1.64 0.31 1.58 220631 0.300 0.198 0.198 0.304 -3.32 1.65 0.38 1.56 220360
15 0.305 0.198 0.198 0.299 1.61 1.58 1.63 -4.39 247483 0.304 0.198 0.199 0.299 1.60 1.63 1.57 -4.36 246732
16 0.305 0.196 0.198 0.301 1.51 0.25 1.16 -2.73 265737 0.304 0.197 0.198 0.300 1.57 0.40 1.21 -2.96 264585
17 0.305 0.195 0.195 0.304 1.64 -1.56 -2.15 1.57 270812 0.304 0.196 0.196 0.304 1.58 -1.52 -1.99 1.46 269506
18 0.305 0.195 0.195 0.304 1.63 -1.61 -2.15 1.62 275195 0.304 0.195 0.196 0.304 1.64 -2.39 -1.40 1.64 273228
19 0.305 0.194 0.197 0.304 1.64 -2.90 -0.49 1.29 283343 0.304 0.194 0.197 0.304 1.64 -3.70 -0.06 1.61 280895
20 0.305 0.195 0.196 0.304 1.57 -2.06 -1.60 1.60 286013 0.304 0.195 0.197 0.304 1.44 -1.98 -1.23 1.34 285312
21 0.305 0.195 0.196 0.304 1.04 -1.51 -1.02 1.15 289342 0.304 0.196 0.197 0.303 1.57 -1.36 -1.05 0.52 287834
22 0.303 0.196 0.197 0.304 -0.63 -0.15 0.16 0.62 289732 0.303 0.197 0.197 0.303 -0.29 0.95 -1.35 0.63 288120
23 0.303 0.196 0.197 0.304 -0.66 -0.46 -0.53 1.51 289743 0.303 0.197 0.197 0.304 -0.24 -0.24 -0.79 1.13 288120
24 0.304 0.195 0.197 0.303 0.79 -1.09 0.01 0.14 289934 0.304 0.196 0.197 0.303 1.51 -1.58 -0.27 0.09 288160
25 0.303 0.196 0.197 0.304 -0.39 0.19 -0.56 0.71 289940 0.302 0.196 0.198 0.304 -0.89 -0.46 0.27 1.05 288160
26 0.305 0.196 0.195 0.304 1.18 0.02 -2.26 0.75 289983 0.303 0.196 0.196 0.304 0.53 -0.89 -1.38 1.43 288160
27 0.305 0.195 0.198 0.303 1.17 -1.94 1.13 -0.47 290041 0.304 0.195 0.198 0.303 1.19 -2.15 0.44 0.29 288292
28 0.304 0.196 0.197 0.303 0.53 -0.35 -0.30 0.03 290367 0.304 0.196 0.197 0.303 1.07 -1.33 -0.45 0.47 289165
29 0.305 0.196 0.196 0.303 1.51 -0.68 -1.72 0.57 290378 0.304 0.196 0.197 0.304 1.59 -1.57 -1.25 0.85 289165
30 0.305 0.196 0.197 0.303 1.22 -0.53 -0.29 -0.52 290710 0.304 0.196 0.197 0.303 0.82 -0.70 -0.58 0.29 289303
31 0.304 0.196 0.197 0.303 0.06 -0.36 0.57 -0.25 291233 0.303 0.196 0.198 0.303 0.08 -0.37 0.51 -0.21 289868
32 0.304 0.196 0.196 0.303 0.74 -0.04 -0.77 -0.03 291245 0.304 0.197 0.197 0.303 1.25 0.11 -1.36 -0.17 289868
33 0.304 0.196 0.197 0.302 0.43 -0.18 0.40 -0.62 291245 0.304 0.196 0.197 0.303 1.20 -1.18 -0.44 0.21 290054
34 0.303 0.196 0.197 0.304 -0.66 0.16 -0.29 0.77 291296 0.302 0.196 0.197 0.304 -0.85 -0.35 -0.33 1.43 290054
35 0.304 0.196 0.197 0.303 0.83 0.00 -0.50 -0.40 291305 0.303 0.197 0.197 0.303 0.56 -0.14 -0.85 0.30 290173
36 0.304 0.195 0.197 0.303 0.74 -1.18 0.09 0.20 291333 0.304 0.195 0.197 0.304 0.78 -1.68 -0.33 0.97 290179
37 0.304 0.196 0.197 0.303 -0.05 -0.69 0.28 0.40 291380 0.303 0.196 0.198 0.303 -0.10 -0.96 1.09 -0.02 290249
38 0.305 0.196 0.196 0.303 1.36 0.14 -1.62 -0.08 291380 0.304 0.196 0.197 0.303 0.92 -1.01 -0.89 0.72 290268
39 0.305 0.196 0.196 0.303 1.51 -0.43 -1.50 0.16 291599 0.304 0.196 0.197 0.303 1.42 -1.46 -0.30 0.10 290313
40 0.304 0.196 0.197 0.303 0.82 -0.82 0.34 -0.41 292125 0.304 0.196 0.198 0.302 1.27 -0.68 0.42 -1.04 290875
41 0.305 0.196 0.196 0.303 1.52 -0.16 -1.63 0.03 292208 0.304 0.196 0.196 0.303 1.47 -0.49 -1.87 0.58 290966
42 0.305 0.195 0.197 0.303 1.20 -1.30 -0.68 0.51 292504 0.304 0.196 0.197 0.303 1.49 -1.58 -0.72 0.50 291123
43 0.304 0.196 0.198 0.302 -0.08 -0.04 0.91 -0.67 292932 0.303 0.196 0.198 0.302 0.22 -1.02 1.24 -0.42 291588
44 0.304 0.196 0.196 0.303 0.39 0.01 -0.92 0.40 293057 0.304 0.196 0.197 0.303 1.00 -0.77 -0.87 0.42 291836
45 0.304 0.196 0.198 0.302 0.68 -0.89 0.84 -0.64 293111 0.304 0.196 0.197 0.303 1.10 -1.37 -0.70 0.69 291836
46 0.304 0.196 0.196 0.304 -0.22 -0.29 -0.72 1.09 293128 0.303 0.196 0.198 0.303 0.09 -1.23 0.41 0.62 291872
47 0.304 0.197 0.196 0.303 0.56 0.46 -1.70 0.51 293144 0.303 0.198 0.195 0.304 0.32 1.25 -2.83 1.05 291872
48 0.305 0.195 0.197 0.303 1.28 -1.84 0.34 0.01 293209 0.304 0.195 0.198 0.303 1.26 -1.71 0.16 0.07 291892
49 0.304 0.196 0.197 0.303 -0.07 -0.83 0.51 0.34 293493 0.304 0.195 0.198 0.303 0.86 -1.66 0.22 0.39 292163
50 0.305 0.196 0.196 0.303 1.46 -0.85 -1.27 0.38 293507 0.304 0.196 0.196 0.303 1.12 -0.40 -1.53 0.55 292168
51 0.305 0.195 0.197 0.303 1.06 -1.34 0.59 -0.42 293854 0.304 0.196 0.198 0.302 1.51 -1.56 0.65 -0.73 292330
52 0.304 0.196 0.197 0.302 0.44 -0.26 0.60 -0.75 293960 0.303 0.197 0.197 0.303 -0.01 0.56 -0.18 -0.32 292955
53 0.305 0.197 0.196 0.303 0.94 0.61 -1.94 0.21 294046 0.304 0.196 0.196 0.303 1.50 -0.38 -1.49 0.12 292918
54 0.305 0.196 0.197 0.303 0.95 -0.24 -0.23 -0.54 294014 0.304 0.196 0.197 0.303 1.43 -1.39 -0.42 0.14 292934
55 0.303 0.197 0.198 0.303 -0.45 0.36 0.76 -0.52 294051 0.303 0.196 0.198 0.303 0.39 -0.32 0.06 -0.16 292849
56 0.303 0.196 0.197 0.304 -0.70 -0.41 -0.11 1.16 294006 0.303 0.197 0.197 0.304 -0.03 -0.24 -0.84 0.97 292739
57 0.305 0.196 0.197 0.303 0.95 -0.69 -0.56 0.12 294443 0.304 0.196 0.197 0.302 1.47 -0.44 -0.34 -0.80 292626
58 0.304 0.197 0.196 0.303 0.25 0.61 -1.11 0.18 294528 0.302 0.196 0.198 0.304 -0.65 -0.98 0.27 1.27 292630
59 0.305 0.196 0.195 0.304 0.96 -0.15 -2.11 0.99 294364 0.304 0.196 0.196 0.304 0.83 -0.97 -1.85 1.62 292463
60 0.304 0.195 0.197 0.303 0.78 -1.09 -0.12 0.27 294161 0.304 0.195 0.198 0.303 1.49 -2.24 0.29 0.20 292767
61 0.303 0.197 0.197 0.303 -0.45 0.52 -0.12 0.10 293952 0.303 0.196 0.197 0.304 -0.18 -0.50 -0.44 0.99 293036
62 0.305 0.196 0.195 0.304 1.28 -0.23 -2.73 1.28 293802 0.304 0.196 0.195 0.304 1.58 -0.87 -2.75 1.55 292891
63 0.304 0.196 0.197 0.303 0.48 -0.27 -0.38 0.08 293979 0.304 0.196 0.197 0.303 1.58 -1.57 -0.81 0.49 293076
64 0.303 0.197 0.197 0.303 -0.59 0.59 0.65 -0.48 293705 0.303 0.197 0.197 0.303 0.40 0.01 -0.57 0.08 292960
65 0.304 0.196 0.196 0.304 0.34 -0.25 -1.77 1.41 293432 0.303 0.196 0.197 0.304 0.19 -1.07 -0.28 0.99 292700
66 0.305 0.196 0.197 0.302 1.25 -0.14 0.40 -1.48 293322 0.304 0.196 0.197 0.302 1.60 -0.51 -0.91 -0.37 292438
67 0.304 0.197 0.197 0.302 0.04 0.53 0.28 -0.74 293093 0.303 0.197 0.198 0.302 -0.15 0.27 0.45 -0.48 292429
68 0.304 0.197 0.195 0.304 0.43 0.51 -2.59 1.37 292735 0.304 0.196 0.196 0.304 1.46 -0.93 -1.93 1.01 292118
69 0.304 0.196 0.198 0.302 0.34 -0.55 1.36 -1.04 292334 0.304 0.197 0.198 0.302 0.68 0.20 0.43 -1.22 292039
70 0.303 0.196 0.198 0.303 -0.54 -0.09 0.85 -0.12 292028 0.303 0.197 0.197 0.303 0.56 0.24 -0.86 -0.02 291681
71 0.303 0.197 0.196 0.304 -0.32 0.65 -1.41 0.98 291909 0.303 0.197 0.196 0.304 0.25 0.19 -2.17 1.47 291349
72 0.305 0.195 0.197 0.303 1.60 -1.42 -0.48 0.04 291570 0.304 0.195 0.198 0.303 1.46 -2.21 0.08 0.38 291231
73 0.303 0.196 0.197 0.303 -0.42 0.29 -0.30 0.43 291472 0.302 0.197 0.198 0.303 -0.74 0.39 0.26 0.18 290999
74 0.305 0.197 0.195 0.304 0.95 0.43 -2.94 1.22 291116 0.303 0.197 0.196 0.304 0.45 0.71 -2.27 0.90 290745
75 0.304 0.196 0.197 0.302 0.75 -0.78 0.67 -0.65 291002 0.304 0.196 0.198 0.303 0.67 -0.99 0.44 -0.19 290722
76 0.303 0.196 0.197 0.303 -0.37 -0.41 0.54 0.26 290769 0.303 0.196 0.198 0.303 -0.25 -0.78 0.14 0.81 290645
77 0.304 0.196 0.196 0.304 0.51 -0.75 -1.52 1.46 290538 0.303 0.197 0.197 0.303 -0.50 0.41 -0.46 0.54 290305
78 0.304 0.196 0.197 0.303 0.85 -0.34 -0.50 -0.13 290289 0.304 0.196 0.197 0.302 1.24 -0.48 -0.12 -0.71 290029
79 0.304 0.197 0.197 0.302 0.05 0.63 0.55 -1.08 289928 0.303 0.196 0.198 0.303 -0.26 -0.44 0.35 0.34 289981
80 0.305 0.196 0.195 0.304 1.22 0.01 -3.04 1.41 289556 0.304 0.196 0.196 0.304 0.79 -0.60 -1.73 1.23 289580
81 0.304 0.195 0.197 0.303 0.78 -1.35 0.53 -0.07 289327 0.304 0.196 0.197 0.303 0.87 -0.34 -0.41 -0.22 289543
82 0.303 0.197 0.197 0.303 -0.91 1.30 -0.63 0.33 289025 0.302 0.197 0.198 0.303 -0.95 1.02 0.25 -0.15 289528
83 0.304 0.196 0.196 0.304 0.36 -0.12 -1.96 1.45 288778 0.304 0.197 0.195 0.304 0.71 0.11 -2.74 1.57 289285
84 0.305 0.195 0.197 0.304 1.27 -2.03 -0.67 1.06 288418 0.304 0.195 0.198 0.303 0.75 -2.09 0.43 0.68 289517
85 0.303 0.197 0.197 0.303 -0.70 0.63 0.34 -0.13 288134 0.304 0.196 0.198 0.303 1.08 -1.55 0.60 -0.27 289207
86 0.305 0.196 0.195 0.304 1.13 -0.14 -2.07 0.79 287804 0.304 0.196 0.197 0.303 0.63 -0.56 -0.71 0.47 288876
87 0.304 0.197 0.197 0.303 -0.02 0.39 -0.20 -0.14 287560 0.304 0.196 0.198 0.302 0.65 -0.97 0.98 -0.66 288544
88 0.303 0.197 0.198 0.302 -1.11 0.45 1.54 -0.61 287269 0.303 0.196 0.198 0.302 0.40 -0.48 0.54 -0.45 288219
89 0.303 0.197 0.197 0.303 -1.32 1.30 -0.34 0.49 286963 0.303 0.197 0.197 0.303 0.16 0.04 -0.55 0.29 287964
90 0.304 0.196 0.197 0.303 0.51 -1.04 0.10 0.30 286894 0.304 0.195 0.198 0.303 0.84 -1.80 0.87 -0.04 287757
91 0.303 0.197 0.197 0.303 -0.51 0.91 0.12 -0.38 286566 0.303 0.198 0.197 0.302 0.21 1.37 -0.49 -0.98 287639
92 0.304 0.196 0.195 0.304 0.52 0.18 -2.55 1.53 286438 0.304 0.197 0.195 0.304 0.83 0.15 -2.96 1.61 287692
93 0.305 0.195 0.197 0.303 1.42 -1.41 0.07 -0.26 286161 0.304 0.196 0.198 0.302 0.89 -0.32 0.23 -0.81 287724
94 0.304 0.196 0.198 0.302 0.35 -0.16 1.07 -1.13 286011 0.302 0.196 0.198 0.303 -0.64 -0.34 0.47 0.53 287410
95 0.304 0.196 0.196 0.304 0.49 -0.66 -1.76 1.60 285816 0.303 0.198 0.196 0.303 0.49 1.35 -2.25 0.29 287160
96 0.305 0.195 0.197 0.303 1.55 -2.09 0.05 0.21 285737 0.304 0.196 0.198 0.303 1.13 -1.47 0.35 -0.16 286868
97 0.303 0.196 0.198 0.303 -0.81 0.18 1.19 -0.38 285474 0.302 0.197 0.198 0.303 -0.99 0.78 0.48 -0.09 286808
98 0.305 0.197 0.195 0.303 1.01 1.46 -3.08 0.39 285142 0.304 0.197 0.196 0.304 0.62 0.16 -2.15 1.11 286537
99 0.305 0.195 0.197 0.303 1.18 -1.41 -0.24 0.25 284843 0.304 0.195 0.198 0.303 1.08 -1.66 0.39 0.01 286266
100 0.303 0.197 0.198 0.302 -0.94 1.59 1.32 -1.58 284762 0.303 0.197 0.198 0.302 -0.18 1.06 0.58 -1.24 285952
101 0.304 0.197 0.196 0.303 -0.17 0.99 -1.40 0.53 284878 0.304 0.197 0.197 0.303 0.70 0.40 -0.94 -0.23 285628
102 0.305 0.196 0.197 0.303 0.88 -0.52 -0.13 -0.32 284524 0.304 0.196 0.198 0.302 0.64 -0.79 0.66 -0.53 285274
103 0.303 0.196 0.197 0.304 -1.34 -0.10 0.67 0.86 284346 0.302 0.197 0.198 0.302 -0.66 0.80 0.70 -0.64 284978
104 0.305 0.196 0.195 0.304 0.92 -0.55 -2.02 1.30 284151 0.303 0.197 0.197 0.304 -0.43 0.09 -1.10 1.32 284708
105 0.304 0.195 0.197 0.304 0.49 -1.33 -0.49 1.09 284006 0.304 0.195 0.198 0.304 0.76 -2.88 0.49 1.31 284917
106 0.303 0.197 0.197 0.303 -0.80 0.48 -0.08 0.45 283675 0.303 0.197 0.198 0.303 -0.54 0.32 0.41 -0.09 284713
107 0.304 0.196 0.195 0.304 0.24 0.15 -2.27 1.59 283333 0.304 0.196 0.196 0.304 0.96 -0.31 -1.96 1.01 284375
108 0.305 0.196 0.198 0.302 1.12 -0.95 0.81 -1.00 283331 0.304 0.195 0.199 0.302 1.21 -2.39 1.42 -0.38 284129
109 0.302 0.197 0.198 0.303 -1.63 1.16 1.02 -0.25 283228 0.302 0.197 0.198 0.302 -1.12 0.87 1.23 -0.70 284116
110 0.305 0.196 0.196 0.304 1.01 -1.01 -0.92 0.66 283088 0.304 0.197 0.195 0.304 1.06 -0.08 -2.96 1.57 283896
111 0.304 0.196 0.197 0.303 0.74 -0.84 0.10 -0.10 282798 0.304 0.196 0.198 0.302 1.13 -1.43 1.10 -0.85 284151
112 0.303 0.196 0.198 0.303 -1.41 0.04 1.27 0.28 282514 0.303 0.197 0.199 0.302 -0.55 0.76 1.46 -1.37 284213
113 0.304 0.197 0.195 0.304 0.01 1.32 -2.58 1.07 282783 0.302 0.197 0.197 0.304 -0.75 1.00 -1.35 1.06 284000
114 0.305 0.195 0.197 0.303 1.38 -1.51 0.38 -0.41 282711 0.304 0.196 0.198 0.302 1.26 -0.95 0.83 -1.16 283749
115 0.302 0.197 0.198 0.303 -1.82 0.70 1.48 -0.06 282999 0.301 0.197 0.199 0.303 -2.01 0.90 1.55 -0.11 283806
116 0.305 0.197 0.194 0.304 1.07 0.78 -3.57 1.35 283311 0.304 0.197 0.195 0.304 1.56 0.31 -3.69 1.37 283963
117 0.305 0.195 0.197 0.303 1.61 -1.88 0.58 -0.48 283509 0.304 0.194 0.198 0.303 1.60 -3.29 1.11 0.28 284407
118 0.302 0.197 0.197 0.304 -2.55 1.64 0.26 0.90 283567 0.302 0.198 0.196 0.304 -0.93 1.19 -1.39 1.11 284812
119 0.305 0.197 0.194 0.304 0.93 0.93 -3.75 1.52 283576 0.304 0.196 0.196 0.304 1.23 -0.83 -2.26 1.45 284645
120 0.305 0.195 0.197 0.303 1.05 -1.43 0.35 -0.11 283723 0.303 0.196 0.198 0.302 0.26 -0.67 1.12 -0.66 284493
121 0.302 0.197 0.198 0.303 -2.39 1.37 1.67 -0.24 283472 0.302 0.198 0.199 0.302 -1.44 1.59 1.50 -1.23 284479
122 0.305 0.197 0.194 0.304 1.14 1.11 -4.15 1.50 283491 0.303 0.197 0.196 0.304 0.10 0.60 -2.37 1.44 284948
123 0.304 0.194 0.198 0.304 0.00 -2.48 1.49 0.86 283916 0.304 0.195 0.198 0.302 1.50 -1.84 0.70 -0.52 284733
124 0.302 0.197 0.198 0.303 -2.15 1.52 1.46 -0.42 283980 0.303 0.197 0.199 0.302 -0.17 0.17 1.36 -1.15 284587
125 0.304 0.197 0.194 0.304 0.80 1.56 -3.85 1.19 284402 0.304 0.198 0.194 0.304 1.44 1.58 -4.27 0.89 284468
126 0.305 0.192 0.198 0.304 1.60 -5.18 1.61 1.49 284197 0.304 0.194 0.198 0.303 1.27 -3.45 1.08 0.78 284204
127 0.302 0.197 0.198 0.303 -1.74 1.22 1.12 -0.28 282723 0.301 0.197 0.199 0.303 -2.81 0.88 1.49 0.76 283472
128 0.305 0.197 0.193 0.304 0.97 1.62 -4.68 1.68 281313 0.304 0.197 0.196 0.304 1.02 -0.07 -2.59 1.30 282018
129 0.305 0.195 0.197 0.303 1.03 -1.65 0.48 -0.02 280767 0.304 0.196 0.198 0.302 0.89 -1.35 1.10 -0.67 281501
130 0.302 0.196 0.198 0.304 -1.84 0.03 1.26 0.73 280441 0.302 0.197 0.198 0.302 -0.82 0.71 0.98 -0.64 281025
131 0.304 0.197 0.195 0.304 0.07 1.14 -2.96 1.50 279616 0.304 0.197 0.196 0.303 0.62 0.85 -1.77 0.18 280676
132 0.305 0.196 0.197 0.302 1.36 -0.82 0.05 -0.70 279288 0.304 0.196 0.197 0.303 0.58 -0.36 -0.39 0.07 280390
133 0.303 0.197 0.198 0.303 -0.90 0.54 1.02 -0.44 278873 0.302 0.197 0.199 0.303 -1.09 0.01 1.48 -0.20 279883
134 0.303 0.197 0.195 0.304 -0.65 1.52 -2.30 1.33 278324 0.303 0.196 0.197 0.304 -0.53 -0.43 -0.58 1.40 279346
135 0.304 0.195 0.198 0.303 0.40 -1.13 0.98 -0.27 278047 0.304 0.196 0.198 0.302 0.75 -0.97 0.53 -0.37 278969
136 0.302 0.197 0.198 0.303 -1.79 1.42 1.28 -0.54 277507 0.302 0.197 0.198 0.302 -0.79 0.21 1.22 -0.45 278526
137 0.304 0.197 0.195 0.304 -0.11 1.39 -2.25 0.86 277180 0.303 0.197 0.196 0.304 0.11 0.31 -1.90 1.27 278212
138 0.305 0.195 0.197 0.303 1.28 -1.58 -0.35 0.38 276489 0.304 0.196 0.198 0.302 1.34 -1.36 0.26 -0.39 277631
139 0.302 0.196 0.198 0.304 -2.11 -0.04 0.92 1.35 276101 0.301 0.197 0.198 0.304 -2.35 0.48 0.87 1.18 277099
140 0.305 0.196 0.196 0.304 1.62 -0.66 -1.97 0.66 275270 0.304 0.197 0.195 0.304 0.90 0.01 -2.91 1.61 276304
141 0.304 0.195 0.198 0.303 0.15 -1.45 1.36 -0.07 274897 0.304 0.195 0.199 0.302 1.32 -2.15 1.47 -0.74 275806
142 0.303 0.196 0.198 0.302 -0.33 0.20 0.87 -0.60 274398 0.302 0.197 0.199 0.302 -0.77 0.20 1.33 -0.56 275418
143 0.303 0.197 0.196 0.304 -0.77 0.85 -1.14 1.03 274096 0.303 0.197 0.197 0.303 -0.07 0.57 -0.39 -0.09 274895
144 0.304 0.196 0.197 0.303 0.69 -0.99 0.57 -0.32 273757 0.304 0.196 0.197 0.303 1.13 -0.37 -0.64 -0.25 274480
145 0.303 0.197 0.198 0.302 -1.02 0.98 0.87 -0.58 273492 0.302 0.198 0.198 0.303 -1.42 1.28 0.38 -0.02 274154
146 0.303 0.197 0.196 0.304 -1.25 1.15 -1.28 1.36 273235 0.303 0.197 0.197 0.303 0.30 0.55 -1.19 0.25 273902
147 0.305 0.196 0.197 0.302 1.41 -0.49 -0.44 -0.60 272944 0.304 0.195 0.198 0.303 1.56 -2.64 0.80 0.03 273506
148 0.303 0.197 0.198 0.303 -1.38 1.13 0.87 -0.36 272541 0.302 0.197 0.199 0.302 -1.13 0.54 1.31 -0.47 272696
149 0.304 0.197 0.195 0.304 0.50 0.46 -2.51 1.28 272265 0.303 0.197 0.196 0.304 -0.08 0.68 -2.46 1.62 272166
150 0.304 0.196 0.198 0.302 0.37 -0.88 1.13 -0.58 271966 0.304 0.195 0.198 0.302 0.66 -1.62 1.27 -0.36 271763
151 0.303 0.196 0.198 0.303 -1.10 0.12 0.69 0.40 271466 0.303 0.196 0.198 0.302 -0.08 -0.41 1.10 -0.52 271400
152 0.304 0.197 0.195 0.304 0.48 0.99 -2.68 0.99 271179 0.303 0.197 0.196 0.303 0.41 0.43 -1.56 0.58 271090
153 0.304 0.195 0.197 0.304 0.20 -1.05 0.07 0.65 270816 0.303 0.195 0.198 0.304 0.12 -1.63 0.10 1.20 270839
154 0.303 0.197 0.198 0.303 -1.28 0.45 1.39 -0.32 270527 0.303 0.197 0.198 0.302 0.08 0.15 0.20 -0.39 270490
155 0.304 0.196 0.195 0.304 0.76 0.22 -2.89 1.55 270255 0.303 0.197 0.196 0.304 -0.17 0.33 -1.41 1.11 270170
156 0.304 0.196 0.198 0.303 -0.08 -0.32 0.91 -0.43 269996 0.303 0.197 0.198 0.302 0.19 -0.09 0.76 -0.77 269910
157 0.302 0.197 0.198 0.302 -1.46 1.52 1.32 -0.99 269590 0.303 0.198 0.197 0.302 -0.39 1.37 -0.03 -0.77 269328
158 0.304 0.197 0.195 0.304 0.33 0.44 -2.12 1.12 268870 0.303 0.197 0.196 0.304 0.25 0.04 -1.46 0.98 268820
159 0.304 0.196 0.197 0.303 0.03 -0.22 0.51 -0.28 268603 0.303 0.197 0.198 0.302 -0.34 0.38 1.08 -0.93 268470
160 0.303 0.197 0.197 0.302 -0.48 1.47 -0.17 -0.65 268338 0.302 0.197 0.198 0.303 -0.99 0.81 0.11 0.19 268114
161 0.304 0.196 0.195 0.304 0.66 -0.59 -2.05 1.62 267906 0.303 0.197 0.196 0.304 0.17 -0.26 -1.51 1.36 267612
162 0.304 0.197 0.197 0.302 0.56 0.38 0.26 -1.11 267637 0.304 0.197 0.198 0.301 1.24 0.08 0.18 -1.47 267267
163 0.302 0.197 0.198 0.303 -1.43 0.70 1.02 -0.06 267349 0.303 0.197 0.197 0.302 0.30 0.68 -0.16 -0.76 266995
164 0.304 0.196 0.196 0.304 0.07 0.19 -1.30 0.89 267048 0.303 0.197 0.196 0.304 0.05 0.48 -1.69 0.99 266717
165 0.305 0.196 0.198 0.301 1.18 -0.76 1.37 -1.71 266538 0.304 0.197 0.198 0.301 0.70 0.26 1.27 -2.02 266222
166 0.304 0.197 0.197 0.302 0.04 0.36 0.23 -0.55 266212 0.304 0.197 0.198 0.302 0.74 -0.21 0.11 -0.65 265971
167 0.304 0.196 0.196 0.304 -0.02 0.05 -1.74 1.48 265948 0.303 0.197 0.197 0.303 -0.05 -0.03 -0.15 0.21 265710
168 0.304 0.195 0.197 0.304 0.60 -1.07 -0.43 0.70 265692 0.304 0.197 0.197 0.302 0.77 0.01 -0.27 -0.54 265474
169 0.302 0.198 0.198 0.303 -1.85 1.63 1.02 -0.43 265271 0.302 0.197 0.198 0.303 -1.48 0.25 0.73 0.63 265193
170 0.304 0.197 0.195 0.304 0.12 0.95 -1.98 0.77 264879 0.302 0.197 0.197 0.304 -0.62 -0.21 -0.53 1.26 264797
171 0.303 0.196 0.197 0.303 -0.73 0.08 0.14 0.54 264578 0.304 0.196 0.197 0.303 0.60 -0.29 -0.58 0.15 264379
172 0.302 0.197 0.197 0.303 -1.61 1.54 0.25 0.06 264230 0.302 0.197 0.199 0.303 -1.23 0.17 1.59 -0.30 264026
173 0.304 0.196 0.196 0.303 0.47 0.21 -1.26 0.44 263891 0.303 0.197 0.196 0.303 0.07 0.60 -1.58 0.78 263687
174 0.304 0.196 0.197 0.303 0.11 -0.21 0.01 0.06 263575 0.303 0.197 0.198 0.302 0.23 0.17 0.09 -0.46 263445
175 0.303 0.197 0.198 0.301 -0.35 0.96 1.34 -1.65 263255 0.302 0.197 0.198 0.303 -0.81 0.90 0.10 -0.06 262929
176 0.304 0.196 0.196 0.304 0.30 -0.17 -1.01 0.71 262967 0.303 0.197 0.196 0.304 -0.39 0.57 -2.03 1.66 262511
177 0.303 0.196 0.197 0.303 -0.69 -0.45 0.64 0.53 262632 0.303 0.197 0.198 0.303 -0.37 0.28 0.35 -0.17 262212
178 0.302 0.197 0.198 0.303 -2.29 1.18 1.04 0.37 261984 0.302 0.197 0.198 0.303 -0.84 0.29 0.47 0.17 261895
179 0.305 0.196 0.196 0.303 1.23 -0.39 -1.28 0.20 261695 0.304 0.196 0.197 0.303 0.59 -0.55 -0.78 0.56 261608
180 0.304 0.196 0.197 0.303 0.55 -0.90 0.20 0.05 261408 0.303 0.197 0.197 0.303 -0.18 0.48 -0.33 0.04 261358
181 0.303 0.196 0.198 0.303 -0.93 -0.61 1.47 0.18 261138 0.303 0.196 0.198 0.303 -0.12 -0.69 0.95 -0.11 260970
182 0.304 0.196 0.196 0.303 0.21 0.10 -0.93 0.51 260861 0.303 0.197 0.197 0.303 0.24 0.00 -0.62 0.29 260696
183 0.303 0.196 0.198 0.303 -0.37 -0.88 0.69 0.53 260596 0.304 0.196 0.197 0.303 0.72 -1.03 -0.22 0.36 260434
184 0.303 0.197 0.197 0.303 -0.33 0.42 -0.02 -0.01 260189 0.303 0.196 0.198 0.303 -0.11 -0.70 0.53 0.26 260058
185 0.303 0.197 0.196 0.304 -1.10 1.06 -1.40 1.39 259649 0.303 0.197 0.196 0.304 -0.09 0.31 -1.41 1.05 259757
186 0.304 0.196 0.198 0.302 0.12 -0.20 1.00 -0.82 259157 0.303 0.197 0.198 0.302 0.55 0.54 0.34 -1.31 259453
187 0.303 0.197 0.197 0.302 -0.32 1.23 0.02 -0.77 258875 0.303 0.196 0.197 0.303 0.45 -0.53 -0.07 0.07 259152
188 0.304 0.197 0.195 0.304 0.37 1.04 -2.22 0.65 258471 0.303 0.197 0.197 0.303 0.12 0.26 -0.66 0.22 258878
189 0.304 0.196 0.198 0.302 0.75 -0.66 0.97 -1.01 258202 0.303 0.197 0.198 0.302 0.13 0.37 0.88 -1.21 258582
190 0.303 0.197 0.198 0.302 -0.70 0.81 1.12 -0.97 257769 0.302 0.197 0.199 0.302 -0.61 -0.01 1.55 -0.73 258065
191 0.304 0.196 0.197 0.303 0.04 -0.01 -0.21 0.15 257526 0.303 0.196 0.197 0.304 0.36 -0.88 -0.97 1.23 257801
192 0.304 0.196 0.196 0.304 0.51 -0.56 -0.99 0.83 257222 0.302 0.196 0.198 0.304 -0.61 -0.33 0.12 0.79 257393
193 0.303 0.197 0.197 0.302 -0.31 0.93 0.11 -0.59 256896 0.303 0.197 0.198 0.303 -0.54 0.80 0.17 -0.30 257135
194 0.303 0.197 0.196 0.304 -0.94 1.18 -1.12 0.90 256575 0.303 0.197 0.197 0.302 0.43 0.76 -0.41 -0.74 256833
195 0.303 0.196 0.197 0.303 -0.38 -0.13 0.33 0.20 256315 0.303 0.196 0.198 0.303 0.21 -0.34 0.14 -0.03 256583
196 0.303 0.196 0.198 0.303 -1.25 0.27 0.75 0.37 256040 0.302 0.197 0.199 0.302 -0.63 0.21 1.38 -0.75 256316
197 0.304 0.197 0.197 0.303 0.42 0.33 -0.38 -0.38 255761 0.303 0.197 0.197 0.303 -0.48 0.63 -0.87 0.69 256052
198 0.304 0.196 0.198 0.302 0.45 -0.02 0.93 -1.24 255409 0.304 0.197 0.197 0.303 0.87 0.23 -0.88 -0.30 255611
199 0.304 0.196 0.198 0.301 0.62 0.17 1.08 -1.70 255130 0.302 0.197 0.198 0.302 -0.64 0.46 1.01 -0.63 255275
200 0.304 0.197 0.197 0.303 0.36 0.48 -0.40 -0.42 254816 0.303 0.198 0.197 0.302 -0.26 1.56 -0.84 -0.36 255012
201 0.304 0.196 0.197 0.303 0.67 -0.58 0.27 -0.40 254566 0.303 0.197 0.198 0.303 -0.49 0.26 0.63 -0.28 254748
202 0.302 0.197 0.198 0.303 -1.81 0.45 1.14 0.44 254264 0.301 0.197 0.199 0.303 -1.95 0.94 1.38 -0.06 254403
203 0.304 0.197 0.196 0.303 0.23 0.67 -1.19 0.22 253958 0.303 0.197 0.197 0.303 -0.35 0.21 -0.67 0.75 254079
204 0.303 0.196 0.198 0.303 -0.63 -0.07 1.15 -0.30 253693 0.303 0.196 0.198 0.303 0.02 -0.83 1.15 -0.29 253730
205 0.303 0.196 0.197 0.303 -0.47 0.27 0.46 -0.17 253364 0.302 0.197 0.198 0.302 -0.84 0.82 0.88 -0.63 253473
206 0.303 0.197 0.196 0.304 -0.54 0.93 -1.05 0.65 253043 0.303 0.198 0.196 0.303 0.13 1.33 -1.80 0.28 253233
207 0.304 0.196 0.198 0.302 0.33 -0.79 1.07 -0.56 252763 0.303 0.197 0.198 0.302 0.37 -0.09 0.25 -0.51 252945
208 0.303 0.197 0.198 0.302 -1.24 1.02 1.17 -0.65 252442 0.302 0.197 0.198 0.303 -1.52 -0.22 1.21 0.66 252408
209 0.303 0.197 0.197 0.303 -0.55 1.10 -0.07 -0.34 252130 0.302 0.198 0.197 0.303 -0.72 1.21 -0.40 0.02 252118
210 0.303 0.197 0.197 0.302 -0.47 0.92 0.35 -0.62 251857 0.303 0.197 0.198 0.303 -0.42 0.05 0.57 -0.11 251750
211 0.303 0.196 0.198 0.303 -0.66 0.10 0.89 -0.20 251573 0.303 0.197 0.198 0.303 -0.43 -0.07 0.86 -0.25 251471
212 0.303 0.197 0.197 0.303 -0.34 0.72 0.13 -0.39 251261 0.302 0.197 0.197 0.303 -0.66 0.31 -0.04 0.43 251173
213 0.304 0.196 0.197 0.303 0.28 -0.16 0.33 -0.43 250929 0.303 0.198 0.198 0.302 -0.46 1.04 0.27 -0.67 250862
214 0.304 0.196 0.198 0.302 -0.18 0.20 1.40 -1.21 250448 0.303 0.197 0.198 0.302 -0.24 0.80 0.85 -1.19 250348
215 0.304 0.196 0.197 0.303 -0.15 0.18 -0.40 0.35 250131 0.303 0.197 0.197 0.303 -0.48 0.63 -0.63 0.48 250055
216 0.305 0.196 0.197 0.303 0.81 -0.37 -0.06 -0.44 249815 0.303 0.197 0.198 0.302 0.13 0.00 0.38 -0.46 249754
217 0.303 0.197 0.197 0.303 -0.84 1.10 0.09 -0.18 249177 0.303 0.197 0.198 0.302 -0.43 0.87 0.08 -0.39 249483
218 0.304 0.196 0.196 0.303 0.72 -0.01 -1.41 0.51 248826 0.302 0.197 0.197 0.303 -1.12 0.62 -0.22 0.77 249184
219 0.305 0.196 0.198 0.302 0.91 -0.89 0.66 -0.71 248543 0.303 0.196 0.197 0.303 0.20 -0.32 -0.16 0.21 248899
220 0.303 0.197 0.197 0.303 -0.69 0.70 0.34 -0.20 248250 0.302 0.198 0.198 0.302 -1.24 1.28 0.79 -0.56 248620
221 0.303 0.196 0.197 0.304 -0.38 -0.80 0.12 0.96 248007 0.303 0.197 0.197 0.303 -0.19 0.60 -0.62 0.21 248322
222 0.304 0.196 0.198 0.303 -0.12 -0.36 0.96 -0.39 247646 0.303 0.197 0.199 0.302 -0.24 0.61 1.24 -1.35 247944
223 0.303 0.197 0.198 0.302 -1.01 1.12 0.71 -0.57 247360 0.302 0.197 0.199 0.302 -1.07 0.97 1.38 -0.96 247643
224 0.303 0.197 0.196 0.304 -0.32 0.83 -1.44 0.85 247060 0.302 0.197 0.197 0.304 -0.74 0.32 -0.46 0.86 247351
225 0.304 0.197 0.197 0.302 0.29 0.32 0.43 -0.93 246676 0.302 0.197 0.198 0.303 -1.44 0.41 0.70 0.48 247061
226 0.303 0.197 0.198 0.302 -0.75 1.09 0.79 -0.88 246389 0.302 0.197 0.199 0.302 -1.12 0.47 1.61 -0.69 246783
227 0.304 0.196 0.196 0.303 -0.05 0.20 -0.77 0.54 246119 0.303 0.197 0.197 0.302 -0.16 0.74 -0.19 -0.32 246515
228 0.304 0.196 0.198 0.303 -0.05 -0.71 1.22 -0.39 245816 0.302 0.197 0.198 0.303 -1.34 0.76 1.11 -0.28 245989
229 0.303 0.197 0.197 0.303 -1.24 1.52 -0.29 0.18 245479 0.302 0.197 0.198 0.303 -0.92 0.80 0.26 0.00 245624
230 0.304 0.197 0.196 0.304 -0.15 0.32 -1.14 0.86 245226 0.303 0.197 0.197 0.303 -0.45 0.09 -0.42 0.73 245370
231 0.303 0.197 0.197 0.303 -0.85 0.87 0.28 -0.14 244929 0.302 0.198 0.198 0.302 -1.23 1.17 1.21 -0.83 245060
232 0.303 0.197 0.197 0.303 -0.80 0.51 -0.22 0.55 244686 0.302 0.198 0.198 0.303 -1.14 1.25 0.33 -0.23 244776
233 0.303 0.197 0.196 0.304 -0.83 0.80 -0.84 0.86 244385 0.303 0.197 0.198 0.303 -0.47 0.63 0.02 -0.09 244484
234 0.304 0.197 0.197 0.303 -0.18 0.48 -0.45 0.16 244130 0.301 0.198 0.198 0.303 -2.00 1.27 0.66 0.33 244231
235 0.303 0.197 0.198 0.302 -1.28 1.38 1.25 -1.00 243786 0.300 0.198 0.198 0.304 -2.72 1.24 0.59 1.14 243909
236 0.304 0.197 0.195 0.304 -0.25 1.42 -1.87 0.64 243501 0.302 0.197 0.198 0.303 -1.38 0.67 0.87 0.05 243660
237 0.303 0.197 0.197 0.303 -1.06 1.19 0.61 -0.50 243220 0.302 0.197 0.198 0.302 -0.71 0.48 0.99 -0.56 243357
238 0.304 0.196 0.198 0.303 -0.25 -0.34 0.68 -0.04 242911 0.302 0.197 0.198 0.302 -1.13 0.95 0.75 -0.34 242786
239 0.304 0.197 0.196 0.303 0.53 0.70 -1.63 0.27 242637 0.302 0.198 0.197 0.302 -0.67 1.38 -0.12 -0.41 242495
240 0.303 0.197 0.197 0.303 -0.32 1.07 -0.12 -0.50 242295 0.302 0.197 0.198 0.302 -0.87 0.65 1.10 -0.64 242175
241 0.302 0.197 0.198 0.303 -1.36 1.30 0.72 -0.40 241998 0.302 0.197 0.198 0.303 -1.49 0.46 1.15 0.09 241888
242 0.304 0.196 0.196 0.303 0.46 0.01 -0.92 0.33 241728 0.302 0.197 0.198 0.302 -0.57 0.77 0.41 -0.46 241532
243 0.304 0.196 0.198 0.302 0.02 0.14 0.95 -0.96 241467 0.303 0.197 0.198 0.303 -0.46 0.28 0.16 0.08 241270
244 0.303 0.197 0.198 0.302 -0.72 1.02 0.83 -0.88 241200 0.302 0.197 0.198 0.303 -1.36 0.09 1.16 0.28 240938
245 0.303 0.197 0.196 0.304 -0.63 0.29 -1.03 1.26 240848 0.302 0.198 0.197 0.303 -1.45 1.58 -0.09 0.17 240522
246 0.303 0.197 0.197 0.302 -0.50 1.28 0.09 -0.69 240464 0.302 0.198 0.198 0.302 -1.24 1.48 0.55 -0.52 239928
247 0.303 0.196 0.198 0.303 -0.68 0.15 0.71 -0.06 240073 0.302 0.197 0.198 0.303 -0.84 0.79 0.02 0.14 239458
248 0.303 0.196 0.197 0.304 -0.88 -0.33 -0.14 1.29 239088 0.303 0.197 0.197 0.303 -0.29 0.34 -0.48 0.41 238424
249 0.305 0.193 0.198 0.304 1.06 -4.57 1.59 1.50 235006 0.304 0.194 0.199 0.303 0.91 -3.27 1.57 0.56 234203
250 0.305 0.191 0.198 0.305 1.64 -5.41 1.61 1.64 186469 0.305 0.192 0.199 0.304 1.61 -5.34 1.64 1.59 190051
251 0 0 0 0 . . . . 0 0 0 0 0 . . . . 0As expected, read cropped using HCAP induce moderate z-scores (see final results above), assessing homogeneous composition among all positions (see fastq_info reports below).
 |
 |
Note however that the above fastq_info reports do not represent the initial sequenced positions (as some reads were cropped at the 5' end). To not modify the initial lengths of the input reads, HCAP can be run with the option -n to replace each cropped base by the character state N:
HCAP -1 ERR6798852_1.fastq.gz -2 ERR6798852_2.fastq.gz -o ERR6798852n -nThe fastq_info reports derived from these (uncropped) reads are show below. Note that the higher number of errors per read is trivially caused by the Phred score character # associated with each character state N.
 |
 |
In all cases, one can observe that the last position was systematically discarded because of the absence of the base C at position 251.