Starting /dee2/code/volunteer_pipeline.sh ERR4018148
    current disk space = 365850364
    free memory = 199151228 
ERR4018148 SRAfilesize
c652b2703cb47be966162fd0b5cd7fc7  ERR4018148.sra
ERR4018148.sra file validated
ERR4018148 is single end
ERR4018148 is conventional basespace
ERR4018148 read1 length is 51 nt
##FastQC	0.11.5
>>Basic Statistics	pass
#Measure	Value
Filename	ERR4018148_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	31.306	35.0	35.0	35.0	15.0	35.0
2	33.754	35.0	35.0	35.0	27.0	35.0
3	34.172	35.0	35.0	35.0	32.0	35.0
4	34.331	35.0	35.0	35.0	33.0	35.0
5	33.68	35.0	35.0	35.0	31.0	35.0
6	38.29775	39.0	38.0	40.0	35.0	40.0
7	38.941	40.0	39.0	40.0	38.0	40.0
8	39.1	40.0	39.0	40.0	38.0	40.0
9	38.93275	40.0	39.0	40.0	37.0	40.0
10	39.13925	40.0	39.0	40.0	38.0	40.0
11	39.14625	40.0	39.0	40.0	38.0	40.0
12	39.14925	40.0	39.0	40.0	38.0	40.0
13	39.09275	40.0	39.0	40.0	38.0	40.0
14	39.01675	40.0	39.0	40.0	38.0	40.0
15	39.0675	40.0	39.0	40.0	38.0	40.0
16	39.09575	40.0	39.0	40.0	38.0	40.0
17	39.247	40.0	39.0	40.0	38.0	40.0
18	39.154	40.0	39.0	40.0	38.0	40.0
19	39.07525	40.0	39.0	40.0	38.0	40.0
20	39.13725	40.0	39.0	40.0	38.0	40.0
21	39.226	40.0	39.0	40.0	38.0	40.0
22	39.20575	40.0	39.0	40.0	38.0	40.0
23	39.16025	40.0	39.0	40.0	38.0	40.0
24	39.139	40.0	39.0	40.0	38.0	40.0
25	39.193	40.0	39.0	40.0	38.0	40.0
26	39.12525	40.0	39.0	40.0	38.0	40.0
27	39.21575	40.0	39.0	40.0	39.0	40.0
28	39.16275	40.0	39.0	40.0	38.0	40.0
29	39.12625	40.0	39.0	40.0	38.0	40.0
30	39.131	40.0	39.0	40.0	38.0	40.0
31	39.08975	40.0	39.0	40.0	38.0	40.0
32	39.0365	40.0	39.0	40.0	38.0	40.0
33	39.0555	40.0	39.0	40.0	38.0	40.0
34	39.13575	40.0	39.0	40.0	38.0	40.0
35	39.067	40.0	39.0	40.0	38.0	40.0
36	38.988	40.0	39.0	40.0	38.0	40.0
37	38.98325	40.0	39.0	40.0	38.0	40.0
38	38.959	40.0	39.0	40.0	38.0	40.0
39	38.95475	40.0	39.0	40.0	38.0	40.0
40	38.9275	40.0	39.0	40.0	38.0	40.0
41	39.021	40.0	39.0	40.0	38.0	40.0
42	38.94725	40.0	39.0	40.0	38.0	40.0
43	38.9905	40.0	39.0	40.0	38.0	40.0
44	38.99025	40.0	39.0	40.0	38.0	40.0
45	38.98475	40.0	39.0	40.0	38.0	40.0
46	39.0125	40.0	39.0	40.0	38.0	40.0
47	39.05625	40.0	39.0	40.0	38.0	40.0
48	39.01625	40.0	39.0	40.0	38.0	40.0
49	38.7245	40.0	39.0	40.0	37.0	40.0
50	38.9075	40.0	39.0	40.0	37.0	40.0
51	38.34025	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	1.0
20	1.0
21	0.0
22	0.0
23	1.0
24	0.0
25	1.0
26	2.0
27	3.0
28	3.0
29	11.0
30	11.0
31	19.0
32	22.0
33	30.0
34	70.0
35	59.0
36	149.0
37	284.0
38	1219.0
39	2113.0
>>END_MODULE
>>Per base sequence content	fail
#Base	G	A	T	C
1	23.282758620689656	20.413793103448278	33.682758620689654	22.620689655172413
2	44.574999999999996	15.024999999999999	30.25	10.15
3	20.875	37.1	12.775	29.25
4	17.224999999999998	36.5	21.925	24.349999999999998
5	20.25	37.0	15.1	27.650000000000002
6	26.674999999999997	18.825	43.0	11.5
7	55.7	7.875	26.85	9.575
8	26.200000000000003	49.45	13.450000000000001	10.9
9	11.774999999999999	38.75	39.300000000000004	10.174999999999999
10	51.24999999999999	15.049999999999999	24.2	9.5
11	39.025	44.574999999999996	9.625	6.775
12	30.125	29.575000000000003	16.45	23.849999999999998
13	16.625	30.925000000000004	20.275000000000002	32.175
14	11.200000000000001	38.574999999999996	33.125	17.1
15	18.3	21.349999999999998	32.574999999999996	27.775
16	14.649999999999999	34.1	33.225	18.025
17	21.525	24.5	16.725	37.25
18	20.775	31.3	29.65	18.275
19	14.274999999999999	31.65	35.199999999999996	18.875
20	15.325	23.95	31.974999999999998	28.749999999999996
21	22.8	14.274999999999999	52.449999999999996	10.475
22	33.225	22.95	27.425	16.400000000000002
23	25.650000000000002	31.874999999999996	11.25	31.225
24	31.1	31.574999999999996	17.599999999999998	19.725
25	31.45	23.724999999999998	13.450000000000001	31.374999999999996
26	23.400000000000002	21.65	26.450000000000003	28.499999999999996
27	29.675	19.875	12.675	37.775
28	19.35	19.025	26.174999999999997	35.449999999999996
29	28.625	22.3	26.724999999999998	22.35
30	38.9	14.924999999999999	25.35	20.825
31	27.6	25.85	27.750000000000004	18.8
32	26.450000000000003	24.125	15.975	33.45
33	15.475	32.725	20.599999999999998	31.2
34	13.600000000000001	22.1	27.775	36.525
35	33.5	17.775	22.650000000000002	26.075
36	33.550000000000004	37.125	9.950000000000001	19.375
37	11.799999999999999	43.65	15.9	28.65
38	19.925	31.45	22.175	26.450000000000003
39	24.2	33.2	26.075	16.525000000000002
40	29.075	14.374999999999998	21.175	35.375
41	32.525	28.625	21.15	17.7
42	18.625	35.55	34.150000000000006	11.675
43	19.675	27.625	23.35	29.349999999999998
44	27.950000000000003	32.074999999999996	25.2	14.774999999999999
45	23.825	38.824999999999996	17.424999999999997	19.925
46	25.5	29.9	19.7	24.9
47	12.35	28.875	39.85	18.925
48	24.349999999999998	40.125	19.525000000000002	16.0
49	12.35	32.6	35.875	19.175
50	26.1	20.325	19.325	34.25
51	25.624999999999996	20.95	22.45	30.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.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	3.0
26	4.5
27	3.0
28	8.0
29	13.0
30	26.0
31	39.0
32	31.5
33	24.0
34	50.0
35	76.0
36	266.5
37	457.0
38	414.5
39	372.0
40	361.0
41	350.0
42	234.5
43	119.0
44	127.5
45	136.0
46	268.0
47	400.0
48	317.5
49	235.0
50	281.5
51	328.0
52	456.0
53	584.0
54	473.5
55	363.0
56	282.0
57	201.0
58	184.0
59	167.0
60	107.5
61	48.0
62	33.5
63	19.0
64	13.5
65	8.0
66	13.0
67	18.0
68	17.5
69	17.0
70	15.5
71	14.0
72	8.5
73	3.0
74	1.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	9.375
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.8
#Duplication Level	Percentage of deduplicated	Percentage of total
1	83.03072625698324	29.725
2	6.424581005586592	4.6
3	2.653631284916201	2.85
4	1.3966480446927374	2.0
5	0.9078212290502794	1.625
6	0.6983240223463687	1.5
7	0.41899441340782123	1.05
8	0.5586592178770949	1.6
9	0.5586592178770949	1.7999999999999998
>10	2.653631284916201	19.35
>50	0.27932960893854747	6.625
>100	0.41899441340782123	27.275
>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
AGATCGGAAGAGCACACGTCTGAACTCCAGTCACGGCTACATCTCGTATGC	349	8.725	TruSeq Adapter, Index 11 (100% over 50bp)
TGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAACT	191	4.775	No Hit
GTAAATGATGAGATTTCCATTGGTCCGTGTTTCTGAACTCCATGATTTTTG	156	3.9	No Hit
TGCACTGATGAAACTAGAAATAGGAAGTGCCGTCCGAAACGACAACTGACG	156	3.9	No Hit
GGAAGAGCACACGTCTGAACTCCAGTCACGGCTACATCTCGTATGCCGTCT	130	3.25	TruSeq Adapter, Index 11 (100% over 51bp)
GTGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAAC	109	2.725	No Hit
CTCGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGTAGATCG	87	2.175	No Hit
TTGCTGTGATGACTATCTTAGGACACCTTTGGATTAACCGTGAAATAAAAC	71	1.775	No Hit
CGCGACCTCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACGGCTACAT	56	1.4000000000000001	TruSeq Adapter, Index 11 (100% over 41bp)
TACTGTGATGATAACATAGTTCAGCAGACGTGACTCTGGTGAACAATCATG	51	1.275	No Hit
ATACATGATGATCTCACTCCAACTTGAACTCTCTCACTGATTACTTGATGA	44	1.0999999999999999	No Hit
TTGCAATGATGTGAATCTTTTGCTGAAGTGACCTTGAAGTGTAATTCCTGA	41	1.0250000000000001	No Hit
NGATCGGAAGAGCACACGTCTGAACTCCAGTCACGGCTACATCTCGTATGC	40	1.0	TruSeq Adapter, Index 11 (100% over 50bp)
CACAATGATGACTTATGTTTTTGCCGTTTACCCAGCTGAGGGTTTCTTTGA	37	0.9249999999999999	No Hit
GGTCTGTGATGAAACCCATTATCAGTGGACATCTATGGATGAGAAATGCGG	35	0.8750000000000001	No Hit
CACCCTGATGAACTGAATACCGCCCCAGTCTGATAGATGTGGAGAAAGGTA	35	0.8750000000000001	No Hit
TGGAATGATGACATTCTCCGGAATCGCTGTACTGTCTTGATGAAAGTACTT	35	0.8750000000000001	No Hit
TAGCTTATCAGACTGATGTTGACAGATCGGAAGAGCACACGTCTGAACTCC	32	0.8	Illumina Multiplexing PCR Primer 2.01 (100% over 28bp)
CTGCGATGATGTCATTTCTTAGGACACCTTTGGATTAATCATGAAAATAGC	30	0.75	No Hit
NGCACTGATGAAACTAGAAATAGGAAGTGCCGTCCGAAACGACAACTGACG	28	0.7000000000000001	No Hit
GCTAATGATGATACCGCACTAACTGAGCAGTTAGTAGTTGGTCCTTTGATT	25	0.625	No Hit
GTTCTGTGATGAGACTCTGGAATGTGCGCTGGGCACAGCGCCCGAGTGAGA	24	0.6	No Hit
TGCTGTGATGAGGACACCAGTAGGAAGTGCCGTCAGAGTCGATAACTGACG	23	0.575	No Hit
TCGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGTAGATCGG	20	0.5	No Hit
CTTAATGATGATTGTCTATTTGATTGCTTGATGCAATGTGAAAACCACATT	20	0.5	No Hit
TTGAGTGATGACTTGTATTGTCGGATACCCCTTCACCCCGTTCATGGGTGA	20	0.5	No Hit
NGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAACT	17	0.42500000000000004	No Hit
CTTCTGTGATGACTTTACCCAAATGACTTTCGTTCTTCTGAGTTTGCTGAA	17	0.42500000000000004	No Hit
GGTCAATGATGTGTTGGCATGTATTATCTGAATTTGCTGATGTGTCCTAAC	17	0.42500000000000004	No Hit
CGACTCTTAGCGGTGGATCACTCGGCTCGTGCGTCGATGAAGAACGCAGCT	17	0.42500000000000004	No Hit
CGACTCTTAGCGGTGGATCACTCGGCTCGTGAGATCGGAAGAGCACACGTC	16	0.4	No Hit
TGAGGTAGTAGTTTGTGCTGTTAGATCGGAAGAGCACACGTCTGAACTCCA	15	0.375	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
NTAAATGATGAGATTTCCATTGGTCCGTGTTTCTGAACTCCATGATTTTTG	14	0.35000000000000003	No Hit
TCAGTGCACTACAGAACTTTGTAGATCGGAAGAGCACACGTCTGAACTCCA	13	0.325	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
AGTAAATGATGAGATTTCCATTGGTCCGTGTTTCTGAACTCCATGATTTTT	13	0.325	No Hit
CGTTGGATCGATGTGGTGCTGCCGCGTTCTCTTCGGGCCGGGGCCTAAGCC	13	0.325	No Hit
CTACAATGATGGTAATATTTTTCGTCAACAGAGTTCACCTAGTGAGTGTTG	13	0.325	No Hit
AACCCGTAGATCCGATCTTGTAGATCGGAAGAGCACACGTCTGAACTCCAG	12	0.3	Illumina Multiplexing PCR Primer 2.01 (100% over 30bp)
ATCGATGTGGTGCTGCCGCGTTCTCTTCGGGCCGGGGCCTAAGCCGCGCCA	12	0.3	No Hit
CGCGACCTTAGATCGGAAGAGCACACGTCTGAACTCCAGTCACGGCTACAT	12	0.3	TruSeq Adapter, Index 11 (100% over 41bp)
NTGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAAC	12	0.3	No Hit
CGCGACCTCAGATCAGACGTGGCGACCCGCTGAATTTAAGCATATTAGTCA	11	0.27499999999999997	No Hit
CTTCGGTGATGACACGATGACGAGTCAGAAAGGCCACGTCTTGCTCTTGGG	11	0.27499999999999997	No Hit
AATACATGATGATCTCACTCCAACTTGAACTCTCTCACTGATTACTTGATG	10	0.25	No Hit
NTCGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGTAGATCG	10	0.25	No Hit
TTTGAGTGATGACTTGTATTGTCGGATACCCCTTCACCCCGTTCATGGGTG	10	0.25	No Hit
CACAAATGATGAATCCTTTGACGGGCGGGCAGAATCTCTGTGCTGATTGTC	10	0.25	No Hit
CACCCTGATGAACTGAATACCGCCCCAGTCTGATAGCTGTGGAGAAAGGTA	10	0.25	No Hit
AACAGTGATGACCACAACATGCCATCTGAGTTGCTTTGCTGAAATCCAGAG	9	0.22499999999999998	No Hit
TTCAAGTAATCCAGGATAGGCTAGATCGGAAGAGCACACGTCTGAACTCCA	9	0.22499999999999998	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
GACTCTTAGCGGTGGATCACTCGGCTCGTGCGAGATCGGAAGAGCACACGT	9	0.22499999999999998	No Hit
CTCCCATGATGTATAGCACTGGGCTCTGACCACCATGGTGACGTGTTGTGC	9	0.22499999999999998	No Hit
TTCTGTGATGAGACTCTGGAATGTGCGCTGGGCACAGCGCCCGAGTGAGAC	9	0.22499999999999998	No Hit
NGAAGAGCACACGTCTGAACTCCAGTCACGGCTACATCTCGTATGCCGTCT	9	0.22499999999999998	TruSeq Adapter, Index 11 (98% over 51bp)
CACAAGTGATGACTAACAAAGGGACTTAATACTGAAACATCTGGGATTGAA	9	0.22499999999999998	No Hit
TGTCCCATTTGACTACCACATTTCCTACTGATAATAGATTTTGGTGGCTAT	9	0.22499999999999998	No Hit
TCGATGTGGTGCTGCCGCGTTCTCTTCGGGCCGGGGCCTAAGCCGCGCCAG	8	0.2	No Hit
NTGCTGTGATGACTATCTTAGGACACCTTTGGATTAACCGTGAAATAAAAC	8	0.2	No Hit
GCGGATGATGACACCTGGGTATGCTGCACATTCCCGACTGCGTCGTGGGGA	8	0.2	No Hit
ACAAATGATGAATCCTTTGACGGGCGGGCAGAATCTCTGTGCTGATTGTCA	8	0.2	No Hit
GTGGAATGATGACATTCTCCGGAATCGCTGTACTGTCTTGATGAAAGTACT	8	0.2	No Hit
ATACTGTGATGATAACATAGTTCAGCAGACGTGACTCTGGTGAACAATCAT	8	0.2	No Hit
TGCAATGATGAAGGTTATGTTGGTAGGGACAACTGACCTTGTTGATGAATA	8	0.2	No Hit
CGCTGCGATCTATTGAAAGTCAGCCCTCGACACAAGGGTTTGTAGATCGGA	8	0.2	No Hit
NACAATGATGACTTATGTTTTTGCCGTTTACCCAGCTGAGGGTTTCTTTGA	7	0.17500000000000002	No Hit
NTACATGATGATCTCACTCCAACTTGAACTCTCTCACTGATTACTTGATGA	7	0.17500000000000002	No Hit
TAGCTTATCAGACTGATGTTGACAAGATCGGAAGAGCACACGTCTGAACTC	7	0.17500000000000002	Illumina Multiplexing PCR Primer 2.01 (100% over 27bp)
TAGCTTATCAGACTGATGTTGACTAGATCGGAAGAGCACACGTCTGAACTC	7	0.17500000000000002	Illumina Multiplexing PCR Primer 2.01 (100% over 27bp)
CCGCGACCTCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACGGCTACA	7	0.17500000000000002	TruSeq Adapter, Index 11 (100% over 40bp)
TATGGCACTGGTAGAATTCACTAGATCGGAAGAGCACACGTCTGAACTCCA	7	0.17500000000000002	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
TCAAGTGATGGAGACATGACCCTGGGAGGCCCATGACCGAACTGTCTGTAT	6	0.15	No Hit
GGTCAATGATGAGCTGGCATGTATTCTGAATCTAAAGTTGATTAATAACCA	6	0.15	No Hit
TTTGGCAATGGTAGAACTCACACCAGATCGGAAGAGCACACGTCTGAACTC	6	0.15	Illumina Multiplexing PCR Primer 2.01 (100% over 27bp)
GAAGCTGATGATACCACTTCTTTCCCCATCAGATCGACCCTGTTGATCTCA	6	0.15	No Hit
NGCGACCTCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACGGCTACAT	6	0.15	TruSeq Adapter, Index 11 (100% over 41bp)
CACCCGTAGAACCGACCTTGCGAGATCGGAAGAGCACACGTCTGAACTCCA	6	0.15	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
ACTTAATGATGATTGTCTATTTGATTGCTTGATGCAATGTGAAAACCACAT	6	0.15	No Hit
AACTGTGATGAAAGATTTGGTCTATACAGTAATATATTTATTACTAATGAG	6	0.15	No Hit
CGCGACCTCAGATCAGACGTAGATCGGAAGAGCACACGTCTGAACTCCAGT	6	0.15	Illumina Multiplexing PCR Primer 2.01 (100% over 31bp)
GGAGATGATGAACTGACTCAACCTGACCTTCAGAAATGGAGGTGATAAAAC	6	0.15	No Hit
CGCGACCTCAGATCAGAAGATCGGAAGAGCACACGTCTGAACTCCAGTCAC	5	0.125	Illumina Multiplexing PCR Primer 2.01 (100% over 34bp)
GATCGGAAGAGCACACGTCTGAACTCCAGTCACGGCTACATCTCGTATGCC	5	0.125	TruSeq Adapter, Index 11 (100% over 51bp)
GTGCAATGATGTGTGAATTTCTTCACCTGAAATCACCGTGAAGAGTAATTT	5	0.125	No Hit
CGCGACCTCAGATCAGACGTGGCGACCCGCTGAAAGATCGGAAGAGCACAC	5	0.125	No Hit
ACAAGTGATGACTAACAAAGGGACTTAATACTGAAACATCTGGGATTGAAA	5	0.125	No Hit
GACTCTTAGCGGTGGATCACTCGGCTCGTGCGTCGATGAAGAACGCAGCTA	5	0.125	No Hit
TTCTGTGATGACTTTACCCAAATGACTTTCGTTCTTCTGAGTTTGCTGAAG	5	0.125	No Hit
CGCGACCTCAGATCAGACAGATCGGAAGAGCACACGTCTGAACTCCAGTCA	5	0.125	Illumina Multiplexing PCR Primer 2.01 (100% over 33bp)
TGAGGTAGTAGATTGTATAGTTAGATCGGAAGAGCACACGTCTGAACTCCA	5	0.125	Illumina Multiplexing PCR Primer 2.01 (100% over 29bp)
TGAAGTGATGATTCACAATCATGTCTCTTCTCTGAAGAATTCTTGGAGATT	5	0.125	No Hit
TGGAGATGATGAACTGACTCAACCTGACCTTCAGAAATGGAGGTGATAAAA	5	0.125	No Hit
TTTCTGTGATGAATCAACTAGCTCACTATGACTGCTAATGAAAACACTGGA	5	0.125	No Hit
TAGCTTATCAGACTGATGTTGATAGATCGGAAGAGCACACGTCTGAACTCC	5	0.125	Illumina Multiplexing PCR Primer 2.01 (100% over 28bp)
>>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.275	0.0	0.0	0.0	0.0
2	9.35	0.0	0.0	0.0	0.0
3	9.35	0.0	0.0	0.0	0.0
4	9.425	0.0	0.0	0.0	0.0
5	9.475	0.0	0.0	0.0	0.0
6	9.525	0.0	0.0	0.0	0.0
7	9.55	0.0	0.0	0.0	0.0
8	9.775	0.0	0.0	0.0	0.0
9	10.1	0.0	0.0	0.0	0.0
10	12.5	0.0	0.0	0.0	0.0
11	13.125	0.0	0.0	0.0	0.0
12	13.55	0.0	0.0	0.0	0.0
13	14.05	0.0	0.0	0.0	0.0
14	14.425	0.0	0.0	0.0	0.0
15	14.95	0.0	0.0	0.0	0.0
16	15.425	0.0	0.0	0.0	0.0
17	15.75	0.0	0.0	0.0	0.0
18	16.325	0.0	0.0	0.0	0.0
19	17.0	0.0	0.0	0.0	0.0
20	17.55	0.0	0.0	0.0	0.0
21	18.275	0.0	0.0	0.0	0.0
22	20.4	0.0	0.0	0.0	0.0
23	24.075	0.0	0.0	0.0	0.0
24	26.975	0.0	0.0	0.0	0.0
25	28.375	0.0	0.0	0.0	0.0
26	28.8	0.0	0.0	0.0	0.0
27	29.375	0.0	0.0	0.0	0.0
28	29.875	0.0	0.0	0.0	0.0
29	30.275	0.0	0.0	0.0	0.0
30	30.7	0.0	0.0	0.0	0.0
31	31.325	0.0	0.0	0.0	0.0
32	32.3	0.0	0.0	0.0	0.0
33	33.25	0.0	0.0	0.0	0.0
34	33.7	0.0	0.0	0.0	0.0
35	34.0	0.0	0.0	0.0	0.0
36	34.15	0.0	0.0	0.0	0.0
37	34.4	0.0	0.0	0.0	0.0
38	34.625	0.0	0.0	0.0	0.0
39	35.0	0.0	0.0	0.0	0.0
>>END_MODULE
>>Kmer Content	warn
#Sequence	Count	PValue	Obs/Exp Max	Max Obs/Exp Position
GTGAAAT	20	3.6890357E-4	49.888885	1
GAAAACT	45	2.237357E-10	44.9	45
CTGATGA	20	6.2916183E-4	44.9	5
ATTTTTG	20	6.2916183E-4	44.9	45
ATGATGA	65	0.0	41.44615	5
ATCTCGT	75	0.0	35.920002	41
CATCTCG	75	0.0	35.920002	40
ACATCTC	75	0.0	35.920002	39
TACATCT	75	0.0	35.920002	38
TCTCGTA	75	0.0	35.920002	42
CGTATGC	75	0.0	35.920002	45
TCGTATG	75	0.0	35.920002	44
CTCGTAT	75	0.0	35.920002	43
CTACATC	75	0.0	35.920002	37
TGATGAG	25	0.0018869689	35.92	6
TGGTCCG	25	0.0018869689	35.92	21
GTGTTTC	25	0.0018869689	35.92	27
CCGTGTT	25	0.0018869689	35.92	25
CATTGGT	25	0.0018869689	35.92	18
TAAATGA	25	0.0018869689	35.92	2
>>END_MODULE
Rejected 2771430 READS because READLEN < 1
Read 2771430 spots for ERR4018148.sra
Written 2771430 spots for ERR4018148.sra
Rejected 2771430 READS because READLEN < 1
Read 2771430 spots for ERR4018148.sra
Written 2771430 spots for ERR4018148.sra
Rejected 2771430 READS because READLEN < 1
Read 2771430 spots for ERR4018148.sra
Written 2771430 spots for ERR4018148.sra
Rejected 2771430 READS because READLEN < 1
Read 2771430 spots for ERR4018148.sra
Written 2771430 spots for ERR4018148.sra
Rejected 2771430 READS because READLEN < 1
Read 2771430 spots for ERR4018148.sra
Written 2771430 spots for ERR4018148.sra
Rejected 2771430 READS because READLEN < 1
Read 2771430 spots for ERR4018148.sra
Written 2771430 spots for ERR4018148.sra
Rejected 2771432 READS because READLEN < 1
Read 2771432 spots for ERR4018148.sra
Written 2771432 spots for ERR4018148.sra
Rejected 2771430 READS because READLEN < 1
Read 2771430 spots for ERR4018148.sra
Written 2771430 spots for ERR4018148.sra
SRR ids: ['ERR4018148.sra']
extra args: ['--split-files', '--defline-qual', '+']
tempdir: /tmp/pfd_spxk2i50
ERR4018148.sra spots: 22171442
blocks: [[1, 2771430], [2771431, 5542860], [5542861, 8314290], [8314291, 11085720], [11085721, 13857150], [13857151, 16628580], [16628581, 19400010], [19400011, 22171442]]
ERR4018148 file size 3139468
ERR4018148 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 ERR4018148 ERR4018148_1.fastq
Input file:	ERR4018148_1.fastq
trimmed:	ERR4018148-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:02:36 2025 >> started

