Starting /dee2/code/volunteer_pipeline.sh ERR4018149
    current disk space = 375937592
    free memory = 194996296 
ERR4018149 SRAfilesize
55a1bce67d56d71ad11aa7cd98a6cc65  ERR4018149.sra
ERR4018149.sra file validated
ERR4018149 is single end
ERR4018149 is conventional basespace
ERR4018149 read1 length is 51 nt
##FastQC	0.11.5
>>Basic Statistics	pass
#Measure	Value
Filename	ERR4018149_1.fastq
File type	Conventional base calls
Encoding	Sanger / Illumina 1.9
Total Sequences	4000
Sequences flagged as poor quality	0
Sequence length	51
%GC	46
>>END_MODULE
>>Per base sequence quality	pass
#Base	Mean	Median	Lower Quartile	Upper Quartile	10th Percentile	90th Percentile
1	31.4315	35.0	35.0	35.0	27.0	35.0
2	33.78875	35.0	35.0	35.0	27.0	35.0
3	34.23725	35.0	35.0	35.0	33.0	35.0
4	34.43775	35.0	35.0	35.0	33.0	35.0
5	33.72075	35.0	35.0	35.0	31.0	35.0
6	38.38475	39.0	38.0	40.0	36.0	40.0
7	38.99375	40.0	39.0	40.0	38.0	40.0
8	39.102	40.0	39.0	40.0	38.0	40.0
9	39.00975	40.0	39.0	40.0	38.0	40.0
10	39.1855	40.0	39.0	40.0	38.0	40.0
11	39.08125	40.0	39.0	40.0	38.0	40.0
12	39.16675	40.0	39.0	40.0	38.0	40.0
13	39.108	40.0	39.0	40.0	38.0	40.0
14	39.09825	40.0	39.0	40.0	38.0	40.0
15	39.14325	40.0	39.0	40.0	38.0	40.0
16	39.16525	40.0	39.0	40.0	38.0	40.0
17	39.2205	40.0	39.0	40.0	38.0	40.0
18	39.09575	40.0	39.0	40.0	38.0	40.0
19	39.164	40.0	39.0	40.0	38.0	40.0
20	39.17525	40.0	39.0	40.0	38.0	40.0
21	39.2875	40.0	39.0	40.0	38.0	40.0
22	39.243	40.0	39.0	40.0	38.0	40.0
23	39.128	40.0	39.0	40.0	38.0	40.0
24	39.19675	40.0	39.0	40.0	38.0	40.0
25	39.153	40.0	39.0	40.0	38.0	40.0
26	39.16275	40.0	39.0	40.0	38.0	40.0
27	39.234	40.0	39.0	40.0	38.0	40.0
28	39.21675	40.0	39.0	40.0	38.0	40.0
29	39.19925	40.0	39.0	40.0	38.0	40.0
30	39.11	40.0	39.0	40.0	38.0	40.0
31	39.07975	40.0	39.0	40.0	38.0	40.0
32	39.07825	40.0	39.0	40.0	38.0	40.0
33	39.05575	40.0	39.0	40.0	38.0	40.0
34	39.114	40.0	39.0	40.0	38.0	40.0
35	39.045	40.0	39.0	40.0	38.0	40.0
36	38.91475	40.0	39.0	40.0	38.0	40.0
37	39.002	40.0	39.0	40.0	38.0	40.0
38	39.02425	40.0	39.0	40.0	38.0	40.0
39	38.98225	40.0	39.0	40.0	38.0	40.0
40	38.97975	40.0	39.0	40.0	38.0	40.0
41	39.0285	40.0	39.0	40.0	38.0	40.0
42	39.03575	40.0	39.0	40.0	38.0	40.0
43	39.064	40.0	39.0	40.0	38.0	40.0
44	39.0195	40.0	39.0	40.0	38.0	40.0
45	38.93675	40.0	39.0	40.0	37.0	40.0
46	39.01175	40.0	39.0	40.0	38.0	40.0
47	39.05475	40.0	39.0	40.0	38.0	40.0
48	39.0625	40.0	39.0	40.0	38.0	40.0
49	38.6885	40.0	39.0	40.0	37.0	40.0
50	38.87575	40.0	39.0	40.0	37.0	40.0
51	38.4155	39.0	39.0	40.0	36.0	40.0
>>END_MODULE
>>Per sequence quality scores	pass
#Quality	Count
16	1.0
17	0.0
18	0.0
19	0.0
20	0.0
21	0.0
22	0.0
23	1.0
24	2.0
25	1.0
26	2.0
27	5.0
28	3.0
29	5.0
30	13.0
31	20.0
32	20.0
33	34.0
34	53.0
35	79.0
36	133.0
37	291.0
38	1194.0
39	2143.0
>>END_MODULE
>>Per base sequence content	fail
#Base	G	A	T	C
1	23.27372764786795	20.577716643741404	35.378266850068776	20.77028885832187
2	44.4	14.475	31.075000000000003	10.05
3	22.325	36.075	10.925	30.675
4	15.45	38.425	21.175	24.95
5	20.549999999999997	37.925	15.25	26.275
6	26.05	19.5	44.025	10.424999999999999
7	57.275	6.5	27.650000000000002	8.575000000000001
8	27.525	49.875	10.85	11.75
9	10.95	40.125	40.300000000000004	8.625
10	51.225	13.55	25.324999999999996	9.9
11	39.525	44.55	9.45	6.4750000000000005
12	29.175	28.225	16.875	25.724999999999998
13	16.675	31.075000000000003	20.45	31.8
14	10.45	38.125	34.65	16.775000000000002
15	16.725	20.65	33.525	29.099999999999998
16	13.275	34.300000000000004	33.575	18.85
17	23.200000000000003	23.150000000000002	17.175	36.475
18	19.575	32.550000000000004	29.175	18.7
19	14.424999999999999	32.824999999999996	33.800000000000004	18.95
20	15.825	23.45	32.225	28.499999999999996
21	23.125	14.649999999999999	52.7	9.525
22	32.6	22.625	27.900000000000002	16.875
23	25.275	31.85	11.375	31.5
24	31.025000000000002	30.25	18.35	20.375
25	31.374999999999996	23.925	13.475000000000001	31.225
26	23.674999999999997	21.7	26.974999999999998	27.650000000000002
27	29.525000000000002	19.525000000000002	13.0	37.95
28	18.625	20.349999999999998	27.900000000000002	33.125
29	27.85	21.175	27.825	23.150000000000002
30	33.775	14.549999999999999	25.324999999999996	26.35
31	23.849999999999998	25.8	31.65	18.7
32	27.725	21.675	21.075	29.525000000000002
33	17.375	32.324999999999996	17.525	32.775
34	13.55	17.4	33.175	35.875
35	22.175	20.625	23.45	33.75
36	24.099999999999998	37.525	19.85	18.525
37	11.600000000000001	42.8	27.3	18.3
38	30.15	31.374999999999996	12.5	25.974999999999998
39	23.75	23.65	36.375	16.225
40	30.175	23.849999999999998	20.674999999999997	25.3
41	32.9	27.725	22.675	16.7
42	20.25	35.825	32.074999999999996	11.85
43	18.3	29.725	23.150000000000002	28.825
44	24.075	33.725	26.35	15.85
45	21.25	39.5	19.35	19.900000000000002
46	24.45	29.549999999999997	22.925	23.075000000000003
47	14.475	27.275	38.550000000000004	19.7
48	25.6	38.425	21.05	14.924999999999999
49	11.325000000000001	36.0	37.3	15.375
50	25.174999999999997	19.875	19.975	34.975
51	25.0	21.925	24.5	28.575
>>END_MODULE
>>Per sequence GC content	fail
#GC Content	Count
0	0.0
1	0.0
2	0.0
3	0.0
4	0.0
5	0.0
6	0.0
7	0.0
8	0.0
9	0.0
10	0.0
11	0.0
12	0.0
13	0.0
14	0.0
15	0.0
16	0.0
17	0.0
18	0.0
19	0.0
20	0.0
21	0.0
22	0.0
23	0.5
24	1.0
25	1.5
26	2.5
27	3.0
28	14.0
29	25.0
30	36.5
31	48.0
32	29.0
33	10.0
34	46.0
35	82.0
36	279.0
37	476.0
38	429.0
39	382.0
40	383.5
41	385.0
42	247.0
43	109.0
44	124.0
45	139.0
46	275.5
47	412.0
48	528.5
49	645.0
50	555.0
51	465.0
52	348.5
53	232.0
54	210.5
55	189.0
56	170.0
57	151.0
58	134.5
59	118.0
60	85.5
61	53.0
62	34.5
63	16.0
64	11.5
65	7.0
66	12.5
67	18.0
68	18.0
69	18.0
70	14.5
71	11.0
72	6.5
73	2.0
74	1.0
75	0.0
76	0.0
77	0.5
78	1.0
79	0.5
80	0.0
81	0.0
82	0.0
83	0.0
84	0.0
85	0.0
86	0.0
87	0.0
88	0.0
89	0.0
90	0.0
91	0.0
92	0.0
93	0.0
94	0.0
95	0.0
96	0.0
97	0.0
98	0.0
99	0.0
100	0.0
>>END_MODULE
>>Per base N content	warn
#Base	N-Count
1	9.125
2	0.0
3	0.0
4	0.0
5	0.0
6	0.0
7	0.0
8	0.0
9	0.0
10	0.0
11	0.0
12	0.0
13	0.0
14	0.0
15	0.0
16	0.0
17	0.0
18	0.0
19	0.0
20	0.0
21	0.0
22	0.0
23	0.0
24	0.0
25	0.0
26	0.0
27	0.0
28	0.0
29	0.0
30	0.0
31	0.0
32	0.0
33	0.0
34	0.0
35	0.0
36	0.0
37	0.0
38	0.0
39	0.0
40	0.0
41	0.0
42	0.0
43	0.0
44	0.0
45	0.0
46	0.0
47	0.0
48	0.0
49	0.0
50	0.0
51	0.0
>>END_MODULE
>>Sequence Length Distribution	pass
#Length	Count
51	4000.0
>>END_MODULE
>>Sequence Duplication Levels	fail
#Total Deduplicated Percentage	34.25
#Duplication Level	Percentage of deduplicated	Percentage of total
1	80.43795620437956	27.55
2	7.153284671532846	4.9
3	3.722627737226277	3.8249999999999997
4	1.6058394160583942	2.1999999999999997
5	0.8759124087591241	1.5
6	0.8759124087591241	1.7999999999999998
7	0.8029197080291971	1.925
8	0.583941605839416	1.6
9	0.291970802919708	0.8999999999999999
>10	2.9197080291970803	19.7
>50	0.291970802919708	6.950000000000001
>100	0.43795620437956206	27.150000000000002
>500	0.0	0.0
>1k	0.0	0.0
>5k	0.0	0.0
>10k+	0.0	0.0
>>END_MODULE
>>Overrepresented sequences	fail
#Sequence	Count	Percentage	Possible Source
AGATCGGAAGAGCACACGTCTGAACTCCAGTCACCTTGTAATCTCGTATGC	321	8.025	TruSeq Adapter, Index 12 (100% over 50bp)
TGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAACT	185	4.625	No Hit
TGCACTGATGAAACTAGAAATAGGAAGTGCCGTCCGAAACGACAACTGACG	172	4.3	No Hit
GGAAGAGCACACGTCTGAACTCCAGTCACCTTGTAATCTCGTATGCCGTCT	152	3.8	TruSeq Adapter, Index 12 (100% over 51bp)
GTAAATGATGAGATTTCCATTGGTCCGTGTTTCTGAACTCCATGATTTTTG	145	3.6249999999999996	No Hit
GTGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAAC	111	2.775	No Hit
TTGCTGTGATGACTATCTTAGGACACCTTTGGATTAACCGTGAAATAAAAC	83	2.075	No Hit
CTCGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGTAGATCG	81	2.025	No Hit
TACTGTGATGATAACATAGTTCAGCAGACGTGACTCTGGTGAACAATCATG	63	1.575	No Hit
ATACATGATGATCTCACTCCAACTTGAACTCTCTCACTGATTACTTGATGA	51	1.275	No Hit
CACAATGATGACTTATGTTTTTGCCGTTTACCCAGCTGAGGGTTTCTTTGA	45	1.125	No Hit
CGCGACCTCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACCTTGTAAT	45	1.125	TruSeq Adapter, Index 12 (100% over 41bp)
TTGCAATGATGTGAATCTTTTGCTGAAGTGACCTTGAAGTGTAATTCCTGA	44	1.0999999999999999	No Hit
TGGAATGATGACATTCTCCGGAATCGCTGTACTGTCTTGATGAAAGTACTT	44	1.0999999999999999	No Hit
GGTCTGTGATGAAACCCATTATCAGTGGACATCTATGGATGAGAAATGCGG	32	0.8	No Hit
CACCCTGATGAACTGAATACCGCCCCAGTCTGATAGATGTGGAGAAAGGTA	29	0.7250000000000001	No Hit
NGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAACT	28	0.7000000000000001	No Hit
NGATCGGAAGAGCACACGTCTGAACTCCAGTCACCTTGTAATCTCGTATGC	28	0.7000000000000001	TruSeq Adapter, Index 12 (100% over 50bp)
TGCTGTGATGAGGACACCAGTAGGAAGTGCCGTCAGAGTCGATAACTGACG	27	0.675	No Hit
CTGCGATGATGTCATTTCTTAGGACACCTTTGGATTAATCATGAAAATAGC	27	0.675	No Hit
GCTAATGATGATACCGCACTAACTGAGCAGTTAGTAGTTGGTCCTTTGATT	24	0.6	No Hit
TAGCTTATCAGACTGATGTTGACAGATCGGAAGAGCACACGTCTGAACTCC	23	0.575	Illumina Multiplexing PCR Primer 2.01 (100% over 28bp)
NGCACTGATGAAACTAGAAATAGGAAGTGCCGTCCGAAACGACAACTGACG	20	0.5	No Hit
NGAAGAGCACACGTCTGAACTCCAGTCACCTTGTAATCTCGTATGCCGTCT	20	0.5	TruSeq Adapter, Index 12 (98% over 51bp)
TTGAGTGATGACTTGTATTGTCGGATACCCCTTCACCCCGTTCATGGGTGA	20	0.5	No Hit
GTTCTGTGATGAGACTCTGGAATGTGCGCTGGGCACAGCGCCCGAGTGAGA	20	0.5	No Hit
TCGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGTAGATCGG	18	0.44999999999999996	No Hit
CACCCTGATGAACTGAATACCGCCCCAGTCTGATAGCTGTGGAGAAAGGTA	18	0.44999999999999996	No Hit
GGTCAATGATGTGTTGGCATGTATTATCTGAATTTGCTGATGTGTCCTAAC	18	0.44999999999999996	No Hit
CTTCTGTGATGACTTTACCCAAATGACTTTCGTTCTTCTGAGTTTGCTGAA	16	0.4	No Hit
CTTAATGATGATTGTCTATTTGATTGCTTGATGCAATGTGAAAACCACATT	16	0.4	No Hit
CTCGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGCAGATCG	15	0.375	No Hit
ATCGATGTGGTGCTGCCGCGTTCTCTTCGGGCCGGGGCCTAAGCCGCGCCA	15	0.375	No Hit
AGTAAATGATGAGATTTCCATTGGTCCGTGTTTCTGAACTCCATGATTTTT	15	0.375	No Hit
ACAAGTGATGACTAACAAAGGGACTTAATACTGAAACATCTGGGATTGAAA	14	0.35000000000000003	No Hit
NTGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAAC	14	0.35000000000000003	No Hit
TTCAAGTAATCCAGGATAGGCTAGATCGGAAGAGCACACGTCTGAACTCCA	13	0.325	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
GCGGATGATGACACCTGGGTATGCTGCACATTCCCGACTGCGTCGTGGGGA	12	0.3	No Hit
CGCGACCTTAGATCGGAAGAGCACACGTCTGAACTCCAGTCACCTTGTAAT	12	0.3	TruSeq Adapter, Index 12 (100% over 41bp)
CGACTCTTAGCGGTGGATCACTCGGCTCGTGCGTCGATGAAGAACGCAGCT	12	0.3	No Hit
AACAGTGATGACCACAACATGCCATCTGAGTTGCTTTGCTGAAATCCAGAG	11	0.27499999999999997	No Hit
NTAAATGATGAGATTTCCATTGGTCCGTGTTTCTGAACTCCATGATTTTTG	11	0.27499999999999997	No Hit
CTCCCATGATGTATAGCACTGGGCTCTGACCACCATGGTGACGTGTTGTGC	11	0.27499999999999997	No Hit
CTGCAGTGATGACTTTCTTAGGACACCTTTGGATTTACCGTGAAAATAAGT	11	0.27499999999999997	No Hit
AACCCGTAGATCCGATCTTGTAGATCGGAAGAGCACACGTCTGAACTCCAG	10	0.25	Illumina Multiplexing PCR Primer 2.01 (100% over 30bp)
TTCTGTGATGAGACTCTGGAATGTGCGCTGGGCACAGCGCCCGAGTGAGAC	10	0.25	No Hit
CCGCGACCTCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACCTTGTAA	10	0.25	TruSeq Adapter, Index 12 (100% over 40bp)
CGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGTAGATCGGA	10	0.25	No Hit
TAGCTTATCAGACTGATGTTGACTAGATCGGAAGAGCACACGTCTGAACTC	10	0.25	Illumina Multiplexing PCR Primer 2.01 (100% over 27bp)
ACCCCTGTGATGTTCAGCACTGGGCTCTGATCACCCTGAGGACATGGAGCT	10	0.25	No Hit
TCGATGTGGTGCTGCCGCGTTCTCTTCGGGCCGGGGCCTAAGCCGCGCCAG	9	0.22499999999999998	No Hit
TGAGGTAGTAGTTTGTGCTGTTAGATCGGAAGAGCACACGTCTGAACTCCA	9	0.22499999999999998	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
NTGCTGTGATGACTATCTTAGGACACCTTTGGATTAACCGTGAAATAAAAC	9	0.22499999999999998	No Hit
TATCGATGATGTCCTAAAAGTAAATGTCTGAACATATGAATGCTATACTGA	9	0.22499999999999998	No Hit
ATATATGATGACTTTCATAGAATCTCGTTCGGCTGATGATTGCTGTTGAGA	8	0.2	No Hit
AATACATGATGATCTCACTCCAACTTGAACTCTCTCACTGATTACTTGATG	8	0.2	No Hit
CGCGACCTCAGATCAGACGTGGCGACCCGCTGAATTTAAGCATATTAGTCA	8	0.2	No Hit
NTCGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGTAGATCG	8	0.2	No Hit
TTTGAGTGATGACTTGTATTGTCGGATACCCCTTCACCCCGTTCATGGGTG	8	0.2	No Hit
ATCTATGATGATCCTATCCCGAATCTGAATTCCTGTTGAAAAAATATAATT	8	0.2	No Hit
CTACAATGATGGTAATATTTTTCGTCAACAGAGTTCACCTAGTGAGTGTTG	8	0.2	No Hit
GAGCAATGATGACCCCAACATGCCATCTGAGTTGCTTTGCTGAAATCCAGA	8	0.2	No Hit
CTAGACTGAGGCTCCTTGAGGAGATCGGAAGAGCACACGTCTGAACTCCAG	7	0.17500000000000002	Illumina Multiplexing PCR Primer 2.01 (100% over 30bp)
GTGCAATGATGTGTGAATTTCTTCACCTGAAATCACCGTGAAGAGTAATTT	7	0.17500000000000002	No Hit
GACTCTTAGCGGTGGATCACTCGGCTCGTGCGAGATCGGAAGAGCACACGT	7	0.17500000000000002	No Hit
TGCTGTGATGACTATCTTAGGACACCTTTGGATTAACCGTGAAATAAAACA	7	0.17500000000000002	No Hit
ACAAATGATGAATCCTTTGACGGGCGGGCAGAATCTCTGTGCTGATTGTCA	7	0.17500000000000002	No Hit
TAATACTGCCGGGTAATGATGGAAGATCGGAAGAGCACACGTCTGAACTCC	7	0.17500000000000002	Illumina Multiplexing PCR Primer 2.01 (100% over 28bp)
AAATATGATGAGCATACAAAATGTTGATTACAGTGAAAAACAGTAAGATCC	7	0.17500000000000002	No Hit
CGACTCTTAGCGGTGGATCACTCGGCTCGTGAGATCGGAAGAGCACACGTC	7	0.17500000000000002	No Hit
GCACAATGATGACTTATGTTTTTGCCGTTTACCCAGCTGAGGGTTTCTTTG	7	0.17500000000000002	No Hit
CACAAGTGATGACTAACAAAGGGACTTAATACTGAAACATCTGGGATTGAA	7	0.17500000000000002	No Hit
CTTCGGTGATGACACGATGACGAGTCAGAAAGGCCACGTCTTGCTCTTGGG	7	0.17500000000000002	No Hit
TCAGTGCACTACAGAACTTTGTAGATCGGAAGAGCACACGTCTGAACTCCA	6	0.15	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
GACTCTTAGCGGTGGATCACTCGGCTCGTGCGTCGATGAAGAACGCAGCTA	6	0.15	No Hit
TACTTGGATGGGAGACCGCCTGGGAATACCGGGTGCTGTAGGCTTTAGATC	6	0.15	No Hit
TCAAATGATGATATCACCTTAAATAGCTGGAATTACCGGCAGATTGCTAGT	6	0.15	No Hit
ACTTAATGATGATTGTCTATTTGATTGCTTGATGCAATGTGAAAACCACAT	6	0.15	No Hit
GTGGAATGATGACATTCTCCGGAATCGCTGTACTGTCTTGATGAAAGTACT	6	0.15	No Hit
NGCTGTGATGAGGACACCAGTAGGAAGTGCCGTCAGAGTCGATAACTGACG	6	0.15	No Hit
CACAAATGATGAATCCTTTGACGGGCGGGCAGAATCTCTGTGCTGATTGTC	6	0.15	No Hit
CCCCTGTGATGTTCAGCACTGGGCTCTGATCACCCTGAGGACATGGAGCTA	6	0.15	No Hit
CGCGACCTCAGATCAGACGTGGCGACCCGCTGAATTTAGATCGGAAGAGCA	6	0.15	No Hit
CGTTGGATCGATGTGGTGCTGCCGCGTTCTCTTCGGGCCGGGGCCTAAGCC	6	0.15	No Hit
GATCTATGATGATCCTATCCCGAATCTGAATTCCTGTTGAAAAAATATAAT	6	0.15	No Hit
CGCGACCTCAGATCAGAAGATCGGAAGAGCACACGTCTGAACTCCAGTCAC	5	0.125	Illumina Multiplexing PCR Primer 2.01 (100% over 34bp)
GGAGATGATGAACTGACTCAGCCTGACCTTCAGAAATGGAGGTGATAAAAC	5	0.125	No Hit
TGAGGTAGTAGGTTGTATAGTTAGATCGGAAGAGCACACGTCTGAACTCCA	5	0.125	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
TTTGGCAATGGTAGAACTCACACCAGATCGGAAGAGCACACGTCTGAACTC	5	0.125	Illumina Multiplexing PCR Primer 2.01 (100% over 27bp)
TTCCGCTGATGACGTCTGTTAGTGCCACGTGTCTGGGCCACTGAGAGCACA	5	0.125	No Hit
TGAGGTAGTAGATTGTATAGTTAGATCGGAAGAGCACACGTCTGAACTCCA	5	0.125	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
CACCCGTAGAACCGACCTTGCGAGATCGGAAGAGCACACGTCTGAACTCCA	5	0.125	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
TGAAGTGATGATTCACAATCATGTCTCTTCTCTGAAGAATTCTTGGAGATT	5	0.125	No Hit
NGGAATGATGACATTCTCCGGAATCGCTGTACTGTCTTGATGAAAGTACTT	5	0.125	No Hit
CGGGTAGGATGACTCACATTTACAGACTAGAGTTTCTGATGTCCACGTGAC	5	0.125	No Hit
TGTCCCATTTGACTACCACATTTCCTACTGATAATAGATTTTGGTGGCTAT	5	0.125	No Hit
CGCGACCTCAGATCAGACGTAGATCGGAAGAGCACACGTCTGAACTCCAGT	5	0.125	Illumina Multiplexing PCR Primer 2.01 (100% over 31bp)
>>END_MODULE
>>Adapter Content	fail
#Position	Illumina Universal Adapter	Illumina Small RNA 3' Adapter	Illumina Small RNA 5' Adapter	Nextera Transposase Sequence	SOLID Small RNA Adapter
1	8.6	0.0	0.0	0.0	0.0
2	8.775	0.0	0.0	0.0	0.0
3	8.8	0.0	0.0	0.0	0.0
4	8.8	0.0	0.0	0.0	0.0
5	8.8	0.0	0.0	0.0	0.0
6	8.875	0.0	0.0	0.0	0.0
7	8.9	0.0	0.0	0.0	0.0
8	9.125	0.0	0.0	0.0	0.0
9	9.3	0.0	0.0	0.0	0.0
10	11.225	0.0	0.0	0.0	0.0
11	12.075	0.0	0.0	0.0	0.0
12	12.45	0.0	0.0	0.0	0.0
13	12.825	0.0	0.0	0.0	0.0
14	13.375	0.0	0.0	0.0	0.0
15	13.775	0.0	0.0	0.0	0.0
16	14.225	0.0	0.0	0.0	0.0
17	14.525	0.0	0.0	0.0	0.0
18	15.0	0.0	0.0	0.0	0.0
19	15.675	0.0	0.0	0.0	0.0
20	16.275	0.0	0.0	0.0	0.0
21	17.025	0.0	0.0	0.0	0.0
22	18.475	0.0	0.0	0.0	0.0
23	22.15	0.0	0.0	0.0	0.0
24	24.275	0.0	0.0	0.0	0.0
25	25.6	0.0	0.0	0.0	0.0
26	26.025	0.0	0.0	0.0	0.0
27	26.35	0.0	0.0	0.0	0.0
28	26.6	0.0	0.0	0.0	0.0
29	27.2	0.0	0.0	0.0	0.0
30	27.425	0.0	0.0	0.0	0.0
31	28.1	0.0	0.0	0.0	0.0
32	28.8	0.0	0.0	0.0	0.0
33	29.45	0.0	0.0	0.0	0.0
34	29.775	0.0	0.0	0.0	0.0
35	29.925	0.0	0.0	0.0	0.0
36	30.225	0.0	0.0	0.0	0.0
37	30.65	0.0	0.0	0.0	0.0
38	31.125	0.0	0.0	0.0	0.0
39	31.3	0.0	0.0	0.0	0.0
>>END_MODULE
>>Kmer Content	fail
#Sequence	Count	PValue	Obs/Exp Max	Max Obs/Exp Position
TGCACTG	55	0.0	49.888893	1
TAGCTTA	15	0.007837298	49.88889	1
TGAGGTA	15	0.007837298	49.88889	1
CGACAAC	55	0.0	44.9	40
TAGAAAT	55	0.0	44.9	15
GAGGTAG	20	6.2916183E-4	44.9	2
GTTTGTG	20	6.2916183E-4	44.9	11
GTATGCC	30	1.6973481E-6	44.9	41
AAACTAG	55	0.0	44.9	11
GAAACTA	55	0.0	44.9	10
TAGGAAG	55	0.0	44.9	21
GAAACGA	55	0.0	44.9	36
AGTAGTT	20	6.2916183E-4	44.9	7
GAAATAG	55	0.0	44.9	17
AACGACA	55	0.0	44.9	38
ATGATGT	20	6.2916183E-4	44.9	6
AACTGAC	55	0.0	44.9	44
CAACTGA	55	0.0	44.9	43
TATGCCG	30	1.6973481E-6	44.9	42
GAAGTGC	55	0.0	44.9	24
>>END_MODULE
Rejected 2633546 READS because READLEN < 1
Read 2633546 spots for ERR4018149.sra
Written 2633546 spots for ERR4018149.sra
Rejected 2633546 READS because READLEN < 1
Read 2633546 spots for ERR4018149.sra
Written 2633546 spots for ERR4018149.sra
Rejected 2633546 READS because READLEN < 1
Read 2633546 spots for ERR4018149.sra
Written 2633546 spots for ERR4018149.sra
Rejected 2633546 READS because READLEN < 1
Read 2633546 spots for ERR4018149.sra
Written 2633546 spots for ERR4018149.sra
Rejected 2633546 READS because READLEN < 1
Read 2633546 spots for ERR4018149.sra
Written 2633546 spots for ERR4018149.sra
Rejected 2633546 READS because READLEN < 1
Read 2633546 spots for ERR4018149.sra
Written 2633546 spots for ERR4018149.sra
Rejected 2633546 READS because READLEN < 1
Read 2633546 spots for ERR4018149.sra
Written 2633546 spots for ERR4018149.sra
Rejected 2633550 READS because READLEN < 1
Read 2633550 spots for ERR4018149.sra
Written 2633550 spots for ERR4018149.sra
SRR ids: ['ERR4018149.sra']
extra args: ['--split-files', '--defline-qual', '+']
tempdir: /tmp/pfd_7mgg1h8a
ERR4018149.sra spots: 21068372
blocks: [[1, 2633546], [2633547, 5267092], [5267093, 7900638], [7900639, 10534184], [10534185, 13167730], [13167731, 15801276], [15801277, 18434822], [18434823, 21068372]]
ERR4018149 file size 2982188
ERR4018149 completed basic pipeline successfully
skewer v0.2.2 [April 4, 2016]
COMMAND LINE:	/usr/local/bin/skewer -f sanger -l 18 -q 10 -k inf -t 8 -o ERR4018149 ERR4018149_1.fastq
Input file:	ERR4018149_1.fastq
trimmed:	ERR4018149-trimmed.fastq

