Take in a fastq file, output a tab separated table with basic statistics of the quality score at each base position.
Output with ex.fastq:
BP Position Mean Score Variance StDev Median
0 33.8 0.16 0.4 34.0
1 33.4 1.44 1.2 34.0
2 34.0 0.0 0.0 34.0
3 37.0 0.0 0.0 37.0
4 36.6 0.64 0.8 37.0
5 36.6 0.64 0.8 37.0
6 36.6 0.64 0.8 37.0
7 37.0 0.0 0.0 37.0
8 38.8 0.16 0.4 39.0
9 38.8 0.16 0.4 39.0
10 39.0 0.0 0.0 39.0
11 39.0 0.0 0.0 39.0
12 38.2 2.56 1.6 39.0
13 40.2 1.36 1.16619037897 41.0
14 40.6 0.24 0.489897948557 41.0
15 40.8 0.16 0.4 41.0
16 41.0 0.0 0.0 41.0
17 41.0 0.0 0.0 41.0
18 40.6 0.24 0.489897948557 41.0
19 41.0 0.0 0.0 41.0
20 40.2 1.36 1.16619037897 41.0
21 40.4 1.44 1.2 41.0
22 40.4 0.64 0.8 41.0
23 39.8 3.76 1.93907194297 41.0
24 40.0 1.2 1.09544511501 40.0
25 39.4 3.04 1.74355957742 40.0
26 40.6 0.64 0.8 41.0
27 40.0 4.0 2.0 41.0
28 40.2 0.56 0.748331477355 40.0
29 38.8 2.56 1.6 40.0
30 38.4 15.04 3.87814388593 41.0
31 39.2 4.96 2.22710574513 40.0
32 39.4 2.64 1.62480768093 40.0
33 40.0 1.2 1.09544511501 40.0
34 39.2 9.76 3.12409987036 41.0
35 40.2 1.36 1.16619037897 41.0
36 40.8 0.16 0.4 41.0
37 40.8 0.16 0.4 41.0
38 41.0 0.0 0.0 41.0
39 40.4 1.44 1.2 41.0
40 40.6 0.24 0.489897948557 41.0
41 40.0 1.2 1.09544511501 40.0
42 40.8 0.16 0.4 41.0
43 40.0 4.0 2.0 41.0
44 40.2 0.56 0.748331477355 40.0
45 40.8 0.16 0.4 41.0
46 40.2 0.16 0.4 40.0
47 39.4 1.44 1.2 40.0
48 40.6 0.24 0.489897948557 41.0
49 40.2 1.36 1.16619037897 41.0
50 39.8 3.76 1.93907194297 41.0
51 39.2 7.36 2.71293199325 41.0
52 38.4 22.24 4.71593044902 41.0
53 39.6 3.44 1.8547236991 40.0
54 39.8 3.76 1.93907194297 41.0
55 39.0 4.4 2.09761769634 40.0
56 40.0 1.6 1.26491106407 41.0
57 38.8 2.96 1.72046505341 39.0
58 38.6 1.04 1.01980390272 39.0
59 38.2 4.16 2.03960780544 39.0
60 36.8 29.76 5.45527267879 39.0
61 37.4 8.24 2.87054001888 39.0
62 38.0 4.0 2.0 38.0
63 36.2 28.56 5.34415568635 37.0
64 36.6 16.24 4.02988833592 37.0
65 37.6 7.84 2.8 37.0
66 36.6 11.84 3.44093010682 37.0
67 37.0 5.2 2.2803508502 37.0
68 35.6 9.84 3.13687742827 36.0
69 36.0 8.8 2.96647939484 36.0
70 36.2 4.56 2.13541565041 36.0
71 35.4 6.64 2.57681974535 36.0
72 36.4 5.44 2.33238075794 36.0
73 35.8 5.36 2.31516738056 36.0
74 35.8 1.36 1.16619037897 36.0
75 34.0 2.4 1.54919333848 35.0
76 34.4 7.84 2.8 35.0
77 35.0 4.8 2.19089023002 35.0
78 35.4 2.24 1.49666295471 35.0
79 35.6 1.44 1.2 35.0
80 35.6 1.04 1.01980390272 36.0
81 35.0 0.4 0.632455532034 35.0
82 32.8 7.36 2.71293199325 35.0
83 33.2 2.56 1.6 33.0
84 34.6 1.84 1.35646599663 35.0
85 35.0 0.4 0.632455532034 35.0
86 34.2 2.16 1.46969384567 35.0
87 34.8 2.56 1.6 35.0
88 34.2 2.96 1.72046505341 35.0
89 35.0 0.4 0.632455532034 35.0
90 35.0 0.0 0.0 35.0
91 34.8 0.16 0.4 35.0
92 34.8 0.16 0.4 35.0
93 34.4 0.64 0.8 35.0
94 32.8 15.36 3.91918358845 35.0
95 34.2 2.16 1.46969384567 35.0
96 34.8 0.16 0.4 35.0
97 35.2 0.16 0.4 35.0
98 33.2 7.36 2.71293199325 35.0
99 32.0 18.4 4.28952211791 35.0
100 30.4 16.64 4.07921561087 30.0