Starting /dee2/code/volunteer_pipeline.sh ERR4018146
    current disk space = 370819820
    free memory = 197288588 
ERR4018146 SRAfilesize
979dcc79168b9e37223053766b80e48d  ERR4018146.sra
ERR4018146.sra file validated
ERR4018146 is single end
ERR4018146 is conventional basespace
ERR4018146 read1 length is 51 nt
##FastQC	0.11.5
>>Basic Statistics	pass
#Measure	Value
Filename	ERR4018146_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.84025	35.0	35.0	35.0	27.0	35.0
2	33.95225	35.0	35.0	35.0	31.0	35.0
3	34.186	35.0	35.0	35.0	32.0	35.0
4	34.355	35.0	35.0	35.0	33.0	35.0
5	33.72775	35.0	35.0	35.0	31.0	35.0
6	38.34075	39.0	39.0	40.0	35.0	40.0
7	38.961	40.0	39.0	40.0	38.0	40.0
8	39.0235	40.0	39.0	40.0	38.0	40.0
9	38.922	40.0	39.0	40.0	37.0	40.0
10	39.1425	40.0	39.0	40.0	38.0	40.0
11	39.08725	40.0	39.0	40.0	38.0	40.0
12	39.1095	40.0	39.0	40.0	38.0	40.0
13	39.08875	40.0	39.0	40.0	38.0	40.0
14	39.02	40.0	39.0	40.0	38.0	40.0
15	39.09675	40.0	39.0	40.0	38.0	40.0
16	39.074	40.0	39.0	40.0	38.0	40.0
17	39.22875	40.0	39.0	40.0	38.0	40.0
18	39.1245	40.0	39.0	40.0	38.0	40.0
19	39.1875	40.0	39.0	40.0	38.0	40.0
20	39.1825	40.0	39.0	40.0	38.0	40.0
21	39.22825	40.0	39.0	40.0	38.0	40.0
22	39.17675	40.0	39.0	40.0	38.0	40.0
23	39.1195	40.0	39.0	40.0	38.0	40.0
24	39.14775	40.0	39.0	40.0	38.0	40.0
25	39.155	40.0	39.0	40.0	38.0	40.0
26	39.1405	40.0	39.0	40.0	38.0	40.0
27	39.21575	40.0	39.0	40.0	38.0	40.0
28	39.18675	40.0	39.0	40.0	38.0	40.0
29	39.1365	40.0	39.0	40.0	38.0	40.0
30	39.19375	40.0	39.0	40.0	38.0	40.0
31	38.97925	40.0	39.0	40.0	38.0	40.0
32	39.012	40.0	39.0	40.0	38.0	40.0
33	39.1495	40.0	39.0	40.0	38.0	40.0
34	39.089	40.0	39.0	40.0	38.0	40.0
35	39.10775	40.0	39.0	40.0	38.0	40.0
36	38.908	40.0	39.0	40.0	38.0	40.0
37	39.0265	40.0	39.0	40.0	38.0	40.0
38	39.01325	40.0	39.0	40.0	38.0	40.0
39	38.94175	40.0	39.0	40.0	38.0	40.0
40	38.9775	40.0	39.0	40.0	38.0	40.0
41	38.8855	40.0	39.0	40.0	37.0	40.0
42	38.857	40.0	39.0	40.0	38.0	40.0
43	38.991	40.0	39.0	40.0	38.0	40.0
44	38.94925	40.0	39.0	40.0	38.0	40.0
45	38.9115	40.0	39.0	40.0	38.0	40.0
46	39.0055	40.0	39.0	40.0	38.0	40.0
47	39.014	40.0	39.0	40.0	38.0	40.0
48	38.9995	40.0	39.0	40.0	38.0	40.0
49	38.726	40.0	39.0	40.0	36.0	40.0
50	38.81375	40.0	39.0	40.0	37.0	40.0
51	38.2585	39.0	39.0	40.0	35.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	1.0
23	0.0
24	1.0
25	1.0
26	4.0
27	2.0
28	2.0
29	8.0
30	14.0
31	25.0
32	23.0
33	34.0
34	67.0
35	79.0
36	126.0
37	297.0
38	1153.0
39	2162.0
>>END_MODULE
>>Per base sequence content	fail
#Base	G	A	T	C
1	22.192078133478024	20.184481823114485	36.299511665762346	21.323928377645142
2	45.95	13.725000000000001	30.275000000000002	10.05
3	22.675	36.65	11.675	28.999999999999996
4	16.475	38.275	20.225	25.025
5	20.025000000000002	37.4	15.5	27.075
6	26.3	19.7	43.55	10.45
7	55.65	7.199999999999999	27.750000000000004	9.4
8	27.6	47.949999999999996	13.0	11.450000000000001
9	10.575	40.75	39.525	9.15
10	50.975	14.35	25.900000000000002	8.774999999999999
11	40.075	43.3	8.725	7.9
12	28.199999999999996	30.55	17.275	23.974999999999998
13	17.95	30.575000000000003	19.900000000000002	31.574999999999996
14	9.925	38.324999999999996	33.45	18.3
15	17.325	21.125	34.0	27.55
16	14.924999999999999	34.449999999999996	33.2	17.424999999999997
17	22.975	23.875	16.5	36.65
18	19.15	33.175	29.75	17.925
19	15.225	31.85	34.575	18.35
20	15.725	25.624999999999996	32.05	26.6
21	23.1	13.600000000000001	52.75	10.549999999999999
22	31.55	23.65	28.475	16.325
23	26.150000000000002	33.125	10.325	30.4
24	32.2	30.075000000000003	18.15	19.575
25	31.075000000000003	26.55	13.350000000000001	29.025000000000002
26	23.674999999999997	22.2	25.974999999999998	28.15
27	29.075	20.3	13.025	37.6
28	19.275000000000002	18.875	27.075	34.775
29	28.000000000000004	22.075	27.125	22.8
30	38.625	15.475	25.45	20.45
31	24.575	30.7	25.374999999999996	19.35
32	27.400000000000002	23.275000000000002	20.175	29.15
33	13.775	32.2	18.3	35.725
34	14.000000000000002	22.075	27.450000000000003	36.475
35	35.675000000000004	15.925	24.15	24.25
36	24.325	48.275	9.325	18.075
37	11.85	43.6	25.674999999999997	18.875
38	19.6	33.475	10.875	36.05
39	21.8	33.275	26.924999999999997	18.0
40	39.675	14.499999999999998	19.75	26.075
41	32.425	28.299999999999997	22.3	16.975
42	18.55	36.525	32.7	12.225
43	19.075	27.775	23.674999999999997	29.475
44	25.775	34.325	25.95	13.950000000000001
45	20.925	41.949999999999996	17.325	19.8
46	25.1	29.4	22.225	23.275000000000002
47	12.525	28.325	37.525	21.625
48	24.15	41.075	19.55	15.225
49	12.9	35.825	35.699999999999996	15.575
50	25.525	20.875	18.2	35.4
51	25.3	21.3	23.200000000000003	30.2
>>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.5
12	1.0
13	0.5
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	2.5
26	3.5
27	2.0
28	18.0
29	34.0
30	38.5
31	43.0
32	32.5
33	22.0
34	51.5
35	81.0
36	278.0
37	475.0
38	430.0
39	385.0
40	363.5
41	342.0
42	231.5
43	121.0
44	124.5
45	128.0
46	267.5
47	407.0
48	333.0
49	259.0
50	473.0
51	687.0
52	530.0
53	373.0
54	309.5
55	246.0
56	197.5
57	149.0
58	144.0
59	139.0
60	92.5
61	46.0
62	29.5
63	13.0
64	8.5
65	4.0
66	8.5
67	13.0
68	14.0
69	15.0
70	12.0
71	9.0
72	5.0
73	1.0
74	0.5
75	0.0
76	0.0
77	0.0
78	0.0
79	0.0
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	7.85
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	33.475
#Duplication Level	Percentage of deduplicated	Percentage of total
1	78.79014189693801	26.375
2	8.812546676624347	5.8999999999999995
3	2.9873039581777445	3.0
4	1.8670649738610903	2.5
5	1.194921583271098	2.0
6	0.9708737864077669	1.95
7	0.597460791635549	1.4000000000000001
8	0.7468259895444361	2.0
9	0.07468259895444362	0.22499999999999998
>10	3.286034353995519	20.45
>50	0.22404779686333084	5.525
>100	0.4480955937266617	28.675
>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
AGATCGGAAGAGCACACGTCTGAACTCCAGTCACGATCAGATCTCGTATGC	342	8.55	TruSeq Adapter, Index 9 (100% over 50bp)
TGCACTGATGAAACTAGAAATAGGAAGTGCCGTCCGAAACGACAACTGACG	205	5.125	No Hit
TGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAACT	202	5.050000000000001	No Hit
GTAAATGATGAGATTTCCATTGGTCCGTGTTTCTGAACTCCATGATTTTTG	141	3.5249999999999995	No Hit
GGAAGAGCACACGTCTGAACTCCAGTCACGATCAGATCTCGTATGCCGTCT	140	3.5000000000000004	TruSeq Adapter, Index 9 (100% over 51bp)
GTGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAAC	117	2.9250000000000003	No Hit
TTGCTGTGATGACTATCTTAGGACACCTTTGGATTAACCGTGAAATAAAAC	96	2.4	No Hit
CTCGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGTAGATCG	72	1.7999999999999998	No Hit
TTGCAATGATGTGAATCTTTTGCTGAAGTGACCTTGAAGTGTAATTCCTGA	53	1.325	No Hit
CGCGACCTCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACGATCAGAT	44	1.0999999999999999	TruSeq Adapter, Index 9 (100% over 41bp)
TACTGTGATGATAACATAGTTCAGCAGACGTGACTCTGGTGAACAATCATG	43	1.075	No Hit
CACAATGATGACTTATGTTTTTGCCGTTTACCCAGCTGAGGGTTTCTTTGA	42	1.05	No Hit
ATACATGATGATCTCACTCCAACTTGAACTCTCTCACTGATTACTTGATGA	41	1.0250000000000001	No Hit
TGGAATGATGACATTCTCCGGAATCGCTGTACTGTCTTGATGAAAGTACTT	33	0.8250000000000001	No Hit
TTGAGTGATGACTTGTATTGTCGGATACCCCTTCACCCCGTTCATGGGTGA	31	0.775	No Hit
TGCTGTGATGAGGACACCAGTAGGAAGTGCCGTCAGAGTCGATAACTGACG	30	0.75	No Hit
CTGCGATGATGTCATTTCTTAGGACACCTTTGGATTAATCATGAAAATAGC	28	0.7000000000000001	No Hit
GCTAATGATGATACCGCACTAACTGAGCAGTTAGTAGTTGGTCCTTTGATT	27	0.675	No Hit
GGTCTGTGATGAAACCCATTATCAGTGGACATCTATGGATGAGAAATGCGG	25	0.625	No Hit
GGTCAATGATGTGTTGGCATGTATTATCTGAATTTGCTGATGTGTCCTAAC	21	0.525	No Hit
CGACTCTTAGCGGTGGATCACTCGGCTCGTGAGATCGGAAGAGCACACGTC	20	0.5	No Hit
NGATCGGAAGAGCACACGTCTGAACTCCAGTCACGATCAGATCTCGTATGC	20	0.5	TruSeq Adapter, Index 9 (100% over 50bp)
TAGCTTATCAGACTGATGTTGACAGATCGGAAGAGCACACGTCTGAACTCC	19	0.475	Illumina Multiplexing PCR Primer 2.01 (100% over 28bp)
TCGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGTAGATCGG	19	0.475	No Hit
GTTCTGTGATGAGACTCTGGAATGTGCGCTGGGCACAGCGCCCGAGTGAGA	18	0.44999999999999996	No Hit
TGAGGTAGTAGTTTGTGCTGTTAGATCGGAAGAGCACACGTCTGAACTCCA	17	0.42500000000000004	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
CTTCGGTGATGACACGATGACGAGTCAGAAAGGCCACGTCTTGCTCTTGGG	17	0.42500000000000004	No Hit
NGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAACT	16	0.4	No Hit
CTCCCATGATGTATAGCACTGGGCTCTGACCACCATGGTGACGTGTTGTGC	15	0.375	No Hit
CGCGACCTTAGATCGGAAGAGCACACGTCTGAACTCCAGTCACGATCAGAT	15	0.375	TruSeq Adapter, Index 9 (100% over 41bp)
CTTAATGATGATTGTCTATTTGATTGCTTGATGCAATGTGAAAACCACATT	15	0.375	No Hit
NGCACTGATGAAACTAGAAATAGGAAGTGCCGTCCGAAACGACAACTGACG	14	0.35000000000000003	No Hit
CACCCTGATGAACTGAATACCGCCCCAGTCTGATAGATGTGGAGAAAGGTA	14	0.35000000000000003	No Hit
CGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGTAGATCGGA	14	0.35000000000000003	No Hit
CTTCTGTGATGACTTTACCCAAATGACTTTCGTTCTTCTGAGTTTGCTGAA	13	0.325	No Hit
TTCAAGTAATCCAGGATAGGCTAGATCGGAAGAGCACACGTCTGAACTCCA	13	0.325	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
CTCGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGCAGATCG	13	0.325	No Hit
ATCGATGTGGTGCTGCCGCGTTCTCTTCGGGCCGGGGCCTAAGCCGCGCCA	13	0.325	No Hit
AGTAAATGATGAGATTTCCATTGGTCCGTGTTTCTGAACTCCATGATTTTT	13	0.325	No Hit
AACAGTGATGACCACAACATGCCATCTGAGTTGCTTTGCTGAAATCCAGAG	12	0.3	No Hit
TATCGATGATGTCCTAAAAGTAAATGTCTGAACATATGAATGCTATACTGA	12	0.3	No Hit
NGAAGAGCACACGTCTGAACTCCAGTCACGATCAGATCTCGTATGCCGTCT	12	0.3	TruSeq Adapter, Index 9 (98% over 51bp)
CACCCTGATGAACTGAATACCGCCCCAGTCTGATAGCTGTGGAGAAAGGTA	12	0.3	No Hit
CTACAATGATGGTAATATTTTTCGTCAACAGAGTTCACCTAGTGAGTGTTG	12	0.3	No Hit
AACCCGTAGATCCGATCTTGTAGATCGGAAGAGCACACGTCTGAACTCCAG	11	0.27499999999999997	Illumina Multiplexing PCR Primer 2.01 (100% over 30bp)
TTTGAGTGATGACTTGTATTGTCGGATACCCCTTCACCCCGTTCATGGGTG	11	0.27499999999999997	No Hit
AAATATGATGAGCATACAAAATGTTGATTACAGTGAAAAACAGTAAGATCC	11	0.27499999999999997	No Hit
NTGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAAC	11	0.27499999999999997	No Hit
TAGCTTATCAGACTGATGTTGACTAGATCGGAAGAGCACACGTCTGAACTC	11	0.27499999999999997	Illumina Multiplexing PCR Primer 2.01 (100% over 27bp)
NTAAATGATGAGATTTCCATTGGTCCGTGTTTCTGAACTCCATGATTTTTG	10	0.25	No Hit
TGAGGTAGTAGGTTGTATGGTTAGATCGGAAGAGCACACGTCTGAACTCCA	10	0.25	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
CGCGACCTCAGATCAGACGTAGATCGGAAGAGCACACGTCTGAACTCCAGT	10	0.25	Illumina Multiplexing PCR Primer 2.01 (100% over 31bp)
CGACTCTTAGCGGTGGATCACTCGGCTCGTGCGTCGATGAAGAACGCAGCT	10	0.25	No Hit
TCAGTGCACTACAGAACTTTGTAGATCGGAAGAGCACACGTCTGAACTCCA	9	0.22499999999999998	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
ACAAGTGATGACTAACAAAGGGACTTAATACTGAAACATCTGGGATTGAAA	8	0.2	No Hit
GGTCAATGATGAGCTGGCATGTATTCTGAATCTAAAGTTGATTAATAACCA	8	0.2	No Hit
CGCGACCTCAGATCAGACGTGGCGACCCGCTGAATTTAAGCATATTAGTCA	8	0.2	No Hit
NTCGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGTAGATCG	8	0.2	No Hit
TGCTGTGATGACTATCTTAGGACACCTTTGGATTAACCGTGAAATAAAACA	8	0.2	No Hit
AGATCGGAAGAGCACACGTCTAAACTCCAGTCACGATCAGATCTCGTATGC	8	0.2	TruSeq Adapter, Index 9 (98% over 50bp)
ACTTAATGATGATTGTCTATTTGATTGCTTGATGCAATGTGAAAACCACAT	8	0.2	No Hit
GAAATATGATGAGCATACAAAATGTTGATTACAGTGAAAAACAGTAAGATC	8	0.2	No Hit
CTGCAGTGATGACTTTCTTAGGACACCTTTGGATTTACCGTGAAAATAAGT	8	0.2	No Hit
CCGCGACCTCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACGATCAGA	8	0.2	TruSeq Adapter, Index 9 (100% over 40bp)
NTGCTGTGATGACTATCTTAGGACACCTTTGGATTAACCGTGAAATAAAAC	7	0.17500000000000002	No Hit
GCGGATGATGACACCTGGGTATGCTGCACATTCCCGACTGCGTCGTGGGGA	7	0.17500000000000002	No Hit
TTCTGTGATGAGACTCTGGAATGTGCGCTGGGCACAGCGCCCGAGTGAGAC	7	0.17500000000000002	No Hit
ATCTATGATGATCCTATCCCGAATCTGAATTCCTGTTGAAAAAATATAATT	7	0.17500000000000002	No Hit
CACAAGTGATGACTAACAAAGGGACTTAATACTGAAACATCTGGGATTGAA	7	0.17500000000000002	No Hit
TAGGGTGATGAAAAAGTATCCTTAGGCGTGGTTCTGGCCGTCTTGGTCACC	7	0.17500000000000002	No Hit
ACCCCTGTGATGTTCAGCACTGGGCTCTGATCACCCTGAGGACATGGAGCT	7	0.17500000000000002	No Hit
TATGGCACTGGTAGAATTCACTAGATCGGAAGAGCACACGTCTGAACTCCA	7	0.17500000000000002	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
GTGCAATGATGTGTGAATTTCTTCACCTGAAATCACCGTGAAGAGTAATTT	6	0.15	No Hit
TGGCGGTGATGACCCACTTGCCCTCACTGAGGACAAAGTTCGTTGGTGAGA	6	0.15	No Hit
NAGCTTATCAGACTGATGTTGACAGATCGGAAGAGCACACGTCTGAACTCC	6	0.15	Illumina Multiplexing PCR Primer 2.01 (100% over 28bp)
TTTGGCAATGGTAGAACTCACACCAGATCGGAAGAGCACACGTCTGAACTC	6	0.15	Illumina Multiplexing PCR Primer 2.01 (100% over 27bp)
ACAAATGATGAATCCTTTGACGGGCGGGCAGAATCTCTGTGCTGATTGTCA	6	0.15	No Hit
GTGGAATGATGACATTCTCCGGAATCGCTGTACTGTCTTGATGAAAGTACT	6	0.15	No Hit
CACAAATGATGAATCCTTTGACGGGCGGGCAGAATCTCTGTGCTGATTGTC	6	0.15	No Hit
TGCAATGATGAAGGTTATGTTGGTAGGGACAACTGACCTTGTTGATGAATA	6	0.15	No Hit
NTGCAATGATGTGAATCTTTTGCTGAAGTGACCTTGAAGTGTAATTCCTGA	6	0.15	No Hit
CGCGACCTCAGATCGGAAGAGCACACGTCTAAACTCCAGTCACGATCAGAT	6	0.15	TruSeq Adapter, Index 9 (97% over 41bp)
CGCGACCTCAGATCAGACGTGGCGACCCGCTGAATTTAGATCGGAAGAGCA	6	0.15	No Hit
CGTTGGATCGATGTGGTGCTGCCGCGTTCTCTTCGGGCCGGGGCCTAAGCC	6	0.15	No Hit
TGCAATGATGTGTGAATTTCTTCACCTGAAATCACCGTGAAGAGTAATTTA	6	0.15	No Hit
AGTCTGTGATGAGCAATAATCACCATCTTTCGTTTGAGTCTCACGACTGTG	5	0.125	No Hit
CGACTCTTAGCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACGATCAG	5	0.125	TruSeq Adapter, Index 9 (100% over 39bp)
AATACATGATGATCTCACTCCAACTTGAACTCTCTCACTGATTACTTGATG	5	0.125	No Hit
NGACTCTTAGCGGTGGATCACTCGGCTCGTGCGTCGATGAAGAACGCAGCT	5	0.125	No Hit
NGTCTGTGATGAAACCCATTATCAGTGGACATCTATGGATGAGAAATGCGG	5	0.125	No Hit
TTCTGTGATGACTTTACCCAAATGACTTTCGTTCTTCTGAGTTTGCTGAAG	5	0.125	No Hit
GACTCTTAGCGGTGGATCACTCGGCTCGTGCGAGATCGGAAGAGCACACGT	5	0.125	No Hit
CTTGATGATGAATGTATTATAAATGACAAGCATATGGCTGAACTTAGGTGA	5	0.125	No Hit
GAAGCTGATGATACCACTTCTTTCCCCATCAGATCGACCCTGTTGATCTCA	5	0.125	No Hit
CTACGGGGATGATTTAGAAGCTGAACTCTCTCTTTCTGATGGTTTAGTGGA	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)
NGCTGTGATGAGGACACCAGTAGGAAGTGCCGTCAGAGTCGATAACTGACG	5	0.125	No Hit
TCGGATCCGTCTGAGCTTGGCTAGATCGGAAGAGCACACGTCTGAACTCCA	5	0.125	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
GAGCAATGATGACCCCAACATGCCATCTGAGTTGCTTTGCTGAAATCCAGA	5	0.125	No Hit
GGAGATGATGAACTGACTCAACCTGACCTTCAGAAATGGAGGTGATAAAAC	5	0.125	No Hit
>>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	9.2	0.0	0.0	0.0	0.0
2	9.3	0.0	0.0	0.0	0.0
3	9.4	0.0	0.0	0.0	0.0
4	9.4	0.0	0.0	0.0	0.0
5	9.425	0.0	0.0	0.0	0.0
6	9.45	0.0	0.0	0.0	0.0
7	9.6	0.0	0.0	0.0	0.0
8	9.775	0.0	0.0	0.0	0.0
9	10.05	0.0	0.0	0.0	0.0
10	12.45	0.0	0.0	0.0	0.0
11	13.125	0.0	0.0	0.0	0.0
12	13.7	0.0	0.0	0.0	0.0
13	14.125	0.0	0.0	0.0	0.0
14	14.875	0.0	0.0	0.0	0.0
15	15.35	0.0	0.0	0.0	0.0
16	15.775	0.0	0.0	0.0	0.0
17	16.275	0.0	0.0	0.0	0.0
18	16.875	0.0	0.0	0.0	0.0
19	17.575	0.0	0.0	0.0	0.0
20	18.2	0.0	0.0	0.0	0.0
21	19.1	0.0	0.0	0.0	0.0
22	20.65	0.0	0.0	0.0	0.0
23	24.975	0.0	0.0	0.0	0.0
24	27.075	0.0	0.0	0.0	0.0
25	28.575	0.0	0.0	0.0	0.0
26	29.175	0.0	0.0	0.0	0.0
27	29.625	0.0	0.0	0.0	0.0
28	29.9	0.0	0.0	0.0	0.0
29	30.55	0.0	0.0	0.0	0.0
30	30.875	0.0	0.0	0.0	0.0
31	31.5	0.0	0.0	0.0	0.0
32	32.475	0.0	0.0	0.0	0.0
33	32.8	0.0	0.0	0.0	0.0
34	33.175	0.0	0.0	0.0	0.0
35	33.45	0.0	0.0	0.0	0.0
36	33.625	0.0	0.0	0.0	0.0
37	33.775	0.0	0.0	0.0	0.0
38	34.15	0.0	0.0	0.0	0.0
39	34.325	0.0	0.0	0.0	0.0
>>END_MODULE
>>Kmer Content	warn
#Sequence	Count	PValue	Obs/Exp Max	Max Obs/Exp Position
TGCACTG	15	0.008276992	49.219177	1
ATACATG	30	8.7783155E-7	49.219177	1
TGATGAT	40	4.432877E-9	44.9125	6
CCAACTT	30	1.694094E-6	44.9125	19
CATGATG	30	1.694094E-6	44.9125	4
GTATGCC	25	3.273556E-5	44.9125	41
AAACTAG	20	6.2829873E-4	44.9125	11
ATCTCAC	30	1.694094E-6	44.9125	11
GAAACTA	20	6.2829873E-4	44.9125	10
ACTGATT	30	1.694094E-6	44.9125	36
ACTGATG	20	6.2829873E-4	44.9125	4
TTACTTG	30	1.694094E-6	44.9125	41
ACATGAT	30	1.694094E-6	44.9125	3
ACTGACG	20	6.2829873E-4	44.9125	45
ACTTGAT	30	1.694094E-6	44.9125	43
ACTTGAA	30	1.694094E-6	44.9125	22
CAACTTG	30	1.694094E-6	44.9125	20
ATGATCT	30	1.694094E-6	44.9125	8
TGATCTC	30	1.694094E-6	44.9125	9
TATGCCG	25	3.273556E-5	44.9125	42
>>END_MODULE
Rejected 2301475 READS because READLEN < 1
Read 2301475 spots for ERR4018146.sra
Written 2301475 spots for ERR4018146.sra
Rejected 2301475 READS because READLEN < 1
Read 2301475 spots for ERR4018146.sra
Written 2301475 spots for ERR4018146.sra
Rejected 2301475 READS because READLEN < 1
Read 2301475 spots for ERR4018146.sra
Written 2301475 spots for ERR4018146.sra
Rejected 2301475 READS because READLEN < 1
Read 2301475 spots for ERR4018146.sra
Written 2301475 spots for ERR4018146.sra
Rejected 2301480 READS because READLEN < 1
Read 2301480 spots for ERR4018146.sra
Written 2301480 spots for ERR4018146.sra
Rejected 2301475 READS because READLEN < 1
Read 2301475 spots for ERR4018146.sra
Written 2301475 spots for ERR4018146.sra
Rejected 2301475 READS because READLEN < 1
Read 2301475 spots for ERR4018146.sra
Written 2301475 spots for ERR4018146.sra
Rejected 2301475 READS because READLEN < 1
Read 2301475 spots for ERR4018146.sra
Written 2301475 spots for ERR4018146.sra
SRR ids: ['ERR4018146.sra']
extra args: ['--split-files', '--defline-qual', '+']
tempdir: /tmp/pfd_t9kgj22a
ERR4018146.sra spots: 18411805
blocks: [[1, 2301475], [2301476, 4602950], [4602951, 6904425], [6904426, 9205900], [9205901, 11507375], [11507376, 13808850], [13808851, 16110325], [16110326, 18411805]]
ERR4018146 file size 2603420
ERR4018146 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 ERR4018146 ERR4018146_1.fastq
Input file:	ERR4018146_1.fastq
trimmed:	ERR4018146-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:04:11 2025 >> started

