Starting /dee2/code/volunteer_pipeline.sh ERR6133556 current disk space = 1544482127872 free memory = 1598270452 ERR6133556 SRAfilesize 95479c5f7df4b370c05f743ef87913a5 ERR6133556.sra ERR6133556.sra file validated ERR6133556 is single end ERR6133556 is conventional basespace ERR6133556 read1 length is 70-93 nt ##FastQC 0.11.5 >>Basic Statistics pass #Measure Value Filename ERR6133556_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 48 >>END_MODULE >>Per base sequence quality pass #Base Mean Median Lower Quartile Upper Quartile 10th Percentile 90th Percentile 1 35.43675 37.0 33.0 37.0 33.0 37.0 2 36.436 37.0 37.0 37.0 37.0 37.0 3 35.84475 37.0 37.0 37.0 33.0 37.0 4 35.26175 37.0 37.0 37.0 33.0 37.0 5 35.0675 37.0 37.0 37.0 33.0 37.0 6 35.546 37.0 37.0 37.0 33.0 37.0 7 37.29125 37.0 37.0 40.0 33.0 40.0 8 37.24425 37.0 37.0 40.0 33.0 40.0 9 37.36025 37.0 37.0 40.0 33.0 40.0 10-11 37.279125 37.0 37.0 40.0 33.0 40.0 12-13 37.2615 37.0 37.0 40.0 33.0 40.0 14-15 37.17575 37.0 37.0 40.0 33.0 40.0 16-17 37.151125 37.0 37.0 40.0 33.0 40.0 18-19 37.119625 37.0 37.0 40.0 33.0 40.0 20-21 36.8705 37.0 37.0 40.0 33.0 40.0 22-23 36.84075 37.0 37.0 40.0 33.0 40.0 24-25 36.9965 37.0 37.0 40.0 33.0 40.0 26-27 36.925250000000005 37.0 37.0 40.0 33.0 40.0 28-29 36.876 37.0 37.0 40.0 33.0 40.0 30-31 36.7175 37.0 37.0 40.0 33.0 40.0 32-33 36.496375 37.0 37.0 40.0 33.0 40.0 34-35 36.342 37.0 37.0 40.0 33.0 40.0 36-37 36.18325 37.0 37.0 40.0 33.0 40.0 38-39 36.00375 37.0 37.0 40.0 33.0 40.0 40-41 35.878125 37.0 35.0 40.0 33.0 40.0 42-43 35.756125 37.0 33.0 40.0 33.0 40.0 44-45 35.503625 37.0 33.0 37.0 33.0 40.0 46-47 35.287625 37.0 33.0 37.0 30.0 40.0 48-49 35.16175 37.0 33.0 37.0 27.0 40.0 50-51 34.92275 37.0 33.0 37.0 27.0 40.0 52-53 34.574875000000006 37.0 33.0 37.0 27.0 40.0 54-55 34.499875 37.0 33.0 37.0 27.0 37.0 56-57 34.398624999999996 37.0 33.0 37.0 27.0 37.0 58-59 33.193625 35.0 33.0 37.0 27.0 37.0 60-61 33.854749999999996 37.0 33.0 37.0 27.0 37.0 62-63 33.770624999999995 37.0 33.0 37.0 27.0 37.0 64-65 33.798 37.0 33.0 37.0 27.0 37.0 66-67 33.89575 37.0 33.0 37.0 27.0 37.0 68-69 33.261125 35.0 33.0 37.0 27.0 37.0 70-71 33.292496871871876 35.0 33.0 37.0 27.0 37.0 72-73 33.711929481104406 37.0 33.0 37.0 27.0 37.0 74-75 33.389256869830035 37.0 33.0 37.0 27.0 37.0 76-77 33.59907036868679 37.0 33.0 37.0 27.0 37.0 78-79 33.599821882873904 37.0 33.0 37.0 27.0 37.0 80-81 33.39917857421342 37.0 33.0 37.0 27.0 37.0 82-83 33.263912748891485 37.0 33.0 37.0 27.0 37.0 84-85 33.28040506513635 37.0 33.0 37.0 27.0 37.0 86-87 33.03105906313645 37.0 33.0 37.0 27.0 37.0 88-89 33.18622708757638 37.0 33.0 37.0 27.0 37.0 90-91 33.05689918533605 35.0 33.0 37.0 27.0 37.0 92-93 32.88505600814664 35.0 33.0 37.0 27.0 37.0 >>END_MODULE >>Per sequence quality scores pass #Quality Count 20 23.0 21 16.0 22 31.0 23 22.0 24 28.0 25 47.0 26 49.0 27 53.0 28 55.0 29 91.0 30 79.0 31 116.0 32 143.0 33 212.0 34 262.0 35 463.0 36 896.0 37 886.0 38 517.0 39 11.0 >>END_MODULE >>Per base sequence content fail #Base G A T C 1 88.64999999999999 3.075 2.45 5.825 2 73.45 16.275000000000002 5.65 4.625 3 39.225 38.275 11.825 10.674999999999999 4 34.825 28.95 16.125 20.1 5 26.625 33.074999999999996 23.025000000000002 17.275 6 21.099999999999998 39.475 22.7 16.725 7 35.675000000000004 29.975 18.25 16.1 8 30.825000000000003 27.950000000000003 23.325000000000003 17.9 9 26.224999999999998 28.199999999999996 26.5 19.075 10-11 26.85 26.924999999999997 26.187500000000004 20.0375 12-13 28.925 25.5625 24.9375 20.575 14-15 23.674999999999997 30.4 26.3125 19.6125 16-17 26.187500000000004 29.8875 23.4125 20.5125 18-19 25.474999999999998 26.825 24.587500000000002 23.1125 20-21 25.85 26.125 26.375 21.65 22-23 28.012500000000003 22.95 25.8625 23.175 24-25 27.8625 23.3875 25.937500000000004 22.8125 26-27 26.387500000000003 25.9625 26.5375 21.1125 28-29 27.0125 27.35 24.25 21.3875 30-31 29.049999999999997 25.912499999999998 22.975 22.0625 32-33 26.737499999999997 28.8625 23.25 21.15 34-35 26.787499999999998 25.0625 25.637500000000003 22.5125 36-37 25.912499999999998 23.75 27.1 23.2375 38-39 26.975987993997 26.013006503251624 26.663331665832917 20.34767383691846 40-41 27.633224918689013 25.11883912934701 25.99449587190393 21.253440080060045 42-43 25.99074884360545 28.816102012751593 24.028003500437556 21.165145643205403 44-45 24.637500000000003 25.337500000000002 27.9375 22.0875 46-47 26.4625 22.8125 25.724999999999998 25.0 48-49 27.1375 23.5625 28.125 21.175 50-51 26.375 26.1125 26.75 20.7625 52-53 26.143627020929944 26.757739065045744 24.326356686301544 22.772277227722775 54-55 24.1625 27.85 26.8375 21.15 56-57 26.387500000000003 26.0125 26.025 21.575 58-59 25.0 25.924999999999997 26.525 22.55 60-61 26.7625 25.0125 25.7375 22.4875 62-63 23.9875 26.7625 29.875 19.375 64-65 25.4 27.487499999999997 25.775 21.337500000000002 66-67 25.374999999999996 27.3875 26.4125 20.825 68-69 24.1375 26.737499999999997 26.525 22.6 70-71 25.46273136568284 25.46273136568284 25.975487743871934 23.09904952476238 72-73 26.323049912214696 23.952846751943817 27.727614747930772 21.996488587910708 74-75 24.28697072496545 27.880386983289355 27.3903756753361 20.442266616409096 76-77 23.608662805338707 26.202467892218586 27.209770838579704 22.979098463863007 78-79 24.738236407215844 25.747445439636685 27.601867036709976 21.91245111643749 80-81 24.90503925044315 27.969106102810837 28.247657634844263 18.878197011901747 82-83 25.669671194617237 24.14624857179129 28.10714739113876 22.07693284245271 84-85 24.16995293219692 23.991858542170206 28.558707543569522 23.27948098206335 86-87 23.19246435845214 26.883910386965375 30.346232179226067 19.577393075356415 88-89 23.039714867617107 27.724032586558046 29.531568228105908 19.70468431771894 90-91 26.985743380855396 26.845723014256617 26.40020366598778 19.7683299389002 92-93 23.20519348268839 28.385947046843178 27.864052953156822 20.54480651731161 >>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.5 16 0.5 17 4.0 18 5.0 19 1.0 20 3.5 21 6.0 22 4.0 23 6.5 24 8.0 25 5.5 26 8.5 27 13.0 28 16.5 29 18.0 30 22.0 31 23.5 32 21.0 33 34.0 34 51.0 35 58.5 36 67.5 37 95.5 38 127.5 39 131.5 40 139.0 41 160.5 42 171.5 43 174.0 44 177.0 45 173.5 46 215.0 47 228.5 48 180.5 49 175.5 50 165.0 51 160.0 52 160.5 53 161.0 54 163.5 55 116.5 56 85.5 57 90.5 58 82.0 59 63.5 60 68.0 61 79.0 62 69.0 63 57.0 64 46.5 65 44.0 66 42.0 67 35.5 68 27.5 69 21.0 70 16.0 71 12.0 72 9.5 73 4.5 74 3.0 75 2.5 76 1.0 77 1.0 78 1.0 79 0.5 80 0.0 81 0.0 82 0.0 83 0.0 84 0.0 85 0.0 86 0.0 87 0.0 88 0.0 89 0.0 90 0.0 91 0.0 92 0.0 93 0.0 94 0.0 95 0.0 96 0.0 97 0.0 98 0.0 99 0.0 100 0.0 >>END_MODULE >>Per base N content 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.0 20-21 0.0 22-23 0.0 24-25 0.0 26-27 0.0 28-29 0.0 30-31 0.0 32-33 0.0 34-35 0.0 36-37 0.0 38-39 0.05 40-41 0.075 42-43 0.0125 44-45 0.0 46-47 0.0 48-49 0.0 50-51 0.0 52-53 0.2625 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 7.0 72 4.0 73 4.0 74 3.0 75 4.0 76 6.0 77 2.0 78 5.0 79 9.0 80 6.0 81 6.0 82 3.0 83 4.0 84 5.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 3928.0 >>END_MODULE >>Sequence Duplication Levels pass #Total Deduplicated Percentage 84.275 #Duplication Level Percentage of deduplicated Percentage of total 1 94.09670720854346 79.3 2 3.8564224265796496 6.5 3 0.5932957579353307 1.5 4 0.5339661821417977 1.7999999999999998 5 0.20765351527736575 0.8750000000000001 6 0.14832393948383268 0.75 7 0.08899436369029962 0.525 8 0.02966478789676654 0.2 9 0.02966478789676654 0.22499999999999998 >10 0.35597745476119846 5.525 >50 0.05932957579353308 2.8000000000000003 >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 57 1.425 No Hit GGGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTGG 55 1.375 No Hit GGGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGT 47 1.175 No Hit GGGGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTG 24 0.6 No Hit GGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTGGG 22 0.5499999999999999 No Hit GGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGTG 21 0.525 No Hit GGAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAA 20 0.5 No Hit GGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTT 18 0.44999999999999996 No Hit GGCCTGTAGTAGGAATCTGGTTCACTGCTTTAGGTATTAGTACTATGGCG 12 0.3 No Hit TACCTAGGCACCCAGAGACGAGGAAGGGCGTAGCAAGCGACGAAATGCTT 12 0.3 No Hit GGATTCAATTTCAACCAATCTGTAGTTGATAGTCAAGGTCGCGTTATTAA 12 0.3 No Hit GGAAAAGAAAGCAAAAGCGATTCCCGTAGTAGCGGCGAGCGAAATGGGAG 12 0.3 No Hit GGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAATCG 11 0.27499999999999997 No Hit GGGGAAGTACACCAGCGACGGCGAGGCCGCCGCCGCCAAGGAAGGCATGT 10 0.25 No Hit GGTGGAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGC 9 0.22499999999999998 No Hit GGAGTATCCTTGATATATAAATATATAGATAAATAGATTTAAATGGAAAA 8 0.2 No Hit GGGAGATGCTAACTCAGATGCCATGAAGACTGGTTCCTTCTACGGTTAGA 7 0.17500000000000002 No Hit GGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGTGTACAAGCTCGTAA 7 0.17500000000000002 No Hit GGGCGCGATCTTGCTCGTGAAGGTAATGAAATTATCCGAGCAGCTTGCAA 7 0.17500000000000002 No Hit GACCACAGCGCATTCCCAATTAGCATCTAAGCTCTCGTTGACATTTCCTT 6 0.15 No Hit GGAGCGGAAGAGGAACGGCGCTGCGACCGCCGCGCCCCAGATGAACATGA 6 0.15 No Hit GGGAGTATTATGAAAGGAGATTAATCATAGATTATCAAAAACCCTAGAAA 6 0.15 No Hit GGCGTTCAACCTAAATGGATTCAATTTCAACCAATCTGTAGTTGATAGTC 6 0.15 No Hit GGGCTGATATCATCAACCGTGCTAATCTTGGTATGGAAGTAATGCACGAA 6 0.15 No Hit GTAGTAGGAATCTGGTTCACTGCTTTAGGTATTAGTACTATGGCGTTCAA 5 0.125 No Hit GAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAAT 5 0.125 No Hit GGTCGCGTTATTAATACTTGGGCTGATATCATCAACCGTGCTAATCTTGG 5 0.125 No Hit GGAGAGCAATACAAGCGTTGTGCTGCTAGGCGAAGCGGTTGAGTGCCGCA 5 0.125 No Hit GGGGGTCGCAGTGACCAGGCCCGGGCGACTGTTTACCAAAAACACAGGTC 5 0.125 No Hit GGAGTTGAAAATAAGCATAGATCCGGAGATTCCCAAATAGGTCAACCTTT 5 0.125 No Hit GAAATCTTTGGGGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGG 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.025 0.0 0.0 0.0 0.0 16-17 0.025 0.0 0.0 0.0 0.0 18-19 0.025 0.0 0.0 0.0 0.0 20-21 0.025 0.0 0.0 0.0 0.0 22-23 0.025 0.0 0.0 0.0 0.0 24-25 0.025 0.0 0.0 0.0 0.0 26-27 0.025 0.0 0.0 0.0 0.0 28-29 0.025 0.0 0.0 0.0 0.0 30-31 0.025 0.0 0.0 0.0 0.0 32-33 0.025 0.0 0.0 0.0 0.0 34-35 0.025 0.0 0.0 0.0 0.0 36-37 0.025 0.0 0.0 0.0 0.0 38-39 0.025 0.0 0.0 0.0 0.0 40-41 0.025 0.0 0.0 0.0 0.0 42-43 0.025 0.0 0.0 0.0 0.0 44-45 0.025 0.0 0.0 0.0 0.0 46-47 0.025 0.0 0.0 0.0 0.0 48-49 0.025 0.0 0.0 0.0 0.0 50-51 0.025 0.0 0.0 0.0 0.0 52-53 0.025 0.0 0.0 0.0 0.0 54-55 0.025 0.0 0.0 0.0 0.0 56-57 0.025 0.0 0.0 0.0 0.0 58-59 0.025 0.0 0.0 0.0 0.0 60-61 0.025 0.0 0.0 0.0 0.0 62-63 0.025 0.0 0.0 0.0 0.0 64-65 0.025 0.0 0.0 0.0 0.0 66-67 0.025 0.0 0.0 0.0 0.0 68-69 0.025 0.0 0.0 0.0 0.0 70-71 0.025 0.0 0.0 0.0 0.0 72-73 0.025 0.0 0.0 0.0 0.0 74-75 0.025 0.0 0.0 0.0 0.0 76-77 0.025 0.0 0.0 0.0 0.0 78-79 0.025 0.0 0.0 0.0 0.0 80-81 0.025 0.0 0.0 0.0 0.0 >>END_MODULE >>Kmer Content warn #Sequence Count PValue Obs/Exp Max Max Obs/Exp Position ATCTTGC 20 7.824886E-4 43.5 76-77 GGCGCGA 20 7.824886E-4 43.5 70-71 CTCGTGA 20 7.824886E-4 43.5 82-83 CGATCTT 20 7.824886E-4 43.5 74-75 CGCGATC 20 7.824886E-4 43.5 72-73 AGGGCGC 25 0.0023441382 34.8 68-69 TGCTCGT 25 0.0023441382 34.8 80-81 CTTGCTC 25 0.0023441382 34.8 78-79 TAGAAGC 30 0.0057259556 29.0 40-41 AGCTCGT 30 0.0057259556 29.0 56-57 CCTGTGT 30 0.0057259556 29.0 46-47 GAAGCCT 30 0.0057259556 29.0 42-43 CGTAACG 30 0.0057259556 29.0 60-61 CGTGAAG 30 0.0057259556 29.0 84-85 TGTACAA 30 0.0057259556 29.0 50-51 ACGAAGG 30 0.0057259556 29.0 64-65 CTCGTAA 30 0.0057259556 29.0 58-59 >>END_MODULE Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125997 READS because READLEN < 1 Read 125997 spots for ERR6133556.sra Written 125997 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra Rejected 125986 READS because READLEN < 1 Read 125986 spots for ERR6133556.sra Written 125986 spots for ERR6133556.sra SRR ids: ['ERR6133556.sra'] extra args: ['--split-files', '--defline-qual', '+'] tempdir: /tmp/pfd_4i418mdn ERR6133556.sra spots: 2519731 blocks: [[1, 125986], [125987, 251972], [251973, 377958], [377959, 503944], [503945, 629930], [629931, 755916], [755917, 881902], [881903, 1007888], [1007889, 1133874], [1133875, 1259860], [1259861, 1385846], [1385847, 1511832], [1511833, 1637818], [1637819, 1763804], [1763805, 1889790], [1889791, 2015776], [2015777, 2141762], [2141763, 2267748], [2267749, 2393734], [2393735, 2519731]] ERR6133556 file size 557089 ERR6133556 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 ERR6133556 ERR6133556_1.fastq Input file: ERR6133556_1.fastq trimmed: ERR6133556-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:33:43 2024 >> started Sat Dec 7 08:33:45 2024 >> done (1.513s) 2519731 reads processed; of these: 139 ( 0.01%) short reads filtered out after trimming by size control 33 ( 0.00%) empty reads filtered out after trimming by size control 2519559 (99.99%) reads available; of these: 60351 ( 2.40%) trimmed reads available after processing 2459208 (97.60%) untrimmed reads available after processing Length distribution of reads after trimming: length count percentage 18 19 0.00% 19 19 0.00% 20 16 0.00% 21 13 0.00% 22 20 0.00% 23 10 0.00% 24 8 0.00% 25 12 0.00% 26 11 0.00% 27 10 0.00% 28 26 0.00% 29 27 0.00% 30 10 0.00% 31 13 0.00% 32 14 0.00% 33 18 0.00% 34 17 0.00% 35 98 0.00% 36 520 0.02% 37 25 0.00% 38 19 0.00% 39 63 0.00% 40 53 0.00% 41 32 0.00% 42 15 0.00% 43 11 0.00% 44 14 0.00% 45 11 0.00% 46 7 0.00% 47 12 0.00% 48 5 0.00% 49 6 0.00% 50 3 0.00% 51 18 0.00% 52 4 0.00% 53 5 0.00% 54 6 0.00% 55 5 0.00% 56 6 0.00% 57 3 0.00% 58 8 0.00% 59 8 0.00% 60 17 0.00% 61 12 0.00% 62 0 0.00% 63 1 0.00% 64 2 0.00% 65 5 0.00% 66 8 0.00% 67 9 0.00% 68 13 0.00% 69 46 0.00% 70 4043 0.16% 71 3614 0.14% 72 3777 0.15% 73 3649 0.14% 74 3717 0.15% 75 3553 0.14% 76 3418 0.14% 77 3476 0.14% 78 3753 0.15% 79 4095 0.16% 80 3770 0.15% 81 3934 0.16% 82 4210 0.17% 83 4456 0.18% 84 4148 0.16% 85 175 0.01% 86 280 0.01% 87 380 0.02% 88 765 0.03% 89 1406 0.06% 90 2870 0.11% 91 8529 0.34% 92 43604 1.73% 93 2402604 95.36% 2519559 reads passed initial QC criterion=sequence-density sequence-density=0.50 sequence-density-rank=1 fanout-score=5.63 fanout-score-rank=21 prefix-density=1.14 prefix-fanout=2.5 sequence=GCCGCCGCCGCCAAGGAAGGCATGTTCGTCAAGAACTACAGCTACTGATCCTAATCGCATCAAGCTTCAACGCCTGTGAGTGAAAACCAGTGATGAGAGTGCTGCTGCTAGCTAGCGCCGGCATTGATGAGCTTGAGAGGGCACTGTAGCCAGTGTGTCAGTCGTTGTTAAATTACAGGTTGAGATCATCAGCGTACTCCGATGGGAGGTGGACATCAGAAAGTATACTGTGTTTTACCACCCTAATTAAGTAAACAACTTTTGGAATGGTGATCATCTTACTGTTTTAAATAAATCTGTTTC criterion=fanout-score sequence-density=0.03 sequence-density-rank=21 fanout-score=44.09 fanout-score-rank=1 prefix-density=0.18 prefix-fanout=8.2 sequence=AAAAGAAGGGGTGTTCCATCTCCGGACGACGATCCTGCCTGCAGAGGAAGACATGCCGGCGATCATGTCGAGCTTCAAGAAGTTCAACGACTCATTCATGGAGCAATACCAAGACTACTCCAGGCTGTGATGTGAAGAGGGAAACAACGAGGTCATCATCGACATGATATATTGCTGCTATTTTCCACCAGCGATTAAAAGTTAAAAAATTTAGCTGTAAGCTGTAACTATCTTGAAGAAACT Started job on | Dec 07 08:33:59 Started mapping on | Dec 07 08:34:00 Finished on | Dec 07 08:34:05 Mapping speed, Million of reads per hour | 1814.08 Number of input reads | 2519559 Average input read length | 92 UNIQUE READS: Uniquely mapped reads number | 1877536 Uniquely mapped reads % | 74.52% Average mapped length | 92.19 Number of splices: Total | 154772 Number of splices: Annotated (sjdb) | 132648 Number of splices: GT/AG | 150580 Number of splices: GC/AG | 2195 Number of splices: AT/AC | 51 Number of splices: Non-canonical | 1946 Mismatch rate per base, % | 0.51% Deletion rate per base | 0.05% Deletion average length | 2.00 Insertion rate per base | 0.02% Insertion average length | 2.12 MULTI-MAPPING READS: Number of reads mapped to multiple loci | 559053 % of reads mapped to multiple loci | 22.19% Number of reads mapped to too many loci | 15059 % of reads mapped to too many loci | 0.60% UNMAPPED READS: % of reads unmapped: too many mismatches | 0.00% % of reads unmapped: too short | 2.66% % of reads unmapped: other | 0.04% CHIMERIC READS: Number of chimeric reads | 0 % of chimeric reads | 0.00% N_unmapped 82970 82970 82970 N_multimapping 559053 559053 559053 N_noFeature 88358 108175 1799132 N_ambiguous 66780 8301 222 UnstrandedReadsAssigned:1722398 PositiveStrandReadsAssigned:1761060 NegativeStrandReadsAssigned:78182 Dataset is classified positive stranded MeadianReadLen=93 20thPercentileLength=93 echo kmer=89 ERR6133556 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: ERR6133556-trimmed.fastq [quant] finding pseudoalignments for the reads ... done [quant] processed 2,519,559 reads, 2,160,570 reads pseudoaligned [ em] quantifying the abundances ... done [ em] the Expectation-Maximization algorithm ran for 1,078 rounds 52973 ERR6133556.ke.tsv 35125 ERR6133556.se.tsv 88098 total ==> ERR6133556.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 78 34.7328 PNS24243 293 194 0 0 KQK14069 1603 1504 321.962 130.784 KQK14071 474 375 1.03831 1.69158 ==> ERR6133556.se.tsv <== BRADI_1g14170v3 322 BRADI_1g53295v3 12 BRADI_1g59795v3 16 BRADI_1g07683v3 0 BRADI_1g00485v3 0 BRADI_1g20270v3 29 BRADI_1g74790v3 13 BRADI_1g09890v3 0 BRADI_1g77505v3 22 BRADI_1g48960v3 0 ERR6133556 completed mapping pipeline successfully