Starting /dee2/code/volunteer_pipeline.sh ERR6133533 current disk space = 1544484282368 free memory = 1462825432 ERR6133533 SRAfilesize 6b74649ae7c8d5848d84a484c91a0658 ERR6133533.sra ERR6133533.sra file validated ERR6133533 is single end ERR6133533 is conventional basespace ERR6133533 read1 length is 70-93 nt ##FastQC 0.11.5 >>Basic Statistics pass #Measure Value Filename ERR6133533_1.fastq File type Conventional base calls Encoding Sanger / Illumina 1.9 Total Sequences 4000 Sequences flagged as poor quality 0 Sequence length 70-93 %GC 47 >>END_MODULE >>Per base sequence quality pass #Base Mean Median Lower Quartile Upper Quartile 10th Percentile 90th Percentile 1 33.72575 33.0 33.0 37.0 27.0 37.0 2 35.96325 37.0 37.0 37.0 33.0 37.0 3 35.16625 37.0 33.0 37.0 33.0 37.0 4 35.35 37.0 37.0 37.0 33.0 37.0 5 35.30525 37.0 37.0 37.0 33.0 37.0 6 35.63725 37.0 37.0 37.0 33.0 37.0 7 37.4885 37.0 37.0 40.0 33.0 40.0 8 37.395 37.0 37.0 40.0 33.0 40.0 9 37.406 37.0 37.0 40.0 33.0 40.0 10-11 37.444374999999994 37.0 37.0 40.0 33.0 40.0 12-13 37.345 37.0 37.0 40.0 33.0 40.0 14-15 37.326125 37.0 37.0 40.0 33.0 40.0 16-17 37.302375 37.0 37.0 40.0 33.0 40.0 18-19 37.1795 37.0 37.0 40.0 33.0 40.0 20-21 37.095625 37.0 37.0 40.0 33.0 40.0 22-23 37.033375 37.0 37.0 40.0 33.0 40.0 24-25 37.136625 37.0 37.0 40.0 33.0 40.0 26-27 37.039375 37.0 37.0 40.0 33.0 40.0 28-29 36.92375 37.0 37.0 40.0 33.0 40.0 30-31 36.779624999999996 37.0 37.0 40.0 33.0 40.0 32-33 36.72025 37.0 37.0 40.0 33.0 40.0 34-35 36.514375 37.0 37.0 40.0 33.0 40.0 36-37 36.3035 37.0 37.0 40.0 33.0 40.0 38-39 36.298874999999995 37.0 37.0 40.0 33.0 40.0 40-41 36.02425 37.0 37.0 40.0 33.0 40.0 42-43 35.972125 37.0 37.0 40.0 33.0 40.0 44-45 35.65925 37.0 33.0 38.5 33.0 40.0 46-47 35.330375000000004 37.0 33.0 37.0 30.0 40.0 48-49 35.382000000000005 37.0 33.0 37.0 33.0 40.0 50-51 35.25325 37.0 33.0 37.0 30.0 40.0 52-53 34.938375 37.0 33.0 37.0 30.0 40.0 54-55 34.876125 37.0 33.0 37.0 30.0 37.0 56-57 34.81075 37.0 33.0 37.0 33.0 37.0 58-59 34.115625 37.0 33.0 37.0 27.0 37.0 60-61 34.34375 37.0 33.0 37.0 27.0 37.0 62-63 34.171 37.0 33.0 37.0 27.0 37.0 64-65 34.20575 37.0 33.0 37.0 27.0 37.0 66-67 34.293375 37.0 33.0 37.0 27.0 37.0 68-69 33.402625 35.0 33.0 37.0 27.0 37.0 70-71 33.53329429429429 35.0 33.0 37.0 27.0 37.0 72-73 33.996778177529166 37.0 33.0 37.0 27.0 37.0 74-75 34.02547690763052 37.0 33.0 37.0 27.0 37.0 76-77 33.85894338606117 37.0 33.0 37.0 27.0 37.0 78-79 34.01131182427939 37.0 33.0 37.0 27.0 37.0 80-81 33.875278647916915 37.0 33.0 37.0 27.0 37.0 82-83 33.794053063495824 37.0 33.0 37.0 27.0 37.0 84-85 33.76647653652636 37.0 33.0 37.0 27.0 37.0 86-87 33.571645954856706 37.0 33.0 37.0 27.0 37.0 88-89 33.6596500126807 37.0 33.0 37.0 27.0 37.0 90-91 33.46525488206949 37.0 33.0 37.0 27.0 37.0 92-93 33.42264773015471 37.0 33.0 37.0 27.0 37.0 >>END_MODULE >>Per sequence quality scores pass #Quality Count 20 12.0 21 12.0 22 27.0 23 29.0 24 19.0 25 36.0 26 40.0 27 49.0 28 62.0 29 62.0 30 90.0 31 113.0 32 156.0 33 183.0 34 258.0 35 414.0 36 924.0 37 862.0 38 646.0 39 6.0 >>END_MODULE >>Per base sequence content fail #Base G A T C 1 90.64999999999999 2.4 1.675 5.2749999999999995 2 77.225 14.149999999999999 4.45 4.175 3 37.025000000000006 39.75 12.875 10.35 4 34.975 27.800000000000004 18.2 19.025 5 26.924999999999997 32.675 22.5 17.9 6 20.625 40.699999999999996 23.225 15.45 7 41.175 28.075 17.275 13.475000000000001 8 30.525000000000002 31.85 20.974999999999998 16.650000000000002 9 28.725 28.375 24.925 17.974999999999998 10-11 27.3125 27.987499999999997 26.6625 18.0375 12-13 30.5 25.637500000000003 25.474999999999998 18.387500000000003 14-15 24.3625 30.65 26.637499999999996 18.35 16-17 25.8125 31.2375 24.5375 18.4125 18-19 24.34054256782098 26.653331666458307 25.91573946743343 23.090386298287285 20-21 26.147305239464803 25.697136426159812 27.747905464549206 20.407652869826183 22-23 27.5875 23.3 26.875 22.237499999999997 24-25 27.725 23.5 26.937499999999996 21.837500000000002 26-27 26.337500000000002 24.525 28.4125 20.724999999999998 28-29 26.753344168021005 25.453181647705964 27.128391048881113 20.665083135391924 30-31 27.675 24.9875 25.8 21.5375 32-33 25.637500000000003 27.212500000000002 25.7 21.45 34-35 26.773426748404855 24.55898911547604 26.973601901663958 21.69398223445515 36-37 25.7875 24.05 27.875 22.287499999999998 38-39 27.0125 23.6625 29.1875 20.1375 40-41 26.932699524643482 25.381536152114087 26.294721040780583 21.391043282461847 42-43 25.86055826761797 26.736763049192643 26.5114532482163 20.891225434973087 44-45 24.133616914800452 25.284624046040282 29.438258476166645 21.143500562992617 46-47 25.4 22.85 28.6125 23.1375 48-49 25.828228528566072 23.6029503687961 29.466183272909113 21.102637829728714 50-51 26.2125 25.587500000000002 27.3 20.9 52-53 25.672463405479796 25.29713499311898 26.610784436381834 22.419617165019393 54-55 24.1625 26.650000000000002 28.0625 21.125 56-57 27.1 25.1 27.35 20.45 58-59 24.212500000000002 25.5625 28.1625 22.0625 60-61 25.2125 24.4875 28.4125 21.8875 62-63 24.1875 25.45 30.875000000000004 19.4875 64-65 23.8375 27.675 27.8375 20.65 66-67 25.55 25.937500000000004 28.025 20.4875 68-69 23.3375 26.3125 28.225 22.125 70-71 25.15007503751876 25.41270635317659 28.12656328164082 21.310655327663834 72-73 26.059162697417896 23.82802707445475 28.95462521935322 21.15818500877413 74-75 23.94578313253012 26.44327309236948 28.865461847389557 20.745481927710845 76-77 23.31951250157055 25.304686518406832 28.86040959919588 22.51539138082674 78-79 25.04089593557317 25.2422297722411 29.306656599974833 20.410217692210896 80-81 24.48002016891466 27.80789108786084 29.54745997730997 18.164628765914532 82-83 24.822964087000507 25.025290844714217 28.566009104704097 21.585735963581186 84-85 24.175963488843813 23.47870182555781 31.00912778904665 21.336206896551722 86-87 23.9158001521684 25.61501394876997 31.118437737763127 19.350748161298505 88-89 23.218361653563278 26.65483134669034 30.129343139741316 19.99746386000507 90-91 26.388536647222928 25.589652548820695 29.03880294192239 18.983007862033986 92-93 22.939386254121228 28.024346943951308 28.759827542480345 20.27643925944712 >>END_MODULE >>Per sequence GC content warn #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 1.0 17 5.0 18 5.0 19 1.5 20 0.5 21 3.0 22 4.0 23 5.5 24 7.0 25 6.0 26 7.5 27 10.0 28 16.5 29 21.0 30 19.0 31 19.0 32 31.0 33 40.5 34 48.0 35 65.5 36 92.0 37 121.5 38 145.5 39 164.0 40 176.0 41 181.5 42 193.0 43 203.5 44 203.0 45 198.0 46 211.0 47 187.5 48 154.5 49 158.5 50 151.5 51 152.0 52 146.0 53 167.5 54 176.5 55 117.0 56 74.0 57 69.0 58 72.0 59 64.0 60 50.0 61 40.0 62 41.0 63 41.5 64 38.5 65 40.0 66 35.5 67 33.0 68 26.0 69 17.0 70 10.5 71 7.0 72 8.0 73 7.0 74 5.0 75 5.5 76 5.5 77 3.5 78 1.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 pass #Base N-Count 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-11 0.0 12-13 0.0 14-15 0.0 16-17 0.0 18-19 0.0125 20-21 0.0375 22-23 0.0 24-25 0.0 26-27 0.0 28-29 0.0125 30-31 0.0 32-33 0.0 34-35 0.08750000000000001 36-37 0.0 38-39 0.0 40-41 0.075 42-43 0.13749999999999998 44-45 0.08750000000000001 46-47 0.0 48-49 0.0125 50-51 0.0 52-53 0.08750000000000001 54-55 0.0 56-57 0.0 58-59 0.0 60-61 0.0 62-63 0.0 64-65 0.0 66-67 0.0 68-69 0.0 70-71 0.0 72-73 0.0 74-75 0.0 76-77 0.0 78-79 0.0 80-81 0.0 82-83 0.0 84-85 0.0 86-87 0.0 88-89 0.0 90-91 0.0 92-93 0.0 >>END_MODULE >>Sequence Length Distribution warn #Length Count 70 4.0 71 4.0 72 6.0 73 2.0 74 0.0 75 3.0 76 3.0 77 3.0 78 3.0 79 3.0 80 5.0 81 6.0 82 8.0 83 5.0 84 2.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 3943.0 >>END_MODULE >>Sequence Duplication Levels pass #Total Deduplicated Percentage 84.175 #Duplication Level Percentage of deduplicated Percentage of total 1 93.88179388179388 79.025 2 3.890703890703891 6.550000000000001 3 0.8910008910008911 2.25 4 0.3564003564003564 1.2 5 0.2079002079002079 0.8750000000000001 6 0.0891000891000891 0.44999999999999996 7 0.029700029700029697 0.17500000000000002 8 0.11880011880011879 0.8 9 0.11880011880011879 0.8999999999999999 >10 0.3861003861003861 6.25 >50 0.029700029700029697 1.525 >100 0.0 0.0 >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 GGGAGAGCAATACAAGCGTTGTGCTGCTAGGCGAAGCGGTTGAGTGCCGC 61 1.525 No Hit GGGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGT 43 1.075 No Hit GGGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTGG 37 0.9249999999999999 No Hit GGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTT 28 0.7000000000000001 No Hit GGATTCAATTTCAACCAATCTGTAGTTGATAGTCAAGGTCGCGTTATTAA 20 0.5 No Hit GGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGTG 19 0.475 No Hit GGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTGGG 16 0.4 No Hit GGGGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTG 15 0.375 No Hit GGCCTGTAGTAGGAATCTGGTTCACTGCTTTAGGTATTAGTACTATGGCG 14 0.35000000000000003 No Hit GGTCGCGTTATTAATACTTGGGCTGATATCATCAACCGTGCTAATCTTGG 14 0.35000000000000003 No Hit GGCGTTCAACCTAAATGGATTCAATTTCAACCAATCTGTAGTTGATAGTC 12 0.3 No Hit GGGTGTGAGCTGGAGAGACGATCGGGTCTCTCAGCCGGCGTCTTCATCAG 11 0.27499999999999997 No Hit GGAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAA 11 0.27499999999999997 No Hit GTAGTAGGAATCTGGTTCACTGCTTTAGGTATTAGTACTATGGCGTTCAA 10 0.25 No Hit GGTGGAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGC 9 0.22499999999999998 No Hit TACCTAGGCACCCAGAGACGAGGAAGGGCGTAGCAAGCGACGAAATGCTT 9 0.22499999999999998 No Hit GGTTTTGAAAAGGGAATCGATCGTGATTTAGAACCTGTTCTTTACATGAA 9 0.22499999999999998 No Hit GGAAAAGAAAGCAAAAGCGATTCCCGTAGTAGCGGCGAGCGAAATGGGAG 9 0.22499999999999998 No Hit GAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAAT 8 0.2 No Hit GGAAGAGCCCGAGCTGCTGCAGCAGGTTTTGAAAAGGGAATCGATCGTGA 8 0.2 No Hit GGAGTATCCTTGATATATAAATATATAGATAAATAGATTTAAATGGAAAA 8 0.2 No Hit GGTGGGGCAATGGGCATTCCGTTGAGTTTGTATGGTGGGTGCTGTTTTGC 8 0.2 No Hit GGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAATCG 7 0.17500000000000002 No Hit GGGTGAGCATATATATTTATACGACGAATAAAAGGCTGCCACGTGGCGGG 6 0.15 No Hit GGAGAACGGCGGCTGCGGCTGCAGCACCTGCAAGTGCGGCACCAGCTGCG 6 0.15 No Hit GGATTTGAAGAAAAAAAAGACTTCGATTCATTTTCTATTTATTTCGTTAG 6 0.15 No Hit GGGAGTATTATGAAAGGAGATTAATCATAGATTATCAAAAACCCTAGAAA 5 0.125 No Hit GGAAAAAAAGAAAAGAGACGTATATCAAACTGAATATACAGCATTCCATG 5 0.125 No Hit GGAACAAACGAGGATAAGATAAAATTTCTTTAAATTTATTTTGCCCAAGT 5 0.125 No Hit GGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGTGTACAAGCTCGTAA 5 0.125 No Hit GGGCTGATATCATCAACCGTGCTAATCTTGGTATGGAAGTAATGCACGAA 5 0.125 No Hit GGGAATCGATCGTGATTTAGAACCTGTTCTTTACATGAACCCTCTTAACT 5 0.125 No Hit GGGGAAGTACACCAGCGACGGCGAGGCCGCCGCCGCCAAGGAAGGCATGT 5 0.125 No Hit >>END_MODULE >>Adapter Content pass #Position Illumina Universal Adapter Illumina Small RNA 3' Adapter Illumina Small RNA 5' Adapter Nextera Transposase Sequence SOLID Small RNA Adapter 1 0.0 0.0 0.0 0.0 0.0 2 0.0 0.0 0.0 0.0 0.0 3 0.0 0.0 0.0 0.0 0.0 4 0.0 0.0 0.0 0.0 0.0 5 0.0 0.0 0.0 0.0 0.0 6 0.0 0.0 0.0 0.0 0.0 7 0.0 0.0 0.0 0.0 0.0 8 0.0 0.0 0.0 0.0 0.0 9 0.0 0.0 0.0 0.0 0.0 10-11 0.0 0.0 0.0 0.0 0.0 12-13 0.0 0.0 0.0 0.0 0.0 14-15 0.0 0.0 0.0 0.0 0.0 16-17 0.0 0.0 0.0 0.0 0.0 18-19 0.0 0.0 0.0 0.0 0.0 20-21 0.0 0.0 0.0 0.0 0.0 22-23 0.0 0.0 0.0 0.0 0.0 24-25 0.0 0.0 0.0 0.0 0.0 26-27 0.0 0.0 0.0 0.0 0.0 28-29 0.0 0.0 0.0 0.0 0.0 30-31 0.0 0.0 0.0 0.0 0.0 32-33 0.0 0.0 0.0 0.0 0.0 34-35 0.0 0.0 0.0 0.0 0.0 36-37 0.0 0.0 0.0 0.0 0.0 38-39 0.0 0.0 0.0 0.0 0.0 40-41 0.0 0.0 0.0 0.0 0.0 42-43 0.0 0.0 0.0 0.0 0.0 44-45 0.0 0.0 0.0 0.0 0.0 46-47 0.0 0.0 0.0 0.0 0.0 48-49 0.0 0.0 0.0 0.0 0.0 50-51 0.0 0.0 0.0 0.0 0.0 52-53 0.0 0.0 0.0 0.0 0.0 54-55 0.0 0.0 0.0 0.0 0.0 56-57 0.0 0.0 0.0 0.0 0.0 58-59 0.0 0.0 0.0 0.0 0.0 60-61 0.0 0.0 0.0 0.0 0.0 62-63 0.0 0.0 0.0 0.0 0.0 64-65 0.0 0.0 0.0 0.0 0.0 66-67 0.0 0.0 0.0 0.0 0.0 68-69 0.0 0.0 0.0 0.0 0.0 70-71 0.0 0.0 0.0 0.0 0.0 72-73 0.0 0.0 0.0 0.0 0.0 74-75 0.0 0.0 0.0 0.0 0.0 76-77 0.0 0.0 0.0 0.0 0.0 78-79 0.0 0.0 0.0 0.0 0.0 80-81 0.0 0.0 0.0 0.0 0.0 >>END_MODULE >>Kmer Content warn #Sequence Count PValue Obs/Exp Max Max Obs/Exp Position TGCACCT 20 0.0027298967 65.25 7 CACCTGG 20 0.0027298967 65.25 9 ATGCACC 20 0.0027298967 65.25 6 GCACCTG 20 0.0027298967 65.25 8 GGGAAAT 20 0.0027298967 65.25 1 AAATGCA 20 0.0027298967 65.25 4 GGAAATG 20 0.0027298967 65.25 2 AATGCAC 25 0.006606883 52.2 5 CTTGCTC 20 7.824886E-4 43.5 76-77 TCTTGCT 25 0.0023441382 34.8 76-77 >>END_MODULE Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282316 READS because READLEN < 1 Read 282316 spots for ERR6133533.sra Written 282316 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra Rejected 282306 READS because READLEN < 1 Read 282306 spots for ERR6133533.sra Written 282306 spots for ERR6133533.sra SRR ids: ['ERR6133533.sra'] extra args: ['--split-files', '--defline-qual', '+'] tempdir: /tmp/pfd_8bp8vi9x ERR6133533.sra spots: 5646130 blocks: [[1, 282306], [282307, 564612], [564613, 846918], [846919, 1129224], [1129225, 1411530], [1411531, 1693836], [1693837, 1976142], [1976143, 2258448], [2258449, 2540754], [2540755, 2823060], [2823061, 3105366], [3105367, 3387672], [3387673, 3669978], [3669979, 3952284], [3952285, 4234590], [4234591, 4516896], [4516897, 4799202], [4799203, 5081508], [5081509, 5363814], [5363815, 5646130]] ERR6133533 file size 1252428 ERR6133533 completed basic pipeline successfully skewer v0.2.2 [April 4, 2016] COMMAND LINE: skewer -f sanger -l 18 -q 10 -k inf -t 20 -o ERR6133533 ERR6133533_1.fastq Input file: ERR6133533_1.fastq trimmed: ERR6133533-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): 20 Sat Dec 7 08:06:25 2024 >> started Sat Dec 7 08:06:28 2024 >> done (2.873s) 5646130 reads processed; of these: 49 ( 0.00%) short reads filtered out after trimming by size control 10 ( 0.00%) empty reads filtered out after trimming by size control 5646071 (100.00%) reads available; of these: 99182 ( 1.76%) trimmed reads available after processing 5546889 (98.24%) untrimmed reads available after processing Length distribution of reads after trimming: length count percentage 18 5 0.00% 19 8 0.00% 20 8 0.00% 21 3 0.00% 22 9 0.00% 23 2 0.00% 24 1 0.00% 25 3 0.00% 26 3 0.00% 27 3 0.00% 28 11 0.00% 29 28 0.00% 30 7 0.00% 31 7 0.00% 32 4 0.00% 33 5 0.00% 34 15 0.00% 35 45 0.00% 36 408 0.01% 37 10 0.00% 38 9 0.00% 39 20 0.00% 40 19 0.00% 41 12 0.00% 42 8 0.00% 43 11 0.00% 44 10 0.00% 45 9 0.00% 46 6 0.00% 47 8 0.00% 48 9 0.00% 49 6 0.00% 50 7 0.00% 51 10 0.00% 52 4 0.00% 53 1 0.00% 54 5 0.00% 55 5 0.00% 56 1 0.00% 57 2 0.00% 58 5 0.00% 59 3 0.00% 60 6 0.00% 61 6 0.00% 62 0 0.00% 63 5 0.00% 64 6 0.00% 65 10 0.00% 66 11 0.00% 67 18 0.00% 68 51 0.00% 69 75 0.00% 70 5271 0.09% 71 5283 0.09% 72 5424 0.10% 73 4992 0.09% 74 5029 0.09% 75 5205 0.09% 76 4876 0.09% 77 5160 0.09% 78 5554 0.10% 79 5962 0.11% 80 5715 0.10% 81 5996 0.11% 82 6574 0.12% 83 6765 0.12% 84 5948 0.11% 85 334 0.01% 86 530 0.01% 87 914 0.02% 88 1507 0.03% 89 2707 0.05% 90 5754 0.10% 91 16431 0.29% 92 68437 1.21% 93 5464760 96.79% 5646071 reads passed initial QC criterion=sequence-density sequence-density=0.92 sequence-density-rank=1 fanout-score=3.79 fanout-score-rank=25 prefix-density=1.14 prefix-fanout=3.1 sequence=ATGCATGCATGC criterion=fanout-score sequence-density=0.02 sequence-density-rank=24 fanout-score=78.71 fanout-score-rank=1 prefix-density=0.28 prefix-fanout=5.1 sequence=AGGAAGAGGAGGATGCAGTCAGGGCTCACAAACAACCTGCCACTGCCGCCATTGGCCTTCTAAAACAGGAGAGGAGGGGTTAGATAGTTTCGATCTGCAAGGGGGTGGGGCAATGGGCATTCCGTTGAGTTTGTATGGTGGGTGCTGTTTTGCAGAAGCTGTATGCTTATGAGCAGCTATGGTACTTATCGAGGACAGTAGCGTACTTGAGGTGTTGTATTTT Started job on | Dec 07 08:06:47 Started mapping on | Dec 07 08:06:47 Finished on | Dec 07 08:06:54 Mapping speed, Million of reads per hour | 2903.69 Number of input reads | 5646071 Average input read length | 92 UNIQUE READS: Uniquely mapped reads number | 4452287 Uniquely mapped reads % | 78.86% Average mapped length | 92.43 Number of splices: Total | 286269 Number of splices: Annotated (sjdb) | 241042 Number of splices: GT/AG | 278740 Number of splices: GC/AG | 5108 Number of splices: AT/AC | 98 Number of splices: Non-canonical | 2323 Mismatch rate per base, % | 0.41% Deletion rate per base | 0.04% Deletion average length | 1.77 Insertion rate per base | 0.02% Insertion average length | 1.98 MULTI-MAPPING READS: Number of reads mapped to multiple loci | 1099554 % of reads mapped to multiple loci | 19.47% Number of reads mapped to too many loci | 21657 % of reads mapped to too many loci | 0.38% UNMAPPED READS: % of reads unmapped: too many mismatches | 0.00% % of reads unmapped: too short | 1.25% % of reads unmapped: other | 0.03% CHIMERIC READS: Number of chimeric reads | 0 % of chimeric reads | 0.00% N_unmapped 94230 94230 94230 N_multimapping 1099554 1099554 1099554 N_noFeature 228150 267011 4264352 N_ambiguous 166371 17388 384 UnstrandedReadsAssigned:4057766 PositiveStrandReadsAssigned:4167888 NegativeStrandReadsAssigned:187551 Dataset is classified positive stranded MeadianReadLen=93 20thPercentileLength=93 echo kmer=89 ERR6133533 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: 52,972 [index] number of k-mers: 66,720,672 [index] number of equivalence classes: 111,837 [quant] running in single-end mode [quant] will process file 1: ERR6133533-trimmed.fastq [quant] finding pseudoalignments for the reads ... done [quant] processed 5,646,071 reads, 4,930,431 reads pseudoaligned [ em] quantifying the abundances ... done [ em] the Expectation-Maximization algorithm ran for 1,075 rounds 52973 ERR6133533.ke.tsv 35125 ERR6133533.se.tsv 88098 total ==> ERR6133533.ke.tsv <== target_id length eff_length est_counts tpm PNS24245 936 837 0 0 PNS24247 1044 945 0 0 PNS24249 1928 1829 0 0 PNS24246 1044 945 0 0 PNS24248 1044 945 0 0 PNS24244 1471 1372 139 27.0044 PNS24243 293 194 0 0 KQK14069 1603 1504 183 32.4323 KQK14071 474 375 0 0 ==> ERR6133533.se.tsv <== BRADI_1g14170v3 186 BRADI_1g53295v3 44 BRADI_1g59795v3 39 BRADI_1g07683v3 0 BRADI_1g00485v3 1 BRADI_1g20270v3 120 BRADI_1g74790v3 59 BRADI_1g09890v3 0 BRADI_1g77505v3 161 BRADI_1g48960v3 0 ERR6133533 completed mapping pipeline successfully