Starting /dee2/code/volunteer_pipeline.sh ERR6133451 current disk space = 1515896098816 free memory = 1597865080 ERR6133451 SRAfilesize 11d446e9fc23be2312422176ab3b31f2 ERR6133451.sra ERR6133451.sra file validated ERR6133451 is single end ERR6133451 is conventional basespace ERR6133451 read1 length is 70-93 nt ##FastQC 0.11.5 >>Basic Statistics pass #Measure Value Filename ERR6133451_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 44 >>END_MODULE >>Per base sequence quality pass #Base Mean Median Lower Quartile Upper Quartile 10th Percentile 90th Percentile 1 34.74975 37.0 33.0 37.0 33.0 37.0 2 36.0205 37.0 37.0 37.0 33.0 37.0 3 35.359 37.0 33.0 37.0 33.0 37.0 4 34.98725 37.0 37.0 37.0 33.0 37.0 5 34.93975 37.0 37.0 37.0 27.0 37.0 6 35.33 37.0 37.0 37.0 33.0 37.0 7 36.84325 37.0 37.0 40.0 33.0 40.0 8 37.0775 37.0 37.0 40.0 33.0 40.0 9 37.09125 37.0 37.0 40.0 33.0 40.0 10-11 37.008750000000006 37.0 37.0 40.0 33.0 40.0 12-13 36.93775 37.0 37.0 40.0 33.0 40.0 14-15 36.990750000000006 37.0 37.0 40.0 33.0 40.0 16-17 36.768 37.0 37.0 40.0 33.0 40.0 18-19 36.397000000000006 37.0 37.0 40.0 33.0 40.0 20-21 36.35475 37.0 37.0 40.0 33.0 40.0 22-23 35.991875 37.0 35.0 40.0 33.0 40.0 24-25 36.20625 37.0 37.0 40.0 33.0 40.0 26-27 36.35275 37.0 37.0 40.0 33.0 40.0 28-29 36.391999999999996 37.0 37.0 40.0 33.0 40.0 30-31 36.476875 37.0 37.0 40.0 33.0 40.0 32-33 36.33675 37.0 37.0 40.0 33.0 40.0 34-35 36.331125 37.0 37.0 40.0 33.0 40.0 36-37 36.30625 37.0 37.0 40.0 33.0 40.0 38-39 36.244125 37.0 37.0 40.0 33.0 40.0 40-41 36.168375 37.0 37.0 40.0 33.0 40.0 42-43 36.066 37.0 37.0 40.0 33.0 40.0 44-45 35.77825 37.0 35.0 40.0 33.0 40.0 46-47 35.73875 37.0 33.0 40.0 33.0 40.0 48-49 35.870374999999996 37.0 33.0 37.0 33.0 40.0 50-51 35.685500000000005 37.0 33.0 37.0 33.0 40.0 52-53 35.3635 37.0 33.0 37.0 33.0 40.0 54-55 35.355625 37.0 33.0 37.0 33.0 40.0 56-57 35.198625 37.0 33.0 37.0 33.0 40.0 58-59 34.654875000000004 37.0 33.0 37.0 27.0 38.5 60-61 34.664125 37.0 33.0 37.0 27.0 37.0 62-63 34.516625000000005 37.0 33.0 37.0 27.0 37.0 64-65 34.497125 37.0 33.0 37.0 27.0 37.0 66-67 34.235 37.0 33.0 37.0 27.0 37.0 68-69 33.422875000000005 35.0 33.0 37.0 27.0 37.0 70-71 33.60044974874372 35.0 33.0 37.0 27.0 37.0 72-73 33.96549976412175 37.0 33.0 37.0 27.0 37.0 74-75 33.987560045204845 37.0 33.0 37.0 27.0 37.0 76-77 33.89382203424153 37.0 33.0 37.0 27.0 37.0 78-79 33.89827257005854 37.0 33.0 37.0 27.0 37.0 80-81 33.79514228732691 37.0 33.0 37.0 27.0 37.0 82-83 33.60642458125284 37.0 33.0 37.0 27.0 37.0 84-85 33.50803832911147 35.0 33.0 37.0 27.0 37.0 86-87 33.276162790697676 35.0 33.0 37.0 27.0 37.0 88-89 33.37721483942414 37.0 33.0 37.0 27.0 37.0 90-91 33.244739756367665 33.0 33.0 37.0 27.0 37.0 92-93 33.16348283499447 33.0 33.0 37.0 27.0 37.0 >>END_MODULE >>Per sequence quality scores pass #Quality Count 20 7.0 21 13.0 22 12.0 23 27.0 24 20.0 25 28.0 26 40.0 27 61.0 28 54.0 29 86.0 30 108.0 31 131.0 32 147.0 33 194.0 34 274.0 35 521.0 36 794.0 37 872.0 38 569.0 39 42.0 >>END_MODULE >>Per base sequence content fail #Base G A T C 1 81.72500000000001 7.95 2.5250000000000004 7.8 2 66.85 18.875 9.049999999999999 5.225 3 37.2 37.0 14.7 11.1 4 34.175 27.500000000000004 19.85 18.475 5 24.325 31.574999999999996 24.099999999999998 20.0 6 19.400000000000002 37.8 26.700000000000003 16.1 7 34.175 30.0 20.8 15.024999999999999 8 30.725 28.999999999999996 24.775 15.5 9 26.525 26.650000000000002 28.525 18.3 10-11 25.575 27.8625 27.987499999999997 18.575 12-13 28.5625 25.525 29.049999999999997 16.8625 14-15 22.4625 27.275 31.324999999999996 18.9375 16-17 24.075 30.5375 27.400000000000002 17.9875 18-19 23.3375 27.125 28.199999999999996 21.337500000000002 20-21 24.487243621810904 26.23811905952976 28.026513256628316 21.248124062031014 22-23 25.637500000000003 25.687500000000004 28.3625 20.3125 24-25 25.924999999999997 25.6 28.3125 20.1625 26-27 23.9 26.174999999999997 30.4375 19.4875 28-29 24.575 27.05 28.7375 19.6375 30-31 25.45 26.3125 28.037499999999998 20.200000000000003 32-33 23.5125 28.0875 28.812500000000004 19.5875 34-35 25.3 24.85 28.8625 20.9875 36-37 24.8 24.2625 29.6625 21.275 38-39 23.775 25.650000000000002 31.55 19.025 40-41 26.1125 25.624999999999996 28.075 20.1875 42-43 23.76720901126408 27.63454317897372 28.16020025031289 20.438047559449313 44-45 23.7125 26.0125 29.7375 20.5375 46-47 24.425 24.349999999999998 28.749999999999996 22.475 48-49 23.5 26.650000000000002 31.1875 18.6625 50-51 24.175 27.575 29.925 18.325 52-53 24.611723446893787 27.89328657314629 28.269038076152302 19.225951903807616 54-55 23.865483185398176 29.041130141267658 28.091011376422053 19.002375296912113 56-57 24.65 26.400000000000002 29.612500000000004 19.3375 58-59 22.975 26.187500000000004 30.2375 20.599999999999998 60-61 23.799999999999997 25.124999999999996 31.225 19.85 62-63 21.087500000000002 29.975 30.7875 18.15 64-65 23.9 26.525 29.9375 19.6375 66-67 23.8875 28.487499999999997 28.799999999999997 18.825 68-69 22.277784723090384 27.29091136392049 30.24128016002 20.19002375296912 70-71 25.36340852130326 25.062656641604008 29.736842105263158 19.83709273182957 72-73 24.844167408726626 25.54382394097443 29.271085103676377 20.340923546622566 74-75 22.62150982419855 26.719234746639092 30.765253360910034 19.894002068252327 76-77 22.473821989528798 26.034031413612563 31.544502617801047 19.94764397905759 78-79 23.236734153764722 24.970226280269948 31.89096202196639 19.902077543998942 80-81 24.099611728477708 26.911233096800107 30.94122372472888 18.047931449993307 82-83 23.737099402498643 24.74198804997284 31.6811515480717 19.839760999456818 84-85 23.78124568429775 24.209363347603922 31.01781521889242 20.99157574920591 86-87 23.034330011074196 26.6749723145072 32.43355481727575 17.857142857142858 88-89 21.566998892580287 29.512735326688816 31.367663344407532 17.552602436323365 90-91 24.82004429678848 27.380952380952383 29.2358803986711 18.56312292358804 92-93 22.093023255813954 29.983388704318937 29.651162790697676 18.272425249169437 >>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 5.0 18 6.5 19 1.5 20 1.0 21 2.0 22 4.0 23 7.0 24 8.0 25 10.5 26 13.5 27 21.5 28 36.0 29 45.0 30 52.0 31 58.5 32 67.0 33 87.0 34 93.0 35 113.0 36 132.0 37 143.0 38 183.0 39 195.5 40 191.0 41 212.0 42 235.0 43 240.5 44 213.0 45 174.5 46 181.0 47 185.0 48 171.0 49 175.5 50 158.0 51 134.0 52 130.5 53 136.0 54 116.0 55 72.0 56 56.0 57 57.0 58 49.5 59 37.5 60 27.5 61 22.5 62 24.5 63 28.0 64 29.5 65 16.5 66 9.0 67 14.5 68 12.0 69 7.0 70 5.0 71 5.0 72 5.0 73 4.0 74 2.5 75 2.0 76 1.0 77 0.0 78 0.0 79 0.5 80 0.5 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.05 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.0 40-41 0.0 42-43 0.125 44-45 0.0 46-47 0.0 48-49 0.0 50-51 0.0 52-53 0.2 54-55 0.0125 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.0125 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 20.0 71 36.0 72 27.0 73 35.0 74 28.0 75 23.0 76 22.0 77 22.0 78 17.0 79 24.0 80 23.0 81 25.0 82 32.0 83 37.0 84 17.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 3612.0 >>END_MODULE >>Sequence Duplication Levels pass #Total Deduplicated Percentage 86.825 #Duplication Level Percentage of deduplicated Percentage of total 1 94.27008350129572 81.85 2 3.4840195796141664 6.05 3 0.9501871580765908 2.475 4 0.48949035416066805 1.7000000000000002 5 0.2303484019579614 1.0 6 0.14396775122372588 0.75 7 0.05758710048949035 0.35000000000000003 8 0.028793550244745177 0.2 9 0.0 0.0 >10 0.34552260293694215 5.625 >50 0.0 0.0 >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 warn #Sequence Count Percentage Possible Source GGGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGT 39 0.975 No Hit GGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGTG 32 0.8 No Hit GGGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTGG 32 0.8 No Hit GGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTT 18 0.44999999999999996 No Hit GGATTCAATTTCAACCAATCTGTAGTTGATAGTCAAGGTCGCGTTATTAA 17 0.42500000000000004 No Hit GGGCGCGATCTTGCTCGTGAAGGTAATGAAATTATCCGAGCAGCTTGCAA 16 0.4 No Hit GGGGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTG 15 0.375 No Hit GGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGTGTACAAGCTCGTAA 14 0.35000000000000003 No Hit GGAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAA 12 0.3 No Hit GGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAATCG 10 0.25 No Hit GGAGTATCCTTGATATATAAATATATAGATAAATAGATTTAAATGGAAAA 10 0.25 No Hit GGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTGGG 10 0.25 No Hit GGGGATGGAGCGACAGAAGTATGGAAATAGGATAAGGTAGCGGCGAGACG 8 0.2 No Hit GGGCTTTTTTAGTATTTATCTAAAGGAAGGAACAAACGAGGATAAGATAA 7 0.17500000000000002 No Hit GGCGTTCAACCTAAATGGATTCAATTTCAACCAATCTGTAGTTGATAGTC 7 0.17500000000000002 No Hit GGGAGTATTATGAAAGGAGATTAATCATAGATTATCAAAAACCCTAGAAA 6 0.15 No Hit GGGAAACAACGAGGTCATCATCGACATGATATATTGCTGCTATTTTCCAC 6 0.15 No Hit GGTAATGAAATTATCCGAGCAGCTTGCAAATGGAGTCCTGAACTAGCCGC 6 0.15 No Hit GGGAATCGATCGTGATTTAGAACCTGTTCTTTACATGAACCCTCTTAACT 6 0.15 No Hit GGGGGCATTCGTATTTCATAGTCAGAGGTGAAATTCTTGGATTTATGAAA 6 0.15 No Hit AACGTATGAATTTGTGTGATGTGATGTACCCGTGCTATTCTCAGTGAAAT 5 0.125 No Hit GGTGGAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGC 5 0.125 No Hit GGCGCGATCTTGCTCGTGAAGGTAATGAAATTATCCGAGCAGCTTGCAAA 5 0.125 No Hit GGAAGAGCCCGAGCTGCTGCAGCAGGTTTTGAAAAGGGAATCGATCGTGA 5 0.125 No Hit CAACCTAAATGGATTCAATTTCAACCAATCTGTAGTTGATAGTCAAGGTC 5 0.125 No Hit GAAGGTAATGAAATTATCCGAGCAGCTTGCAAATGGAGTCCTGAACTAGC 5 0.125 No Hit CAGCTAATCGAGTGGCTTTAGAAGCCTGTGTACAAGCTCGTAACGAAGGG 5 0.125 No Hit GTATTTAGCCTTGCAAGGTGGTCCTTGCTGATTCACACGGGATTCCACGT 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.0125 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 GGGAAAT 25 0.0071765394 51.12 1 TAGAAGC 20 8.67354E-4 42.6 38-39 GAAGCCT 20 8.67354E-4 42.6 40-41 TGTACAA 20 8.67354E-4 42.6 48-49 AGGGCGC 25 0.0025177265 34.29434 66-67 AGCCTGT 25 0.002597359 34.079998 42-43 AGCTCGT 25 0.002597359 34.079998 54-55 CCTGTGT 25 0.002597359 34.079998 44-45 TACAAGC 25 0.002597359 34.079998 50-51 CGTAACG 25 0.002597359 34.079998 58-59 ACGAAGG 25 0.002597359 34.079998 62-63 CTCGTAA 25 0.002597359 34.079998 56-57 CAAGCTC 25 0.002597359 34.079998 52-53 TGAAGGT 30 0.0030450604 32.927536 84-85 AAGGTAA 30 0.0030450604 32.927536 86-87 CGTGAAG 30 0.0030450604 32.927536 82-83 CTCGTGA 30 0.0030450604 32.927536 80-81 TGCTCGT 30 0.0031563714 32.690647 78-79 GGCGCGA 30 0.0059593343 28.759495 68-69 CGCGATC 30 0.0059593343 28.759495 70-71 >>END_MODULE Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra Rejected 91489 READS because READLEN < 1 Read 91489 spots for ERR6133451.sra Written 91489 spots for ERR6133451.sra SRR ids: ['ERR6133451.sra'] extra args: ['--split-files', '--defline-qual', '+'] tempdir: /tmp/pfd_s2ckte4u ERR6133451.sra spots: 1829780 blocks: [[1, 91489], [91490, 182978], [182979, 274467], [274468, 365956], [365957, 457445], [457446, 548934], [548935, 640423], [640424, 731912], [731913, 823401], [823402, 914890], [914891, 1006379], [1006380, 1097868], [1097869, 1189357], [1189358, 1280846], [1280847, 1372335], [1372336, 1463824], [1463825, 1555313], [1555314, 1646802], [1646803, 1738291], [1738292, 1829780]] ERR6133451 file size 399438 ERR6133451 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 ERR6133451 ERR6133451_1.fastq Input file: ERR6133451_1.fastq trimmed: ERR6133451-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 Thu Dec 12 02:17:16 2024 >> started Thu Dec 12 02:17:17 2024 >> done (1.011s) 1829780 reads processed; of these: 153 ( 0.01%) short reads filtered out after trimming by size control 9 ( 0.00%) empty reads filtered out after trimming by size control 1829618 (99.99%) reads available; of these: 25468 ( 1.39%) trimmed reads available after processing 1804150 (98.61%) untrimmed reads available after processing Length distribution of reads after trimming: length count percentage 18 16 0.00% 19 30 0.00% 20 19 0.00% 21 15 0.00% 22 17 0.00% 23 19 0.00% 24 10 0.00% 25 8 0.00% 26 10 0.00% 27 13 0.00% 28 103 0.01% 29 54 0.00% 30 18 0.00% 31 20 0.00% 32 23 0.00% 33 25 0.00% 34 18 0.00% 35 44 0.00% 36 307 0.02% 37 18 0.00% 38 21 0.00% 39 42 0.00% 40 43 0.00% 41 41 0.00% 42 16 0.00% 43 14 0.00% 44 23 0.00% 45 10 0.00% 46 23 0.00% 47 11 0.00% 48 7 0.00% 49 4 0.00% 50 16 0.00% 51 37 0.00% 52 9 0.00% 53 5 0.00% 54 12 0.00% 55 5 0.00% 56 27 0.00% 57 36 0.00% 58 4 0.00% 59 6 0.00% 60 14 0.00% 61 5 0.00% 62 2 0.00% 63 4 0.00% 64 5 0.00% 65 2 0.00% 66 8 0.00% 67 14 0.00% 68 34 0.00% 69 132 0.01% 70 14106 0.77% 71 14273 0.78% 72 12824 0.70% 73 12764 0.70% 74 12631 0.69% 75 12684 0.69% 76 11187 0.61% 77 10825 0.59% 78 12248 0.67% 79 11301 0.62% 80 11061 0.60% 81 10838 0.59% 82 11861 0.65% 83 13957 0.76% 84 9858 0.54% 85 61 0.00% 86 100 0.01% 87 150 0.01% 88 289 0.02% 89 631 0.03% 90 1255 0.07% 91 4009 0.22% 92 15930 0.87% 93 1623356 88.73% 1829618 reads passed initial QC criterion=sequence-density sequence-density=0.39 sequence-density-rank=1 fanout-score=3.78 fanout-score-rank=27 prefix-density=1.42 prefix-fanout=1.0 sequence=ATTCGAGTTCGAGCCGGTAGATAAACTAGATAGCTAGACTAAGTGGATAAAATTAGATAGAAAAAAGGTCTAAA criterion=fanout-score sequence-density=0.02 sequence-density-rank=27 fanout-score=283.54 fanout-score-rank=1 prefix-density=0.79 prefix-fanout=6.2 sequence=GAAGAAGAAAAGTTTTCTCAACATGGGGAGGAAGTCCCTCCGAAATTTGATTTGTTATTGTATTGTAAGGGGCTTTTTTAGTATTTATCTAAAGGAAGGAACAAACGAGGATAAGATAAAATTTCTTTAAATTTATTTTGCCCAAGTGGGATATATGGCATACTATTCTTTCCATTTTCATCTTTTTTTCTATTCCACTCCATCTAGATATAAGAAAGAACCCAATGCAATGAAATTCCACTAATATACAATACAAAAAAGAAGAA Started job on | Dec 12 02:17:35 Started mapping on | Dec 12 02:17:35 Finished on | Dec 12 02:17:44 Mapping speed, Million of reads per hour | 731.85 Number of input reads | 1829618 Average input read length | 91 UNIQUE READS: Uniquely mapped reads number | 1238170 Uniquely mapped reads % | 67.67% Average mapped length | 91.93 Number of splices: Total | 48199 Number of splices: Annotated (sjdb) | 38523 Number of splices: GT/AG | 45746 Number of splices: GC/AG | 1757 Number of splices: AT/AC | 15 Number of splices: Non-canonical | 681 Mismatch rate per base, % | 0.48% Deletion rate per base | 0.04% Deletion average length | 1.85 Insertion rate per base | 0.02% Insertion average length | 1.73 MULTI-MAPPING READS: Number of reads mapped to multiple loci | 332375 % of reads mapped to multiple loci | 18.17% Number of reads mapped to too many loci | 16949 % of reads mapped to too many loci | 0.93% UNMAPPED READS: % of reads unmapped: too many mismatches | 0.00% % of reads unmapped: too short | 13.12% % of reads unmapped: other | 0.12% CHIMERIC READS: Number of chimeric reads | 0 % of chimeric reads | 0.00% N_unmapped 259073 259073 259073 N_multimapping 332375 332375 332375 N_noFeature 76133 87062 1183417 N_ambiguous 48762 4947 141 UnstrandedReadsAssigned:1113275 PositiveStrandReadsAssigned:1146161 NegativeStrandReadsAssigned:54612 Dataset is classified positive stranded MeadianReadLen=93 20thPercentileLength=93 echo kmer=89 ERR6133451 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: ERR6133451-trimmed.fastq [quant] finding pseudoalignments for the reads ... done [quant] processed 1,829,618 reads, 1,372,934 reads pseudoaligned [ em] quantifying the abundances ... done [ em] the Expectation-Maximization algorithm ran for 1,033 rounds 52973 ERR6133451.ke.tsv 35125 ERR6133451.se.tsv 88098 total ==> ERR6133451.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 46 32.8096 PNS24243 293 194 0 0 KQK14069 1603 1504 11 7.15717 KQK14071 474 375 0 0 ==> ERR6133451.se.tsv <== BRADI_1g14170v3 11 BRADI_1g53295v3 19 BRADI_1g59795v3 16 BRADI_1g07683v3 0 BRADI_1g00485v3 0 BRADI_1g20270v3 22 BRADI_1g74790v3 11 BRADI_1g09890v3 0 BRADI_1g77505v3 33 BRADI_1g48960v3 0 ERR6133451 completed mapping pipeline successfully