Parameters used:
-- 3' end adapter sequence (-x):	AGATCGGAAGAGCACACGTCTGAACTCCAGTCAC
-- maximum error ratio allowed (-r):	0.100
-- maximum indel error ratio allowed (-d):	0.030
-- end quality threshold (-q):		10
-- minimum read length allowed after trimming (-l):	18
-- file format (-f):		Sanger/Illumina 1.8+ FASTQ 
-- minimum overlap length for adapter detection (-k):	inf
-- number of concurrent threads (-t):	8
Sun Dec 28 15:05:53 2025 >> started

Sun Dec 28 15:05:59 2025 >> done (6.229s)
21068372 reads processed; of these:
 1338329 ( 6.35%) short reads filtered out after trimming by size control
 2097227 ( 9.95%) empty reads filtered out after trimming by size control
17632816 (83.69%) reads available; of these:
  203323 ( 1.15%) trimmed reads available after processing
17429493 (98.85%) untrimmed reads available after processing

Length distribution of reads after trimming:
length	count	percentage
 18	  109953	  0.62%
 19	   93275	  0.53%
 20	       0	  0.00%
 21	       0	  0.00%
 22	       0	  0.00%
 23	       0	  0.00%
 24	       0	  0.00%
 25	       0	  0.00%
 26	       0	  0.00%
 27	       0	  0.00%
 28	       0	  0.00%
 29	       0	  0.00%
 30	       0	  0.00%
 31	       0	  0.00%
 32	       0	  0.00%
 33	       0	  0.00%
 34	       0	  0.00%
 35	       0	  0.00%
 36	       0	  0.00%
 37	       0	  0.00%
 38	       0	  0.00%
 39	       0	  0.00%
 40	       0	  0.00%
 41	       0	  0.00%
 42	       0	  0.00%
 43	       0	  0.00%
 44	       0	  0.00%
 45	       0	  0.00%
 46	       0	  0.00%
 47	       0	  0.00%
 48	       0	  0.00%
 49	      31	  0.00%
 50	      64	  0.00%
 51	17429493	 98.85%