Sun Dec 28 15:02:42 2025 >> done (6.014s)
22171442 reads processed; of these:
 1608871 ( 7.26%) short reads filtered out after trimming by size control
 2201726 ( 9.93%) empty reads filtered out after trimming by size control
18360845 (82.81%) reads available; of these:
  250653 ( 1.37%) trimmed reads available after processing
18110192 (98.63%) untrimmed reads available after processing

Length distribution of reads after trimming:
length	count	percentage
 18	  135182	  0.74%
 19	  115368	  0.63%
 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	      45	  0.00%
 50	      58	  0.00%
 51	18110192	 98.63%
18360845 reads passed initial QC


criterion=sequence-density
sequence-density=20.58
sequence-density-rank=1
fanout-score=33.89
fanout-score-rank=1
prefix-density=21.41
prefix-fanout=32.6
sequence=AGATCGGAAGAGCACACGTCTGAACTCCAGTCACGGCTACATCTCGTATGCCGTCTTC


criterion=fanout-score
sequence-density=20.58
sequence-density-rank=1
fanout-score=33.89
fanout-score-rank=1
prefix-density=21.41
prefix-fanout=32.6
sequence=AGATCGGAAGAGCACACGTCTGAACTCCAGTCACGGCTACATCTCGTATGCCGTCTTC
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 AGATCGGAAGAGCACACGTCTGAACTCCAGTCACGGCTACATCTCGTATGCCGTCTTC -o ERR4018148 -
Input file:	STDIN
trimmed:	ERR4018148-trimmed.fastq

