Starting /dee2/code/volunteer_pipeline.sh ERR4018145
    current disk space = 386633856
    free memory = 194500456 
ERR4018145 SRAfilesize
f9f579878408c12aec25d4f2010df770  ERR4018145.sra
ERR4018145.sra file validated
ERR4018145 is single end
ERR4018145 is conventional basespace
ERR4018145 read1 length is 51 nt
##FastQC	0.11.5
>>Basic Statistics	pass
#Measure	Value
Filename	ERR4018145_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	47
>>END_MODULE
>>Per base sequence quality	pass
#Base	Mean	Median	Lower Quartile	Upper Quartile	10th Percentile	90th Percentile
1	30.91725	35.0	35.0	35.0	2.0	35.0
2	33.69525	35.0	35.0	35.0	27.0	35.0
3	34.005	35.0	35.0	35.0	31.0	35.0
4	34.2955	35.0	35.0	35.0	33.0	35.0
5	33.546	35.0	35.0	35.0	31.0	35.0
6	38.21825	39.0	38.0	40.0	35.0	40.0
7	38.883	40.0	39.0	40.0	37.0	40.0
8	39.00025	40.0	39.0	40.0	38.0	40.0
9	38.86575	40.0	39.0	40.0	37.0	40.0
10	39.10025	40.0	39.0	40.0	38.0	40.0
11	39.09225	40.0	39.0	40.0	38.0	40.0
12	39.09975	40.0	39.0	40.0	38.0	40.0
13	38.93	40.0	39.0	40.0	38.0	40.0
14	38.8835	40.0	39.0	40.0	37.0	40.0
15	38.97025	40.0	39.0	40.0	37.0	40.0
16	39.0285	40.0	39.0	40.0	38.0	40.0
17	39.119	40.0	39.0	40.0	38.0	40.0
18	39.06175	40.0	39.0	40.0	38.0	40.0
19	39.112	40.0	39.0	40.0	38.0	40.0
20	39.14725	40.0	39.0	40.0	38.0	40.0
21	39.175	40.0	39.0	40.0	38.0	40.0
22	39.07675	40.0	39.0	40.0	38.0	40.0
23	39.003	40.0	39.0	40.0	38.0	40.0
24	39.07675	40.0	39.0	40.0	38.0	40.0
25	39.122	40.0	39.0	40.0	38.0	40.0
26	39.0915	40.0	39.0	40.0	38.0	40.0
27	39.12475	40.0	39.0	40.0	38.0	40.0
28	39.2	40.0	39.0	40.0	38.0	40.0
29	39.194	40.0	39.0	40.0	38.0	40.0
30	39.1365	40.0	39.0	40.0	38.0	40.0
31	39.0525	40.0	39.0	40.0	38.0	40.0
32	39.0805	40.0	39.0	40.0	38.0	40.0
33	39.08675	40.0	39.0	40.0	38.0	40.0
34	38.993	40.0	39.0	40.0	38.0	40.0
35	39.03125	40.0	39.0	40.0	38.0	40.0
36	38.9835	40.0	39.0	40.0	38.0	40.0
37	38.9975	40.0	39.0	40.0	38.0	40.0
38	38.97	40.0	39.0	40.0	38.0	40.0
39	38.955	40.0	39.0	40.0	38.0	40.0
40	38.94175	40.0	39.0	40.0	38.0	40.0
41	38.96125	40.0	39.0	40.0	38.0	40.0
42	38.89725	40.0	39.0	40.0	38.0	40.0
43	38.955	40.0	39.0	40.0	38.0	40.0
44	38.96175	40.0	39.0	40.0	38.0	40.0
45	38.945	40.0	39.0	40.0	38.0	40.0
46	38.9485	40.0	39.0	40.0	38.0	40.0
47	38.93825	40.0	39.0	40.0	38.0	40.0
48	38.9835	40.0	39.0	40.0	38.0	40.0
49	38.751	40.0	39.0	40.0	37.0	40.0
50	38.89625	40.0	39.0	40.0	37.0	40.0
51	38.39975	39.0	39.0	40.0	36.0	40.0
>>END_MODULE
>>Per sequence quality scores	pass
#Quality	Count
17	1.0
18	0.0
19	0.0
20	0.0
21	0.0
22	0.0
23	1.0
24	0.0
25	4.0
26	1.0
27	2.0
28	3.0
29	8.0
30	13.0
31	26.0
32	32.0
33	46.0
34	55.0
35	88.0
36	138.0
37	331.0
38	1183.0
39	2068.0
>>END_MODULE
>>Per base sequence content	fail
#Base	G	A	T	C
1	24.512534818941504	16.601671309192202	37.68802228412256	21.197771587743734
2	40.975	14.799999999999999	33.875	10.35
3	23.799999999999997	34.55	12.3	29.349999999999998
4	17.375	39.85	18.275	24.5
5	19.475	42.075	15.15	23.3
6	23.7	19.0	47.225	10.075000000000001
7	54.6	7.85	28.625	8.924999999999999
8	28.749999999999996	48.075	12.35	10.825
9	11.95	36.6	41.55	9.9
10	50.525	14.575	26.625	8.275
11	43.2	40.1	10.375	6.325
12	29.45	28.575	17.45	24.525
13	15.65	32.9	21.65	29.799999999999997
14	10.4	38.6	34.4	16.6
15	17.675	21.725	36.625	23.974999999999998
16	13.900000000000002	31.6	34.699999999999996	19.8
17	21.099999999999998	24.9	18.375	35.625
18	17.599999999999998	32.175	29.349999999999998	20.875
19	14.7	33.575	32.725	19.0
20	15.625	23.45	35.475	25.45
21	23.925	14.325	52.125	9.625
22	31.8	22.925	29.549999999999997	15.725
23	27.375	29.2	12.225	31.2
24	32.025	27.725	19.900000000000002	20.349999999999998
25	32.5	25.275	13.375	28.849999999999998
26	25.874999999999996	21.575	22.775000000000002	29.775000000000002
27	31.324999999999996	22.25	12.85	33.575
28	18.975	18.675	29.525000000000002	32.824999999999996
29	29.725	18.5	29.325000000000003	22.45
30	34.175	14.325	26.85	24.65
31	25.05	27.975	25.45	21.525
32	29.599999999999998	24.5	18.175	27.725
33	13.925	29.849999999999998	20.599999999999998	35.625
34	13.775	17.549999999999997	31.324999999999996	37.35
35	26.3	16.85	24.0	32.85
36	24.7	39.2	10.125	25.974999999999998
37	17.825	47.375	15.675	19.125
38	20.65	34.325	17.325	27.700000000000003
39	23.45	23.799999999999997	28.825	23.925
40	31.424999999999997	14.524999999999999	21.875	32.175
41	32.625	26.75	21.9	18.725
42	19.225	40.075	29.325000000000003	11.375
43	20.75	27.925	24.9	26.424999999999997
44	26.825	35.025	22.025	16.125
45	22.125	41.325	17.349999999999998	19.2
46	24.175	32.0	20.9	22.925
47	12.675	31.45	36.8	19.075
48	23.0	37.425000000000004	20.724999999999998	18.85
49	11.125	35.699999999999996	33.550000000000004	19.625
50	21.85	20.674999999999997	20.5	36.975
51	27.650000000000002	20.225	23.150000000000002	28.975
>>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.5
16	1.0
17	0.5
18	0.0
19	0.0
20	0.0
21	0.0
22	0.0
23	0.5
24	1.0
25	2.5
26	3.0
27	2.0
28	17.5
29	33.0
30	39.0
31	45.0
32	33.5
33	22.0
34	55.5
35	89.0
36	308.5
37	528.0
38	479.0
39	430.0
40	409.0
41	388.0
42	256.0
43	124.0
44	130.5
45	137.0
46	252.0
47	367.0
48	306.5
49	246.0
50	272.0
51	298.0
52	240.5
53	183.0
54	321.5
55	460.0
56	376.0
57	292.0
58	245.0
59	198.0
60	131.5
61	65.0
62	44.5
63	24.0
64	16.0
65	8.0
66	11.5
67	15.0
68	16.5
69	18.0
70	19.5
71	21.0
72	10.5
73	0.0
74	0.5
75	0.5
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	10.25
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	35.375
#Duplication Level	Percentage of deduplicated	Percentage of total
1	81.5547703180212	28.849999999999998
2	6.713780918727916	4.75
3	2.685512367491166	2.85
4	1.9081272084805656	2.7
5	1.342756183745583	2.375
6	0.9893992932862191	2.1
7	0.49469964664310956	1.225
8	0.4240282685512367	1.2
9	0.49469964664310956	1.575
>10	2.685512367491166	18.55
>50	0.35335689045936397	9.5
>100	0.35335689045936397	24.325
>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
AGATCGGAAGAGCACACGTCTGAACTCCAGTCACCCGTCCATCTCGTATGC	235	5.875	TruSeq Adapter, Index 16 (98% over 50bp)
TGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAACT	232	5.800000000000001	No Hit
GTAAATGATGAGATTTCCATTGGTCCGTGTTTCTGAACTCCATGATTTTTG	201	5.025	No Hit
TGCACTGATGAAACTAGAAATAGGAAGTGCCGTCCGAAACGACAACTGACG	158	3.95	No Hit
GTGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAAC	147	3.675	No Hit
GGAAGAGCACACGTCTGAACTCCAGTCACCCGTCCATCTCGTATGCCGTCT	98	2.45	TruSeq Adapter, Index 16 (98% over 51bp)
TTGCTGTGATGACTATCTTAGGACACCTTTGGATTAACCGTGAAATAAAAC	90	2.25	No Hit
CTCGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGTAGATCG	87	2.175	No Hit
TACTGTGATGATAACATAGTTCAGCAGACGTGACTCTGGTGAACAATCATG	53	1.325	No Hit
TTGCAATGATGTGAATCTTTTGCTGAAGTGACCTTGAAGTGTAATTCCTGA	52	1.3	No Hit
CGCGACCTCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACCCGTCCAT	50	1.25	TruSeq Adapter, Index 16 (97% over 41bp)
ATACATGATGATCTCACTCCAACTTGAACTCTCTCACTGATTACTTGATGA	37	0.9249999999999999	No Hit
CACAATGATGACTTATGTTTTTGCCGTTTACCCAGCTGAGGGTTTCTTTGA	34	0.8500000000000001	No Hit
TAGCTTATCAGACTGATGTTGACAGATCGGAAGAGCACACGTCTGAACTCC	33	0.8250000000000001	Illumina Multiplexing PCR Primer 2.01 (100% over 28bp)
TGGAATGATGACATTCTCCGGAATCGCTGTACTGTCTTGATGAAAGTACTT	33	0.8250000000000001	No Hit
GGTCTGTGATGAAACCCATTATCAGTGGACATCTATGGATGAGAAATGCGG	31	0.775	No Hit
NGATCGGAAGAGCACACGTCTGAACTCCAGTCACCCGTCCATCTCGTATGC	28	0.7000000000000001	TruSeq Adapter, Index 16 (98% over 50bp)
CACCCTGATGAACTGAATACCGCCCCAGTCTGATAGATGTGGAGAAAGGTA	28	0.7000000000000001	No Hit
TGCTGTGATGAGGACACCAGTAGGAAGTGCCGTCAGAGTCGATAACTGACG	27	0.675	No Hit
CTGCGATGATGTCATTTCTTAGGACACCTTTGGATTAATCATGAAAATAGC	25	0.625	No Hit
TTGAGTGATGACTTGTATTGTCGGATACCCCTTCACCCCGTTCATGGGTGA	25	0.625	No Hit
NGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAACT	24	0.6	No Hit
GCTAATGATGATACCGCACTAACTGAGCAGTTAGTAGTTGGTCCTTTGATT	23	0.575	No Hit
NTAAATGATGAGATTTCCATTGGTCCGTGTTTCTGAACTCCATGATTTTTG	22	0.5499999999999999	No Hit
CTTAATGATGATTGTCTATTTGATTGCTTGATGCAATGTGAAAACCACATT	21	0.525	No Hit
NGCACTGATGAAACTAGAAATAGGAAGTGCCGTCCGAAACGACAACTGACG	19	0.475	No Hit
GGTCAATGATGTGTTGGCATGTATTATCTGAATTTGCTGATGTGTCCTAAC	19	0.475	No Hit
CTTCTGTGATGACTTTACCCAAATGACTTTCGTTCTTCTGAGTTTGCTGAA	17	0.42500000000000004	No Hit
TCGATGTGGTGCTGCCGCGTTCTCTTCGGGCCGGGGCCTAAGCCGCGCCAG	16	0.4	No Hit
TGAGGTAGTAGTTTGTGCTGTTAGATCGGAAGAGCACACGTCTGAACTCCA	16	0.4	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
TCGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGTAGATCGG	16	0.4	No Hit
CACCCTGATGAACTGAATACCGCCCCAGTCTGATAGCTGTGGAGAAAGGTA	15	0.375	No Hit
ATCGATGTGGTGCTGCCGCGTTCTCTTCGGGCCGGGGCCTAAGCCGCGCCA	14	0.35000000000000003	No Hit
NTGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAAC	14	0.35000000000000003	No Hit
GTTCTGTGATGAGACTCTGGAATGTGCGCTGGGCACAGCGCCCGAGTGAGA	13	0.325	No Hit
ACCCCTGTGATGTTCAGCACTGGGCTCTGATCACCCTGAGGACATGGAGCT	13	0.325	No Hit
TTCAAGTAATCCAGGATAGGCTAGATCGGAAGAGCACACGTCTGAACTCCA	12	0.3	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
AACAGTGATGACCACAACATGCCATCTGAGTTGCTTTGCTGAAATCCAGAG	11	0.27499999999999997	No Hit
NTGCTGTGATGACTATCTTAGGACACCTTTGGATTAACCGTGAAATAAAAC	11	0.27499999999999997	No Hit
GACTCTTAGCGGTGGATCACTCGGCTCGTGCGAGATCGGAAGAGCACACGT	11	0.27499999999999997	No Hit
NTCGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGTAGATCG	11	0.27499999999999997	No Hit
CACAAATGATGAATCCTTTGACGGGCGGGCAGAATCTCTGTGCTGATTGTC	11	0.27499999999999997	No Hit
AAATATGATGAGCATACAAAATGTTGATTACAGTGAAAAACAGTAAGATCC	11	0.27499999999999997	No Hit
CGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGTAGATCGGA	11	0.27499999999999997	No Hit
GCGGATGATGACACCTGGGTATGCTGCACATTCCCGACTGCGTCGTGGGGA	10	0.25	No Hit
CTCGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGAAGATCG	10	0.25	No Hit
AGTAAATGATGAGATTTCCATTGGTCCGTGTTTCTGAACTCCATGATTTTT	10	0.25	No Hit
CGACTCTTAGCGGTGGATCACTCGGCTCGTGCGTCGATGAAGAACGCAGCT	10	0.25	No Hit
GTGCAATGATGTGTGAATTTCTTCACCTGAAATCACCGTGAAGAGTAATTT	9	0.22499999999999998	No Hit
TCAAGTGATGGAGACATGACCCTGGGAGGCCCATGACCGAACTGTCTGTAT	9	0.22499999999999998	No Hit
TTCTGTGATGACTTTACCCAAATGACTTTCGTTCTTCTGAGTTTGCTGAAG	9	0.22499999999999998	No Hit
CTCCCATGATGTATAGCACTGGGCTCTGACCACCATGGTGACGTGTTGTGC	9	0.22499999999999998	No Hit
TTTGAGTGATGACTTGTATTGTCGGATACCCCTTCACCCCGTTCATGGGTG	9	0.22499999999999998	No Hit
CGACTCTTAGCGGTGGATCACTCGGCTCGTGAGATCGGAAGAGCACACGTC	9	0.22499999999999998	No Hit
CACAAGTGATGACTAACAAAGGGACTTAATACTGAAACATCTGGGATTGAA	9	0.22499999999999998	No Hit
GACTCTTAGCGGTGGATCACTCGGCTCGTGCGTCGATGAAGAACGCAGCTA	8	0.2	No Hit
CGCGACCTCAGATCAGACGTGGCGACCCGCTGAATTTAAGCATATTAGTCA	8	0.2	No Hit
GCACAATGATGACTTATGTTTTTGCCGTTTACCCAGCTGAGGGTTTCTTTG	8	0.2	No Hit
CTGCAGTGATGACTTTCTTAGGACACCTTTGGATTTACCGTGAAAATAAGT	8	0.2	No Hit
CTACAATGATGGTAATATTTTTCGTCAACAGAGTTCACCTAGTGAGTGTTG	8	0.2	No Hit
CGCGACCTTAGATCGGAAGAGCACACGTCTGAACTCCAGTCACCCGTCCAT	8	0.2	TruSeq Adapter, Index 16 (97% over 41bp)
CGCGACCTCAGATCAGACGTGGCGACCCGCTGAATTAGATCGGAAGAGCAC	7	0.17500000000000002	No Hit
GTGGAATGATGACATTCTCCGGAATCGCTGTACTGTCTTGATGAAAGTACT	7	0.17500000000000002	No Hit
CCGCGACCTCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACCCGTCCA	7	0.17500000000000002	TruSeq Adapter, Index 16 (97% over 39bp)
ATCTATGATGATCCTATCCCGAATCTGAATTCCTGTTGAAAAAATATAATT	7	0.17500000000000002	No Hit
TGAAGTGATGATTCACAATCATGTCTCTTCTCTGAAGAATTCTTGGAGATT	7	0.17500000000000002	No Hit
TAGCTTATCAGACTGATGTTGACTAGATCGGAAGAGCACACGTCTGAACTC	7	0.17500000000000002	Illumina Multiplexing PCR Primer 2.01 (100% over 27bp)
CGTTGGATCGATGTGGTGCTGCCGCGTTCTCTTCGGGCCGGGGCCTAAGCC	7	0.17500000000000002	No Hit
CCGATACGACTCTTAGCGGTGGATCACTCGGCTCGTGCGTCGATGAAGAAC	6	0.15	No Hit
NGCGACCTCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACCCGTCCAT	6	0.15	TruSeq Adapter, Index 16 (97% over 41bp)
GACTCTTAGCGGTGGATCACTCGGCTCGTGAGATCGGAAGAGCACACGTCT	6	0.15	Illumina Multiplexing PCR Primer 2.01 (100% over 21bp)
CGCGACCTCAGATCAGACAGATCGGAAGAGCACACGTCTGAACTCCAGTCA	6	0.15	Illumina Multiplexing PCR Primer 2.01 (100% over 33bp)
TTCTGTGATGAGACTCTGGAATGTGCGCTGGGCACAGCGCCCGAGTGAGAC	6	0.15	No Hit
NTTAATGATGATTGTCTATTTGATTGCTTGATGCAATGTGAAAACCACATT	6	0.15	No Hit
ACTTAATGATGATTGTCTATTTGATTGCTTGATGCAATGTGAAAACCACAT	6	0.15	No Hit
NACTGTGATGATAACATAGTTCAGCAGACGTGACTCTGGTGAACAATCATG	6	0.15	No Hit
TAGCTTATCAGACTGATGTTGACAAGATCGGAAGAGCACACGTCTGAACTC	6	0.15	Illumina Multiplexing PCR Primer 2.01 (100% over 27bp)
NGAAGAGCACACGTCTGAACTCCAGTCACCCGTCCATCTCGTATGCCGTCT	6	0.15	TruSeq Adapter, Index 16 (97% over 35bp)
CTTCGGTGATGACACGATGACGAGTCAGAAAGGCCACGTCTTGCTCTTGGG	6	0.15	No Hit
CCTAGTGATGACAAGGAGACATTGTCAACCAATCCCCACAAGGGAATGAGG	6	0.15	No Hit
CGCGACCTCAGATCAGACGTAGATCGGAAGAGCACACGTCTGAACTCCAGT	6	0.15	Illumina Multiplexing PCR Primer 2.01 (100% over 31bp)
TGCAATGATGTGTGAATTTCTTCACCTGAAATCACCGTGAAGAGTAATTTA	6	0.15	No Hit
AGTCTGTGATGAGCAATAATCACCATCTTTCGTTTGAGTCTCACGACTGTG	5	0.125	No Hit
CGCGACCTCAGATCAGAAGATCGGAAGAGCACACGTCTGAACTCCAGTCAC	5	0.125	Illumina Multiplexing PCR Primer 2.01 (100% over 34bp)
CGCGACCTCAGATCAGACGTGGCGACCCGCTGAAAGATCGGAAGAGCACAC	5	0.125	No Hit
ACAAGTGATGACTAACAAAGGGACTTAATACTGAAACATCTGGGATTGAAA	5	0.125	No Hit
NTGAGTGATGACTTGTATTGTCGGATACCCCTTCACCCCGTTCATGGGTGA	5	0.125	No Hit
AATACATGATGATCTCACTCCAACTTGAACTCTCTCACTGATTACTTGATG	5	0.125	No Hit
AACCCGTAGATCCGATCTTGTAGATCGGAAGAGCACACGTCTGAACTCCAG	5	0.125	Illumina Multiplexing PCR Primer 2.01 (100% over 30bp)
CTCGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGCAGATCG	5	0.125	No Hit
TGAGGTAGTAGGTTGTATAGTTAGATCGGAAGAGCACACGTCTGAACTCCA	5	0.125	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
ATTGATCATCGACACTTCGAACGCACTTGCGGCCCCGGGTTCCTCCCGGGG	5	0.125	No Hit
TGCTGTGATGACTATCTTAGGACACCTTTGGATTAACCGTGAAATAAAACA	5	0.125	No Hit
ACAAATGATGAATCCTTTGACGGGCGGGCAGAATCTCTGTGCTGATTGTCA	5	0.125	No Hit
ATACTGTGATGATAACATAGTTCAGCAGACGTGACTCTGGTGAACAATCAT	5	0.125	No Hit
NTGCGATGATGTCATTTCTTAGGACACCTTTGGATTAATCATGAAAATAGC	5	0.125	No Hit
TGCTAATGATGATACCGCACTAACTGAGCAGTTAGTAGTTGGTCCTTTGAT	5	0.125	No Hit
TCGGATCCGTCTGAGCTTGGCTAGATCGGAAGAGCACACGTCTGAACTCCA	5	0.125	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
CGCGACCTCAGATCAGACGTGGCGACCCGCTGAATTTAGATCGGAAGAGCA	5	0.125	No Hit
NGGAATGATGACATTCTCCGGAATCGCTGTACTGTCTTGATGAAAGTACTT	5	0.125	No Hit
TGGAGATGATGAACTGACTCAACCTGACCTTCAGAAATGGAGGTGATAAAA	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	6.1	0.0	0.0	0.0	0.0
2	6.15	0.0	0.0	0.0	0.0
3	6.175	0.0	0.0	0.0	0.0
4	6.2	0.0	0.0	0.0	0.0
5	6.275	0.0	0.0	0.0	0.0
6	6.375	0.0	0.0	0.0	0.0
7	6.425	0.0	0.0	0.0	0.0
8	6.525	0.0	0.0	0.0	0.0
9	6.8	0.0	0.0	0.0	0.0
10	8.95	0.0	0.0	0.0	0.0
11	9.525	0.0	0.0	0.0	0.0
12	10.05	0.0	0.0	0.0	0.0
13	10.275	0.0	0.0	0.0	0.0
14	10.875	0.0	0.0	0.0	0.0
15	11.45	0.0	0.0	0.0	0.0
16	11.9	0.0	0.0	0.0	0.0
17	12.325	0.0	0.0	0.0	0.0
18	12.7	0.0	0.0	0.0	0.0
19	13.4	0.0	0.0	0.0	0.0
20	14.0	0.0	0.0	0.0	0.0
21	14.725	0.0	0.0	0.0	0.0
22	16.475	0.0	0.0	0.0	0.0
23	20.3	0.0	0.0	0.0	0.0
24	23.025	0.0	0.0	0.0	0.0
25	24.675	0.0	0.0	0.0	0.0
26	25.15	0.0	0.0	0.0	0.0
27	25.425	0.0	0.0	0.0	0.0
28	25.7	0.0	0.0	0.0	0.0
29	26.05	0.0	0.0	0.0	0.0
30	26.35	0.0	0.0	0.0	0.0
31	27.0	0.0	0.0	0.0	0.0
32	27.8	0.0	0.0	0.0	0.0
33	28.475	0.0	0.0	0.0	0.0
34	28.975	0.0	0.0	0.0	0.0
35	29.175	0.0	0.0	0.0	0.0
36	29.4	0.0	0.0	0.0	0.0
37	29.775	0.0	0.0	0.0	0.0
38	30.25	0.0	0.0	0.0	0.0
39	30.45	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.007837298	49.88889	1
GTAAATG	25	1.721153E-5	49.888885	1
TGATGAG	45	2.237357E-10	44.9	6
CATGATT	40	4.443791E-9	44.9	41
TGGTCCG	30	1.6973481E-6	44.9	21
GTGTTTC	40	4.443791E-9	44.9	27
CCGTGTT	35	8.720781E-8	44.9	25
CATTGGT	35	8.720781E-8	44.9	18
TAAATGA	30	1.6973481E-6	44.9	2
TAGGAAG	20	6.2916183E-4	44.9	21
ATGATTT	40	4.443791E-9	44.9	42
ACTGACG	20	6.2916183E-4	44.9	45
TTTCTGA	40	4.443791E-9	44.9	30
AACTGAC	20	6.2916183E-4	44.9	44
TTTCCAT	40	4.443791E-9	44.9	14
GTCCGTG	30	1.6973481E-6	44.9	23
GGTCCGT	30	1.6973481E-6	44.9	22
TGATTTT	40	4.443791E-9	44.9	43
GAAGTGC	20	6.2916183E-4	44.9	24
ATTGGTC	30	1.6973481E-6	44.9	19
>>END_MODULE
Rejected 2922350 READS because READLEN < 1
Read 2922350 spots for ERR4018145.sra
Written 2922350 spots for ERR4018145.sra
Rejected 2922347 READS because READLEN < 1
Read 2922347 spots for ERR4018145.sra
Written 2922347 spots for ERR4018145.sra
Rejected 2922347 READS because READLEN < 1
Read 2922347 spots for ERR4018145.sra
Written 2922347 spots for ERR4018145.sra
Rejected 2922347 READS because READLEN < 1
Read 2922347 spots for ERR4018145.sra
Written 2922347 spots for ERR4018145.sra
Rejected 2922347 READS because READLEN < 1
Read 2922347 spots for ERR4018145.sra
Written 2922347 spots for ERR4018145.sra
Rejected 2922347 READS because READLEN < 1
Read 2922347 spots for ERR4018145.sra
Written 2922347 spots for ERR4018145.sra
Rejected 2922347 READS because READLEN < 1
Read 2922347 spots for ERR4018145.sra
Written 2922347 spots for ERR4018145.sra
Rejected 2922347 READS because READLEN < 1
Read 2922347 spots for ERR4018145.sra
Written 2922347 spots for ERR4018145.sra
SRR ids: ['ERR4018145.sra']
extra args: ['--split-files', '--defline-qual', '+']
tempdir: /tmp/pfd_ql3cowot
ERR4018145.sra spots: 23378779
blocks: [[1, 2922347], [2922348, 5844694], [5844695, 8767041], [8767042, 11689388], [11689389, 14611735], [14611736, 17534082], [17534083, 20456429], [20456430, 23378779]]
ERR4018145 file size 3311604
ERR4018145 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 ERR4018145 ERR4018145_1.fastq
Input file:	ERR4018145_1.fastq
trimmed:	ERR4018145-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:00:58 2025 >> started