Sun Dec 28 15:04:16 2025 >> done (5.071s)
18411805 reads processed; of these:
 1432785 ( 7.78%) short reads filtered out after trimming by size control
 1879077 (10.21%) empty reads filtered out after trimming by size control
15099943 (82.01%) reads available; of these:
  199738 ( 1.32%) trimmed reads available after processing
14900205 (98.68%) untrimmed reads available after processing

Length distribution of reads after trimming:
length	count	percentage
 18	  108213	  0.72%
 19	   91445	  0.61%
 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	      24	  0.00%
 50	      56	  0.00%
 51	14900205	 98.68%
15099943 reads passed initial QC


criterion=sequence-density
sequence-density=21.28
sequence-density-rank=1
fanout-score=30.92
fanout-score-rank=1
prefix-density=21.97
prefix-fanout=29.9
sequence=AGATCGGAAGAGCACACGTCTGAACTCCAGTCACGATCAGATCTCGTATGCCGTCTTC


criterion=fanout-score
sequence-density=21.28
sequence-density-rank=1
fanout-score=30.92
fanout-score-rank=1
prefix-density=21.97
prefix-fanout=29.9
sequence=AGATCGGAAGAGCACACGTCTGAACTCCAGTCACGATCAGATCTCGTATGCCGTCTTC
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 AGATCGGAAGAGCACACGTCTGAACTCCAGTCACGATCAGATCTCGTATGCCGTCTTC -o ERR4018146 -
Input file:	STDIN
trimmed:	ERR4018146-trimmed.fastq