Parameters used:
-- 3' end adapter sequence (-x):	AGATCGGAAGAGCACACGTCTGAACTCCAGTCACGGCTACATCTCGTATGCCGTCTTC
-- 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:02:56 2025 >> started

Sun Dec 28 15:03:06 2025 >> done (10.255s)
16612193 reads processed; of these:
    7533 ( 0.05%) short reads filtered out after trimming by size control
    3709 ( 0.02%) empty reads filtered out after trimming by size control
16600951 (99.93%) reads available; of these:
 4883523 (29.42%) trimmed reads available after processing
11717428 (70.58%) untrimmed reads available after processing

Length distribution of reads after trimming:
length	count	percentage
 18	  119364	  0.72%
 19	  106414	  0.64%
 20	  168156	  1.01%
 21	  340113	  2.05%
 22	  771421	  4.65%
 23	  517069	  3.11%
 24	  323933	  1.95%
 25	  114039	  0.69%
 26	   83173	  0.50%
 27	   82035	  0.49%
 28	  108423	  0.65%
 29	   70974	  0.43%
 30	  121180	  0.73%
 31	  181994	  1.10%
 32	  151098	  0.91%
 33	   77275	  0.47%
 34	   81106	  0.49%
 35	   58686	  0.35%
 36	   85711	  0.52%
 37	   76908	  0.46%
 38	   49001	  0.30%
 39	   55179	  0.33%
 40	   34409	  0.21%
 41	   34040	  0.21%
 42	   57680	  0.35%
 43	  101167	  0.61%
 44	  147205	  0.89%
 45	  601792	  3.63%
 46	   92653	  0.56%
 47	   44017	  0.27%
 48	  248824	  1.50%
 49	      25	  0.00%
 50	      38	  0.00%
 51	11495849	 69.25%