Sun Dec 28 15:01:05 2025 >> done (6.770s)
23378779 reads processed; of these:
 1566697 ( 6.70%) short reads filtered out after trimming by size control
 1559735 ( 6.67%) empty reads filtered out after trimming by size control
20252347 (86.63%) reads available; of these:
  249947 ( 1.23%) trimmed reads available after processing
20002400 (98.77%) untrimmed reads available after processing

Length distribution of reads after trimming:
length	count	percentage
 18	  134731	  0.67%
 19	  115102	  0.57%
 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	      43	  0.00%
 50	      71	  0.00%
 51	20002400	 98.77%
20252347 reads passed initial QC


criterion=sequence-density
sequence-density=19.13
sequence-density-rank=1
fanout-score=33.56
fanout-score-rank=1
prefix-density=19.93
prefix-fanout=32.2
sequence=AGATCGGAAGAGCACACGTCTGAACTCCAGTCACCCGTCCATCTCGTATGCCGTCTTC


criterion=fanout-score
sequence-density=19.13
sequence-density-rank=1
fanout-score=33.56
fanout-score-rank=1
prefix-density=19.93
prefix-fanout=32.2
sequence=AGATCGGAAGAGCACACGTCTGAACTCCAGTCACCCGTCCATCTCGTATGCCGTCTTC
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 AGATCGGAAGAGCACACGTCTGAACTCCAGTCACCCGTCCATCTCGTATGCCGTCTTC -o ERR4018145 -
Input file:	STDIN
trimmed:	ERR4018145-trimmed.fastq