Parameters used:
-- 3' end adapter sequence (-x):	AGATCGGAAGAGCACACGTCTGAACTCCAGTCACGATCAGATCTCGTATGCCGTCTTC
-- 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:04:31 2025 >> started

Sun Dec 28 15:04:39 2025 >> done (8.274s)
13727221 reads processed; of these:
    5953 ( 0.04%) short reads filtered out after trimming by size control
    3319 ( 0.02%) empty reads filtered out after trimming by size control
13717949 (99.93%) reads available; of these:
 4048949 (29.52%) trimmed reads available after processing
 9669000 (70.48%) untrimmed reads available after processing

Length distribution of reads after trimming:
length	count	percentage
 18	   95924	  0.70%
 19	   86162	  0.63%
 20	  151699	  1.11%
 21	  290710	  2.12%
 22	  738478	  5.38%
 23	  445994	  3.25%
 24	  289008	  2.11%
 25	   90111	  0.66%
 26	   63626	  0.46%
 27	   64361	  0.47%
 28	   85984	  0.63%
 29	   54492	  0.40%
 30	   98641	  0.72%
 31	  156517	  1.14%
 32	  120790	  0.88%
 33	   59285	  0.43%
 34	   53417	  0.39%
 35	   42447	  0.31%
 36	   56710	  0.41%
 37	   54084	  0.39%
 38	   32573	  0.24%
 39	   36017	  0.26%
 40	   22382	  0.16%
 41	   21934	  0.16%
 42	   41831	  0.30%
 43	   74599	  0.54%
 44	  114636	  0.84%
 45	  496296	  3.62%
 46	   59477	  0.43%
 47	   31525	  0.23%
 48	  196683	  1.43%
 49	      13	  0.00%
 50	      36	  0.00%
 51	 9491507	 69.19%