criterion=sequence-density
sequence-density=4.63
sequence-density-rank=1
fanout-score=1.74
fanout-score-rank=12
prefix-density=0.02
prefix-fanout=1.7
sequence=GTGAAATGATGGCAATCATCTTTCGGGACTGACCTGAAATGAAGAGAAAACTCTCTGCTGATC


criterion=fanout-score
sequence-density=0.46
sequence-density-rank=12
fanout-score=41.76
fanout-score-rank=1
prefix-density=6.16
prefix-fanout=3.1
sequence=GTGATGACTTTCTTAGGAC
                                 Started job on |	Dec 28 15:03:18
                             Started mapping on |	Dec 28 15:03:18
                                    Finished on |	Dec 28 15:04:02
       Mapping speed, Million of reads per hour |	1501.33

                          Number of input reads |	18349603
                      Average input read length |	45
                                    UNIQUE READS:
                   Uniquely mapped reads number |	12584017
                        Uniquely mapped reads % |	68.58%
                          Average mapped length |	46.79
                       Number of splices: Total |	15577
            Number of splices: Annotated (sjdb) |	3858
                       Number of splices: GT/AG |	14656
                       Number of splices: GC/AG |	233
                       Number of splices: AT/AC |	1
               Number of splices: Non-canonical |	687
                      Mismatch rate per base, % |	0.38%
                         Deletion rate per base |	0.00%
                        Deletion average length |	1.08
                        Insertion rate per base |	0.00%
                       Insertion average length |	1.24
                             MULTI-MAPPING READS:
        Number of reads mapped to multiple loci |	3600857
             % of reads mapped to multiple loci |	19.62%
        Number of reads mapped to too many loci |	437268
             % of reads mapped to too many loci |	2.38%
                                  UNMAPPED READS:
       % of reads unmapped: too many mismatches |	0.00%
                 % of reads unmapped: too short |	9.31%
                     % of reads unmapped: other |	0.11%
                                  CHIMERIC READS:
                       Number of chimeric reads |	0
                            % of chimeric reads |	0.00%
