Starting /dee2/code/volunteer_pipeline.sh ERR6133426 current disk space = 1546870525952 free memory = 1432964132 ERR6133426 SRAfilesize cd6713711090e033575e1071454d6457 ERR6133426.sra ERR6133426.sra file validated ERR6133426 is single end ERR6133426 is conventional basespace ERR6133426 read1 length is 70-93 nt ##FastQC 0.11.5 >>Basic Statistics pass #Measure Value Filename ERR6133426_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 46 >>END_MODULE >>Per base sequence quality pass #Base Mean Median Lower Quartile Upper Quartile 10th Percentile 90th Percentile 1 35.25925 37.0 33.0 37.0 33.0 37.0 2 36.40625 37.0 37.0 37.0 33.0 37.0 3 35.7785 37.0 37.0 37.0 33.0 37.0 4 35.508 37.0 37.0 37.0 33.0 37.0 5 35.2345 37.0 37.0 37.0 33.0 37.0 6 35.5775 37.0 37.0 37.0 33.0 37.0 7 37.2295 37.0 37.0 40.0 33.0 40.0 8 37.335 37.0 37.0 40.0 33.0 40.0 9 37.3335 37.0 37.0 40.0 33.0 40.0 10-11 37.333749999999995 37.0 37.0 40.0 33.0 40.0 12-13 37.230625 37.0 37.0 40.0 33.0 40.0 14-15 37.066 37.0 37.0 40.0 33.0 40.0 16-17 36.887375 37.0 37.0 40.0 33.0 40.0 18-19 37.068124999999995 37.0 37.0 40.0 33.0 40.0 20-21 37.343 37.0 37.0 40.0 33.0 40.0 22-23 37.047125 37.0 37.0 40.0 33.0 40.0 24-25 37.148875000000004 37.0 37.0 40.0 33.0 40.0 26-27 37.069874999999996 37.0 37.0 40.0 33.0 40.0 28-29 37.121 37.0 37.0 40.0 33.0 40.0 30-31 37.17225 37.0 37.0 40.0 33.0 40.0 32-33 37.01375 37.0 37.0 40.0 33.0 40.0 34-35 36.812 37.0 37.0 40.0 33.0 40.0 36-37 36.68825 37.0 37.0 40.0 33.0 40.0 38-39 36.578625 37.0 37.0 40.0 33.0 40.0 40-41 36.481125 37.0 37.0 40.0 33.0 40.0 42-43 36.441 37.0 37.0 40.0 33.0 40.0 44-45 36.177375 37.0 37.0 38.5 33.0 40.0 46-47 36.06975 37.0 37.0 37.0 33.0 40.0 48-49 35.753874999999994 37.0 35.0 37.0 33.0 40.0 50-51 35.670125 37.0 33.0 37.0 33.0 40.0 52-53 35.6175 37.0 33.0 37.0 33.0 40.0 54-55 35.298 37.0 33.0 37.0 33.0 38.5 56-57 35.252250000000004 37.0 33.0 37.0 33.0 37.0 58-59 34.853750000000005 37.0 33.0 37.0 30.0 37.0 60-61 34.6455 37.0 33.0 37.0 33.0 37.0 62-63 34.756375 37.0 33.0 37.0 33.0 37.0 64-65 34.3535 37.0 33.0 37.0 27.0 37.0 66-67 34.317875 37.0 33.0 37.0 27.0 37.0 68-69 33.045375 35.0 33.0 35.0 30.0 37.0 70-71 33.203706217162875 33.0 33.0 37.0 27.0 37.0 72-73 33.97836769646268 37.0 33.0 37.0 27.0 37.0 74-75 34.179627738637606 37.0 33.0 37.0 33.0 37.0 76-77 34.12308070748114 37.0 33.0 37.0 30.0 37.0 78-79 34.23085546392177 37.0 33.0 37.0 30.0 37.0 80-81 34.17151185068804 37.0 33.0 37.0 33.0 37.0 82-83 34.0805532552243 37.0 33.0 37.0 33.0 37.0 84-85 33.980097747419705 37.0 33.0 37.0 27.0 37.0 86-87 33.760070939954396 37.0 33.0 37.0 27.0 37.0 88-89 33.85583987838865 37.0 33.0 37.0 27.0 37.0 90-91 33.672789460349634 35.0 33.0 37.0 27.0 37.0 92-93 33.43286040030402 33.0 33.0 37.0 27.0 37.0 >>END_MODULE >>Per sequence quality scores pass #Quality Count 20 8.0 21 7.0 22 11.0 23 14.0 24 20.0 25 17.0 26 43.0 27 37.0 28 35.0 29 61.0 30 81.0 31 103.0 32 156.0 33 179.0 34 236.0 35 511.0 36 963.0 37 1118.0 38 398.0 39 2.0 >>END_MODULE >>Per base sequence content fail #Base G A T C 1 89.9 2.125 2.325 5.65 2 76.0440110027507 14.678669667416855 5.601400350087522 3.675918979744936 3 36.825 40.699999999999996 12.125 10.35 4 36.199999999999996 28.4 15.975 19.425 5 24.625 32.15 23.95 19.275000000000002 6 20.7 40.325 22.900000000000002 16.075 7 36.1 29.15 18.65 16.1 8 30.125 29.725 23.0 17.150000000000002 9 26.724999999999998 27.925 26.674999999999997 18.675 10-11 25.837500000000002 26.900000000000002 26.700000000000003 20.5625 12-13 30.425 24.6 25.387500000000003 19.5875 14-15 23.875 28.749999999999996 28.262500000000003 19.112499999999997 16-17 25.900000000000002 31.0375 25.2625 17.8 18-19 24.7875 27.05 25.75 22.412499999999998 20-21 24.5 27.650000000000002 25.8125 22.037499999999998 22-23 27.3 23.4625 26.3625 22.875 24-25 27.075 23.45 27.6375 21.837500000000002 26-27 25.3 26.5375 29.2375 18.925 28-29 25.775 26.8625 26.187500000000004 21.175 30-31 27.425 25.112499999999997 26.224999999999998 21.2375 32-33 25.025 27.800000000000004 27.2625 19.9125 34-35 26.25 24.625 27.775 21.349999999999998 36-37 26.3 22.925 28.299999999999997 22.475 38-39 25.2875 25.324999999999996 29.9375 19.45 40-41 26.625 24.6 26.900000000000002 21.875 42-43 26.05 27.462500000000002 25.25 21.2375 44-45 24.224999999999998 24.9125 28.212500000000002 22.650000000000002 46-47 25.4375 22.3125 27.775 24.474999999999998 48-49 24.175 24.5125 30.362499999999997 20.95 50-51 26.1 26.025 28.799999999999997 19.075 52-53 25.0125 26.5 27.6625 20.825 54-55 23.3375 28.487499999999997 27.825 20.349999999999998 56-57 26.737499999999997 24.675 28.4 20.1875 58-59 23.9875 23.674999999999997 29.5 22.8375 60-61 24.9875 24.0375 28.8375 22.1375 62-63 23.325000000000003 26.8 31.162499999999998 18.712500000000002 64-65 24.7875 26.224999999999998 28.249999999999996 20.7375 66-67 25.674999999999997 26.525 27.825 19.975 68-69 24.5 26.687499999999996 26.3625 22.45 70-71 26.072277103913965 24.559209703638864 27.84794297861698 21.520570213830187 72-73 26.236693800876644 24.358171571696932 27.614276768941764 21.790857858484657 74-75 23.419468138484696 26.693426994480685 29.315102860010033 20.572002007024587 76-77 23.439266423816104 25.22296193945484 29.054138927270444 22.28363270945861 78-79 23.94508124448923 24.373346769114498 30.343872024184403 21.337699962211865 80-81 24.293286219081274 26.70368500757193 30.098435133770824 18.90459363957597 82-83 24.823053589484328 24.0267947421638 29.651162790697676 21.498988877654195 84-85 23.91386953768208 23.10322989233692 30.614312856238122 22.368587713742873 86-87 24.372941474537622 25.0316696224981 31.808968837091463 18.786420065872814 88-89 22.333417785659996 29.288066886242714 30.30149480618191 18.07702052191538 90-91 27.02052191537877 25.361033696478337 28.75601722827464 18.862427159868254 92-93 22.156067899670635 28.78135292627312 29.110717000253356 19.95186217380289 >>END_MODULE >>Per sequence GC content pass #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 10.5 18 11.5 19 1.0 20 1.0 21 1.5 22 1.0 23 4.0 24 4.5 25 2.5 26 4.0 27 6.0 28 10.5 29 14.0 30 17.0 31 22.0 32 44.5 33 54.5 34 56.5 35 74.5 36 89.5 37 114.0 38 141.5 39 154.0 40 170.5 41 176.0 42 172.5 43 180.0 44 180.5 45 188.0 46 228.5 47 227.5 48 187.0 49 163.0 50 162.0 51 173.5 52 158.5 53 177.5 54 152.5 55 90.0 56 79.5 57 78.0 58 76.5 59 68.0 60 54.5 61 43.5 62 42.5 63 46.0 64 44.0 65 37.5 66 31.5 67 20.0 68 11.5 69 9.0 70 5.5 71 6.0 72 5.5 73 3.5 74 2.5 75 2.0 76 0.5 77 0.5 78 0.5 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.025 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.0 40-41 0.0 42-43 0.0 44-45 0.0 46-47 0.0 48-49 0.0 50-51 0.0 52-53 0.0 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.012664640324214793 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 3.0 71 2.0 72 5.0 73 2.0 74 4.0 75 3.0 76 1.0 77 8.0 78 5.0 79 4.0 80 2.0 81 1.0 82 8.0 83 3.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 3947.0 >>END_MODULE >>Sequence Duplication Levels pass #Total Deduplicated Percentage 79.77499999999999 #Duplication Level Percentage of deduplicated Percentage of total 1 92.38483234095895 73.7 2 4.512691946098402 7.199999999999999 3 1.3788780946411783 3.3000000000000003 4 0.4700720777185835 1.5 5 0.21936696960200563 0.8750000000000001 6 0.21936696960200563 1.05 7 0.0940144155437167 0.525 8 0.06267627702914447 0.4 9 0.031338138514572234 0.22499999999999998 >10 0.5954246317768724 9.375 >50 0.031338138514572234 1.8499999999999999 >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 GGGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGT 74 1.8499999999999999 No Hit GGGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTGG 47 1.175 No Hit GGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGTG 35 0.8750000000000001 No Hit GGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTGGG 27 0.675 No Hit GGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTT 26 0.65 No Hit GGAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAA 26 0.65 No Hit GGATTCAATTTCAACCAATCTGTAGTTGATAGTCAAGGTCGCGTTATTAA 24 0.6 No Hit GGGGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTG 22 0.5499999999999999 No Hit GGAGTATCCTTGATATATAAATATATAGATAAATAGATTTAAATGGAAAA 20 0.5 No Hit GGGCGCGATCTTGCTCGTGAAGGTAATGAAATTATCCGAGCAGCTTGCAA 20 0.5 No Hit GGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGTGTACAAGCTCGTAA 16 0.4 No Hit GGGAATCGATCGTGATTTAGAACCTGTTCTTTACATGAACCCTCTTAACT 15 0.375 No Hit GGGAGAGCAATACAAGCGTTGTGCTGCTAGGCGAAGCGGTTGAGTGCCGC 15 0.375 No Hit GGAGTCCCATATATATATGTATAAGATGCCAGCCCGGTTTCGTCAACCAT 14 0.35000000000000003 No Hit GGCGTTCAACCTAAATGGATTCAATTTCAACCAATCTGTAGTTGATAGTC 14 0.35000000000000003 No Hit GGTCGCGTTATTAATACTTGGGCTGATATCATCAACCGTGCTAATCTTGG 12 0.3 No Hit TACCTAGGCACCCAGAGACGAGGAAGGGCGTAGCAAGCGACGAAATGCTT 11 0.27499999999999997 No Hit CAACCTAAATGGATTCAATTTCAACCAATCTGTAGTTGATAGTCAAGGTC 11 0.27499999999999997 No Hit GGTGGAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGC 10 0.25 No Hit GGCCTGTAGTAGGAATCTGGTTCACTGCTTTAGGTATTAGTACTATGGCG 10 0.25 No Hit GGAAGAGCCCGAGCTGCTGCAGCAGGTTTTGAAAAGGGAATCGATCGTGA 9 0.22499999999999998 No Hit GGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAATCG 8 0.2 No Hit GTAGTAGGAATCTGGTTCACTGCTTTAGGTATTAGTACTATGGCGTTCAA 8 0.2 No Hit GGAAGAGACTTATAATATTGTGGCTGCTCATGGTTATTTTGGCCGATTAA 7 0.17500000000000002 No Hit GGGATGGAGCGACAGAAGTATGGAAATAGGATAAGGTAGCGGCGAGACGA 7 0.17500000000000002 No Hit GAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAAT 7 0.17500000000000002 No Hit GGGAGTATTATGAAAGGAGATTAATCATAGATTATCAAAAACCCTAGAAA 6 0.15 No Hit GGATTTGAAGAAAAAAAAGACTTCGATTCATTTTCTATTTATTTCGTTAG 6 0.15 No Hit GGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGTGT 6 0.15 No Hit GGAATCGATCGTGATTTAGAACCTGTTCTTTACATGAACCCTCTTAACTA 6 0.15 No Hit GGACTGGGTATCCATGCCAGGTGTTATACCAGTAGCTTCAGGTGGTATTC 6 0.15 No Hit GGGCTTGGCTACAATTCCGAATACGCTATTACATGTTTGCGCTAGTTTTT 6 0.15 No Hit GGGGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTG 6 0.15 No Hit GAGGTAAAGGAGCACAACAAGGTAATTTGCCCGTCCCAGAAGGTTGCACT 5 0.125 No Hit GGCGCGATCTTGCTCGTGAAGGTAATGAAATTATCCGAGCAGCTTGCAAA 5 0.125 No Hit GGTAATGAAATTATCCGAGCAGCTTGCAAATGGAGTCCTGAACTAGCCGC 5 0.125 No Hit GGGCTGATATCATCAACCGTGCTAATCTTGGTATGGAAGTAATGCACGAA 5 0.125 No Hit GGGCCGGGTATCTCTCTGCATGTACGTATGCCTGTAATGTACGTAGCTAC 5 0.125 No Hit GTACAAGCTCGTAACGAAGGGCGCGATCTTGCTCGTGAAGGTAATGAAAT 5 0.125 No Hit GGGTGTCAATATATGATGATGTGTTGTTATAATGTACGCGCCTGCAAACT 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.0125 0.0 0.0 0.0 66-67 0.025 0.025 0.0 0.0 0.0 68-69 0.025 0.025 0.0 0.0 0.0 70-71 0.025 0.025 0.0 0.0 0.0 72-73 0.025 0.025 0.0 0.0 0.0 74-75 0.025 0.025 0.0 0.0 0.0 76-77 0.025 0.025 0.0 0.0 0.0 78-79 0.025 0.025 0.0 0.0 0.0 80-81 0.025 0.025 0.0 0.0 0.0 >>END_MODULE >>Kmer Content pass >>END_MODULE Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105830 READS because READLEN < 1 Read 105830 spots for ERR6133426.sra Written 105830 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra Rejected 105816 READS because READLEN < 1 Read 105816 spots for ERR6133426.sra Written 105816 spots for ERR6133426.sra SRR ids: ['ERR6133426.sra'] extra args: ['--split-files', '--defline-qual', '+'] tempdir: /tmp/pfd_815dzxw2 ERR6133426.sra spots: 2116334 blocks: [[1, 105816], [105817, 211632], [211633, 317448], [317449, 423264], [423265, 529080], [529081, 634896], [634897, 740712], [740713, 846528], [846529, 952344], [952345, 1058160], [1058161, 1163976], [1163977, 1269792], [1269793, 1375608], [1375609, 1481424], [1481425, 1587240], [1587241, 1693056], [1693057, 1798872], [1798873, 1904688], [1904689, 2010504], [2010505, 2116334]] ERR6133426 file size 468051 ERR6133426 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 ERR6133426 ERR6133426_1.fastq Input file: ERR6133426_1.fastq trimmed: ERR6133426-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 05:02:03 2024 >> started Sat Dec 7 05:02:04 2024 >> done (1.475s) 2116334 reads processed; of these: 1278 ( 0.06%) short reads filtered out after trimming by size control 33 ( 0.00%) empty reads filtered out after trimming by size control 2115023 (99.94%) reads available; of these: 27892 ( 1.32%) trimmed reads available after processing 2087131 (98.68%) untrimmed reads available after processing Length distribution of reads after trimming: length count percentage 18 23 0.00% 19 17 0.00% 20 13 0.00% 21 8 0.00% 22 4 0.00% 23 11 0.00% 24 6 0.00% 25 7 0.00% 26 7 0.00% 27 7 0.00% 28 6 0.00% 29 57 0.00% 30 8 0.00% 31 4 0.00% 32 10 0.00% 33 3 0.00% 34 2 0.00% 35 3 0.00% 36 6 0.00% 37 4 0.00% 38 9 0.00% 39 14 0.00% 40 10 0.00% 41 10 0.00% 42 13 0.00% 43 13 0.00% 44 7 0.00% 45 11 0.00% 46 13 0.00% 47 14 0.00% 48 19 0.00% 49 20 0.00% 50 16 0.00% 51 16 0.00% 52 20 0.00% 53 20 0.00% 54 23 0.00% 55 14 0.00% 56 18 0.00% 57 18 0.00% 58 9 0.00% 59 15 0.00% 60 11 0.00% 61 4 0.00% 62 0 0.00% 63 2 0.00% 64 0 0.00% 65 1 0.00% 66 4 0.00% 67 5 0.00% 68 13 0.00% 69 25 0.00% 70 2204 0.10% 71 1904 0.09% 72 1899 0.09% 73 1779 0.08% 74 1988 0.09% 75 2006 0.09% 76 1776 0.08% 77 1848 0.09% 78 2121 0.10% 79 2484 0.12% 80 2118 0.10% 81 2248 0.11% 82 2492 0.12% 83 2625 0.12% 84 2138 0.10% 85 51 0.00% 86 131 0.01% 87 183 0.01% 88 337 0.02% 89 698 0.03% 90 1625 0.08% 91 4739 0.22% 92 19061 0.90% 93 2055975 97.21% 2115023 reads passed initial QC criterion=sequence-density sequence-density=0.31 sequence-density-rank=1 fanout-score=2.85 fanout-score-rank=25 prefix-density=0.38 prefix-fanout=2.4 sequence=TGTACATTTGAACCCTGACTACACATATACACACATATACATGTAATATTATACAATCTGTCGAGTATGTGTTGGTTCAT criterion=fanout-score sequence-density=0.02 sequence-density-rank=27 fanout-score=49.86 fanout-score-rank=1 prefix-density=0.24 prefix-fanout=4.3 sequence=GCTCAAGGACGTCGGGTACCCGGTGAAGCTGGGGAGCGGCGTGGCAGCTGCCGCCGCGTACCTTTCCAATGCCACGCCCCTCATCCCGTCCAGGATCTGATCTGATCCAGCAAATTCGACGAAATTCATCGGGATTACGACGAGAAAGAAGAAGAA Started job on | Dec 07 05:02:21 Started mapping on | Dec 07 05:02:21 Finished on | Dec 07 05:02:28 Mapping speed, Million of reads per hour | 1087.73 Number of input reads | 2115023 Average input read length | 92 UNIQUE READS: Uniquely mapped reads number | 1537624 Uniquely mapped reads % | 72.70% Average mapped length | 92.44 Number of splices: Total | 95004 Number of splices: Annotated (sjdb) | 81186 Number of splices: GT/AG | 92053 Number of splices: GC/AG | 2161 Number of splices: AT/AC | 32 Number of splices: Non-canonical | 758 Mismatch rate per base, % | 0.37% Deletion rate per base | 0.03% Deletion average length | 1.87 Insertion rate per base | 0.02% Insertion average length | 2.09 MULTI-MAPPING READS: Number of reads mapped to multiple loci | 538814 % of reads mapped to multiple loci | 25.48% Number of reads mapped to too many loci | 5913 % of reads mapped to too many loci | 0.28% UNMAPPED READS: % of reads unmapped: too many mismatches | 0.00% % of reads unmapped: too short | 1.53% % of reads unmapped: other | 0.02% CHIMERIC READS: Number of chimeric reads | 0 % of chimeric reads | 0.00% N_unmapped 38585 38585 38585 N_multimapping 538814 538814 538814 N_noFeature 78911 92826 1469451 N_ambiguous 59969 5767 169 UnstrandedReadsAssigned:1398744 PositiveStrandReadsAssigned:1439031 NegativeStrandReadsAssigned:68004 Dataset is classified positive stranded MeadianReadLen=93 20thPercentileLength=93 echo kmer=89 ERR6133426 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: ERR6133426-trimmed.fastq [quant] finding pseudoalignments for the reads ... done [quant] processed 2,115,023 reads, 1,858,862 reads pseudoaligned [ em] quantifying the abundances ... done [ em] the Expectation-Maximization algorithm ran for 928 rounds 52973 ERR6133426.ke.tsv 35125 ERR6133426.se.tsv 88098 total ==> ERR6133426.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 39 20.3179 PNS24243 293 194 0 0 KQK14069 1603 1504 10 4.75248 KQK14071 474 375 0 0 ==> ERR6133426.se.tsv <== BRADI_1g14170v3 9 BRADI_1g53295v3 10 BRADI_1g59795v3 12 BRADI_1g07683v3 0 BRADI_1g00485v3 0 BRADI_1g20270v3 26 BRADI_1g74790v3 24 BRADI_1g09890v3 0 BRADI_1g77505v3 47 BRADI_1g48960v3 0 ERR6133426 completed mapping pipeline successfully