criterion=sequence-density
sequence-density=4.10
sequence-density-rank=1
fanout-score=1.94
fanout-score-rank=11
prefix-density=0.02
prefix-fanout=1.9
sequence=GTGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAACTCTCTGCTGATC


criterion=fanout-score
sequence-density=0.48
sequence-density-rank=14
fanout-score=39.66
fanout-score-rank=1
prefix-density=6.24
prefix-fanout=3.1
sequence=GTGATGACTTTCTTAGGAC
                                 Started job on |	Dec 28 15:04:48
                             Started mapping on |	Dec 28 15:04:48
                                    Finished on |	Dec 28 15:05:15
       Mapping speed, Million of reads per hour |	2012.09

                          Number of input reads |	15090671
                      Average input read length |	45
                                    UNIQUE READS:
                   Uniquely mapped reads number |	10377456
                        Uniquely mapped reads % |	68.77%
                          Average mapped length |	46.48
                       Number of splices: Total |	10581
            Number of splices: Annotated (sjdb) |	2507
                       Number of splices: GT/AG |	9888
                       Number of splices: GC/AG |	176
                       Number of splices: AT/AC |	3
               Number of splices: Non-canonical |	514
                      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.23
                             MULTI-MAPPING READS:
        Number of reads mapped to multiple loci |	2858011
             % of reads mapped to multiple loci |	18.94%
        Number of reads mapped to too many loci |	352487
             % of reads mapped to too many loci |	2.34%
                                  UNMAPPED READS:
       % of reads unmapped: too many mismatches |	0.00%
                 % of reads unmapped: too short |	9.87%
                     % of reads unmapped: other |	0.09%
                                  CHIMERIC READS:
                       Number of chimeric reads |	0
                            % of chimeric reads |	0.00%