17632816 reads passed initial QC


criterion=sequence-density
sequence-density=17.87
sequence-density-rank=1
fanout-score=33.21
fanout-score-rank=1
prefix-density=18.58
prefix-fanout=31.9
sequence=AGATCGGAAGAGCACACGTCTGAACTCCAGTCACCTTGTAATCTCGTATGCCGTCTTC


criterion=fanout-score
sequence-density=17.87
sequence-density-rank=1
fanout-score=33.21
fanout-score-rank=1
prefix-density=18.58
prefix-fanout=31.9
sequence=AGATCGGAAGAGCACACGTCTGAACTCCAGTCACCTTGTAATCTCGTATGCCGTCTTC
Potential 3prime adapter identified. Now checking if in reference sequence
Warning: gzbuffer added in zlib v1.2.3.5. Unable to change buffer size from default of 8192.
1 reads; of these:
  1 (100.00%) were unpaired; of these:
    1 (100.00%) aligned 0 times
    0 (0.00%) aligned exactly 1 time
    0 (0.00%) aligned >1 times
0.00% overall alignment rate
Adapter seq not found in reference. Now shuffling file before clipping
skewer v0.2.2 [April 4, 2016]
COMMAND LINE:	/usr/local/bin/skewer -f sanger -l 18 -t 8 -x AGATCGGAAGAGCACACGTCTGAACTCCAGTCACCTTGTAATCTCGTATGCCGTCTTC -o ERR4018149 -
Input file:	STDIN
trimmed:	ERR4018149-trimmed.fastq