Parameters used:
-- 3' end adapter sequence (-x):	AGATCGGAAGAGCACACGTCTGAACTCCAGTCACCCGTCCATCTCGTATGCCGTCTTC
-- 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:01:19 2025 >> started

Sun Dec 28 15:01:31 2025 >> done (11.572s)
18227112 reads processed; of these:
    7760 ( 0.04%) short reads filtered out after trimming by size control
    2993 ( 0.02%) empty reads filtered out after trimming by size control
18216359 (99.94%) reads available; of these:
 5148172 (28.26%) trimmed reads available after processing
13068187 (71.74%) untrimmed reads available after processing

Length distribution of reads after trimming:
length	count	percentage
 18	  118239	  0.65%
 19	  105471	  0.58%
 20	  165284	  0.91%
 21	  345840	  1.90%
 22	  797908	  4.38%
 23	  555925	  3.05%
 24	  340630	  1.87%
 25	  114818	  0.63%
 26	   82075	  0.45%
 27	   80408	  0.44%
 28	   99807	  0.55%
 29	   70284	  0.39%
 30	  126043	  0.69%
 31	  197072	  1.08%
 32	  154459	  0.85%
 33	   78956	  0.43%
 34	   76504	  0.42%
 35	   56951	  0.31%
 36	   78607	  0.43%
 37	   74447	  0.41%
 38	   50364	  0.28%
 39	   55832	  0.31%
 40	   36661	  0.20%
 41	   36589	  0.20%
 42	   60279	  0.33%
 43	  112672	  0.62%
 44	  172228	  0.95%
 45	  702358	  3.86%
 46	  116319	  0.64%
 47	   48336	  0.27%
 48	  256601	  1.41%
 49	      28	  0.00%
 50	      47	  0.00%
 51	12848317	 70.53%