N_unmapped	1855204	1855204	1855204
N_multimapping	2858011	2858011	2858011
N_noFeature	304052	422281	9927211
N_ambiguous	822056	489991	50
UnstrandedReadsAssigned:9251348 PositiveStrandReadsAssigned:9465184 NegativeStrandReadsAssigned:450195
Dataset is classified positive stranded
MeadianReadLen=51 20thPercentileLength=32 echo kmer=27
ERR4018146 Starting Kallisto single end mapping to ensembl reference transcriptome. kmer=27

[quant] fragment length distribution is truncated gaussian with mean = 100, sd = 20
[index] k-mer length: 27
[index] number of targets: 31,715
[index] number of k-mers: 49,826,120
[index] number of equivalence classes: 77,844
[quant] running in single-end mode
[quant] will process file 1: ERR4018146-trimmed.fastq
[quant] finding pseudoalignments for the reads ... done
[quant] processed 15,090,671 reads, 106,218 reads pseudoaligned
[   em] quantifying the abundances ... done
[   em] the Expectation-Maximization algorithm ran for 666 rounds

  31716 ERR4018146.ke.tsv
  32883 ERR4018146.se.tsv
  64599 total
==> ERR4018146.ke.tsv <==
target_id	length	eff_length	est_counts	tpm
ENSRNOT00000047550.4	955	856	91	693.422
ENSRNOT00000040993.4	1039	940	167	1158.83
ENSRNOT00000050156.3	1545	1446	807	3640.29
ENSRNOT00000043693.3	684	585	295	3289.25
ENSRNOT00000046201.3	204	105	134	8324.26
ENSRNOT00000046108.3	681	582	95	1064.71
ENSRNOT00000049683.3	784	685	93	885.569
ENSRNOT00000041241.3	348	249	35	916.851
ENSRNOT00000044582.3	297	198	107	3524.92

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