Parameters used:
-- 3' end adapter sequence (-x):	AGATCGGAAGAGCACACGTCTGAACTCCAGTCACCTTGTAATCTCGTATGCCGTCTTC
-- maximum error ratio allowed (-r):	0.100
-- maximum indel error ratio allowed (-d):	0.030
-- minimum read length allowed after trimming (-l):	18
-- file format (-f):		Sanger/Illumina 1.8+ FASTQ 
-- minimum overlap length for adapter detection (-k):	3
-- number of concurrent threads (-t):	8
Sun Dec 28 15:06:11 2025 >> started

Sun Dec 28 15:06:20 2025 >> done (8.811s)
15673614 reads processed; of these:
    5665 ( 0.04%) short reads filtered out after trimming by size control
    3201 ( 0.02%) empty reads filtered out after trimming by size control
15664748 (99.94%) reads available; of these:
 4172612 (26.64%) trimmed reads available after processing
11492136 (73.36%) untrimmed reads available after processing

Length distribution of reads after trimming:
length	count	percentage
 18	   95343	  0.61%
 19	   84823	  0.54%
 20	  139586	  0.89%
 21	  273800	  1.75%
 22	  662386	  4.23%
 23	  406753	  2.60%
 24	  264179	  1.69%
 25	   90399	  0.58%
 26	   64692	  0.41%
 27	   64958	  0.41%
 28	   86208	  0.55%
 29	   55688	  0.36%
 30	  103327	  0.66%
 31	  156174	  1.00%
 32	  127526	  0.81%
 33	   65254	  0.42%
 34	   60878	  0.39%
 35	   48222	  0.31%
 36	   69187	  0.44%
 37	   59791	  0.38%
 38	   39596	  0.25%
 39	   43059	  0.27%
 40	   28190	  0.18%
 41	   27820	  0.18%
 42	   52542	  0.34%
 43	  102702	  0.66%
 44	  156478	  1.00%
 45	  528640	  3.37%
 46	   89327	  0.57%
 47	   40233	  0.26%
 48	  261412	  1.67%
 49	      17	  0.00%
 50	      44	  0.00%
 51	11315514	 72.24%