N_unmapped	2164729	2164729	2164729
N_multimapping	3600857	3600857	3600857
N_noFeature	411467	545337	12074966
N_ambiguous	1028592	653354	58
UnstrandedReadsAssigned:11143958 PositiveStrandReadsAssigned:11385326 NegativeStrandReadsAssigned:508993
Dataset is classified positive stranded
MeadianReadLen=51 20thPercentileLength=33 echo kmer=29
ERR4018148 Starting Kallisto single end mapping to ensembl reference transcriptome. kmer=29

[quant] fragment length distribution is truncated gaussian with mean = 100, sd = 20
[index] k-mer length: 29
[index] number of targets: 31,715
[index] number of k-mers: 49,875,633
[index] number of equivalence classes: 74,027
[quant] running in single-end mode
[quant] will process file 1: ERR4018148-trimmed.fastq
[quant] finding pseudoalignments for the reads ... done
[quant] processed 18,349,603 reads, 132,235 reads pseudoaligned
[   em] quantifying the abundances ... done
[   em] the Expectation-Maximization algorithm ran for 702 rounds

  31716 ERR4018148.ke.tsv
  32883 ERR4018148.se.tsv
  64599 total
==> ERR4018148.ke.tsv <==
target_id	length	eff_length	est_counts	tpm
ENSRNOT00000047550.4	955	856	107	669.467
ENSRNOT00000040993.4	1039	940	140	797.663
ENSRNOT00000050156.3	1545	1446	1044	3866.8
ENSRNOT00000043693.3	684	585	321	2938.79
ENSRNOT00000046201.3	204	105	181	9232.27
ENSRNOT00000046108.3	681	582	85	782.195
ENSRNOT00000049683.3	784	685	97	758.403
ENSRNOT00000041241.3	348	249	63	1355.07
ENSRNOT00000044582.3	297	198	80	2163.93

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