criterion=sequence-density
sequence-density=5.20
sequence-density-rank=1
fanout-score=1.99
fanout-score-rank=9
prefix-density=0.02
prefix-fanout=2.0
sequence=GTGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAACTCTCTGCTGATC


criterion=fanout-score
sequence-density=0.97
sequence-density-rank=5
fanout-score=7.78
fanout-score-rank=1
prefix-density=7.49
prefix-fanout=1.0
sequence=GAAGTGCCGTCAGAGTCGATAACTGACGATA
                                 Started job on |	Dec 28 15:01:40
                             Started mapping on |	Dec 28 15:01:40
                                    Finished on |	Dec 28 15:02:18
       Mapping speed, Million of reads per hour |	1917.62

                          Number of input reads |	20241594
                      Average input read length |	45
                                    UNIQUE READS:
                   Uniquely mapped reads number |	14383843
                        Uniquely mapped reads % |	71.06%
                          Average mapped length |	47.11
                       Number of splices: Total |	14741
            Number of splices: Annotated (sjdb) |	3404
                       Number of splices: GT/AG |	13750
                       Number of splices: GC/AG |	234
                       Number of splices: AT/AC |	5
               Number of splices: Non-canonical |	752
                      Mismatch rate per base, % |	0.39%
                         Deletion rate per base |	0.00%
                        Deletion average length |	1.07
                        Insertion rate per base |	0.00%
                       Insertion average length |	1.25
                             MULTI-MAPPING READS:
        Number of reads mapped to multiple loci |	3873856
             % of reads mapped to multiple loci |	19.14%
        Number of reads mapped to too many loci |	412588
             % of reads mapped to too many loci |	2.04%
                                  UNMAPPED READS:
       % of reads unmapped: too many mismatches |	0.00%
                 % of reads unmapped: too short |	7.66%
                     % of reads unmapped: other |	0.10%
                                  CHIMERIC READS:
                       Number of chimeric reads |	0
                            % of chimeric reads |	0.00%