criterion=sequence-density
sequence-density=4.14
sequence-density-rank=1
fanout-score=1.00
fanout-score-rank=14
prefix-density=0.01
prefix-fanout=1.0
sequence=GTGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAACTCTCTGCTGATC


criterion=fanout-score
sequence-density=0.50
sequence-density-rank=11
fanout-score=40.47
fanout-score-rank=1
prefix-density=6.25
prefix-fanout=3.2
sequence=GTGATGACTTTC
                                 Started job on |	Dec 28 15:06:29
                             Started mapping on |	Dec 28 15:06:29
                                    Finished on |	Dec 28 15:06:59
       Mapping speed, Million of reads per hour |	2114.87

                          Number of input reads |	17623950
                      Average input read length |	46
                                    UNIQUE READS:
                   Uniquely mapped reads number |	12157933
                        Uniquely mapped reads % |	68.99%
                          Average mapped length |	47.39
                       Number of splices: Total |	11686
            Number of splices: Annotated (sjdb) |	3015
                       Number of splices: GT/AG |	10853
                       Number of splices: GC/AG |	184
                       Number of splices: AT/AC |	3
               Number of splices: Non-canonical |	646
                      Mismatch rate per base, % |	0.36%
                         Deletion rate per base |	0.00%
                        Deletion average length |	1.09
                        Insertion rate per base |	0.00%
                       Insertion average length |	1.20
                             MULTI-MAPPING READS:
        Number of reads mapped to multiple loci |	3355319
             % of reads mapped to multiple loci |	19.04%
        Number of reads mapped to too many loci |	353658
             % of reads mapped to too many loci |	2.01%
                                  UNMAPPED READS:
       % of reads unmapped: too many mismatches |	0.00%
                 % of reads unmapped: too short |	9.88%
                     % of reads unmapped: other |	0.09%
                                  CHIMERIC READS:
                       Number of chimeric reads |	0
                            % of chimeric reads |	0.00%
N_unmapped	2110698	2110698	2110698
N_multimapping	3355319	3355319	3355319
N_noFeature	359569	474293	11641528
N_ambiguous	1067224	665512	37
UnstrandedReadsAssigned:10731140 PositiveStrandReadsAssigned:11018128 NegativeStrandReadsAssigned:516368
Dataset is classified positive stranded
MeadianReadLen=51 20thPercentileLength=42 echo kmer=37
ERR4018149 Starting Kallisto single end mapping to ensembl reference transcriptome. kmer=31

[quant] fragment length distribution is truncated gaussian with mean = 100, sd = 20
[index] k-mer length: 31
[index] number of targets: 31,715
[index] number of k-mers: 49,915,185
[index] number of equivalence classes: 70,351
[quant] running in single-end mode
[quant] will process file 1: ERR4018149-trimmed.fastq
[quant] finding pseudoalignments for the reads ... done
[quant] processed 17,623,950 reads, 125,729 reads pseudoaligned
[   em] quantifying the abundances ... done
[   em] the Expectation-Maximization algorithm ran for 620 rounds

  31716 ERR4018149.ke.tsv
  32883 ERR4018149.se.tsv
  64599 total
==> ERR4018149.ke.tsv <==
target_id	length	eff_length	est_counts	tpm
ENSRNOT00000047550.4	955	856	87	556.523
ENSRNOT00000040993.4	1039	940	158	920.38
ENSRNOT00000050156.3	1545	1446	782	2961.26
ENSRNOT00000043693.3	684	585	275	2574.04
ENSRNOT00000046201.3	204	105	165	8604.64
ENSRNOT00000046108.3	681	582	85	799.712
ENSRNOT00000049683.3	784	685	90	719.432
ENSRNOT00000041241.3	348	249	57	1253.47
ENSRNOT00000044582.3	297	198	71	1963.5

==> ERR4018149.se.tsv <==
ENSRNOG00000046319	0
ENSRNOG00000047964	0
ENSRNOG00000050370	0
ENSRNOG00000032365	0
ENSRNOG00000040300	0
ENSRNOG00000058808	0
ENSRNOG00000061316	0
ENSRNOG00000050129	1
ENSRNOG00000040316	0
ENSRNOG00000023659	0
ERR4018149 completed mapping pipeline successfully