N_unmapped	1983895	1983895	1983895
N_multimapping	3873856	3873856	3873856
N_noFeature	437715	575028	13826418
N_ambiguous	1152261	732091	64
UnstrandedReadsAssigned:12793867 PositiveStrandReadsAssigned:13076724 NegativeStrandReadsAssigned:557361
Dataset is classified positive stranded
MeadianReadLen=51 20thPercentileLength=37 echo kmer=33
ERR4018145 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: ERR4018145-trimmed.fastq
[quant] finding pseudoalignments for the reads ... done
[quant] processed 20,241,594 reads, 130,729 reads pseudoaligned
[   em] quantifying the abundances ... done
[   em] the Expectation-Maximization algorithm ran for 635 rounds

  31716 ERR4018145.ke.tsv
  32883 ERR4018145.se.tsv
  64599 total
==> ERR4018145.ke.tsv <==
target_id	length	eff_length	est_counts	tpm
ENSRNOT00000047550.4	955	856	102	645.584
ENSRNOT00000040993.4	1039	940	132	760.803
ENSRNOT00000050156.3	1545	1446	964	3611.89
ENSRNOT00000043693.3	684	585	265	2454.24
ENSRNOT00000046201.3	204	105	146	7533.38
ENSRNOT00000046108.3	681	582	69	642.321
ENSRNOT00000049683.3	784	685	104	822.563
ENSRNOT00000041241.3	348	249	69	1501.33
ENSRNOT00000044582.3	297	198	85	2325.84

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