Starting /dee2/code/volunteer_pipeline.sh ERR6133422 current disk space = 1546945880064 free memory = 1421009236 ERR6133422 SRAfilesize 8652a4b8a4c02f2d7487b321cdf3c2c3 ERR6133422.sra ERR6133422.sra file validated ERR6133422 is single end ERR6133422 is conventional basespace ERR6133422 read1 length is 70-93 nt ##FastQC 0.11.5 >>Basic Statistics pass #Measure Value Filename ERR6133422_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 35.23275 37.0 33.0 37.0 33.0 37.0 2 36.45375 37.0 37.0 37.0 37.0 37.0 3 35.69 37.0 37.0 37.0 33.0 37.0 4 35.506 37.0 37.0 37.0 33.0 37.0 5 35.32275 37.0 37.0 37.0 33.0 37.0 6 35.60275 37.0 37.0 37.0 33.0 37.0 7 37.1595 37.0 37.0 40.0 33.0 40.0 8 37.1985 37.0 37.0 40.0 33.0 40.0 9 37.23975 37.0 37.0 40.0 33.0 40.0 10-11 37.236125 37.0 37.0 40.0 33.0 40.0 12-13 37.096125 37.0 37.0 40.0 33.0 40.0 14-15 37.05175 37.0 37.0 40.0 33.0 40.0 16-17 36.9735 37.0 37.0 40.0 33.0 40.0 18-19 37.142250000000004 37.0 37.0 40.0 33.0 40.0 20-21 37.33425 37.0 37.0 40.0 33.0 40.0 22-23 37.0415 37.0 37.0 40.0 33.0 40.0 24-25 37.088499999999996 37.0 37.0 40.0 33.0 40.0 26-27 36.926375 37.0 37.0 40.0 33.0 40.0 28-29 37.099374999999995 37.0 37.0 40.0 33.0 40.0 30-31 37.083875000000006 37.0 37.0 40.0 33.0 40.0 32-33 36.957125000000005 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.839625 37.0 37.0 40.0 33.0 40.0 38-39 36.738749999999996 37.0 37.0 40.0 33.0 40.0 40-41 36.565 37.0 37.0 40.0 33.0 40.0 42-43 36.48075 37.0 37.0 40.0 33.0 40.0 44-45 36.184250000000006 37.0 37.0 38.5 33.0 40.0 46-47 36.166 37.0 37.0 37.0 33.0 40.0 48-49 35.926625 37.0 35.0 37.0 33.0 40.0 50-51 35.732375000000005 37.0 33.0 37.0 33.0 40.0 52-53 35.649375 37.0 33.0 37.0 33.0 40.0 54-55 35.473625 37.0 33.0 37.0 33.0 37.0 56-57 35.40625 37.0 33.0 37.0 33.0 37.0 58-59 34.953500000000005 37.0 33.0 37.0 33.0 37.0 60-61 34.747375000000005 37.0 33.0 37.0 33.0 37.0 62-63 34.888999999999996 37.0 33.0 37.0 33.0 37.0 64-65 34.361625000000004 37.0 33.0 37.0 30.0 37.0 66-67 34.358999999999995 37.0 33.0 37.0 33.0 37.0 68-69 32.963 35.0 33.0 35.0 30.0 37.0 70-71 33.25995781210908 33.0 33.0 37.0 27.0 37.0 72-73 34.08786893281664 37.0 33.0 37.0 27.0 37.0 74-75 34.33541686419519 37.0 33.0 37.0 33.0 37.0 76-77 34.24143836040791 37.0 33.0 37.0 30.0 37.0 78-79 34.25957309173459 37.0 33.0 37.0 33.0 37.0 80-81 34.295262990841884 37.0 33.0 37.0 33.0 37.0 82-83 34.03632002995453 37.0 33.0 37.0 30.0 37.0 84-85 34.0313806125538 37.0 33.0 37.0 30.0 37.0 86-87 33.83447749809306 37.0 33.0 37.0 27.0 37.0 88-89 33.876048817696415 37.0 33.0 37.0 27.0 37.0 90-91 33.65280956013221 35.0 33.0 37.0 27.0 37.0 92-93 33.27739638952454 33.0 33.0 37.0 27.0 37.0 >>END_MODULE >>Per sequence quality scores pass #Quality Count 20 7.0 21 10.0 22 4.0 23 11.0 24 20.0 25 33.0 26 20.0 27 27.0 28 31.0 29 69.0 30 79.0 31 114.0 32 134.0 33 158.0 34 271.0 35 552.0 36 1024.0 37 1059.0 38 374.0 39 3.0 >>END_MODULE >>Per base sequence content fail #Base G A T C 1 90.8 2.0 1.875 5.325 2 76.47058823529412 13.967459324155193 5.657071339173967 3.9048811013767213 3 41.175 37.574999999999996 11.15 10.100000000000001 4 37.275000000000006 28.7 15.2 18.825 5 25.575 32.925 22.7 18.8 6 20.0 39.550000000000004 24.65 15.8 7 37.275000000000006 28.725 18.95 15.049999999999999 8 30.925000000000004 29.849999999999998 22.7 16.525000000000002 9 27.325 26.825 26.775 19.075 10-11 25.937500000000004 28.3875 25.974999999999998 19.7 12-13 29.4375 25.8625 26.0125 18.6875 14-15 24.025 28.849999999999998 28.375 18.75 16-17 24.7875 31.775 25.0 18.4375 18-19 24.712500000000002 26.174999999999997 25.587500000000002 23.525 20-21 25.275 28.225 25.6 20.9 22-23 28.425 23.925 25.275 22.375 24-25 27.1125 23.3875 28.1625 21.337500000000002 26-27 25.087500000000002 26.950000000000003 29.3875 18.575 28-29 25.724999999999998 27.3375 26.3625 20.575 30-31 28.9875 26.187500000000004 23.275000000000002 21.55 32-33 23.925 27.925 27.175 20.974999999999998 34-35 26.8125 24.2625 27.375 21.55 36-37 25.112499999999997 22.900000000000002 28.8875 23.1 38-39 26.8125 24.0375 30.975 18.175 40-41 27.0875 26.187500000000004 25.2625 21.462500000000002 42-43 25.5375 28.5625 24.875 21.025 44-45 24.587500000000002 25.362499999999997 28.287499999999998 21.762500000000003 46-47 24.5625 23.575 26.8625 25.0 48-49 24.85 25.374999999999996 30.112499999999997 19.662499999999998 50-51 25.224999999999998 26.337500000000002 29.062500000000004 19.375 52-53 25.3 27.237499999999997 25.95 21.512500000000003 54-55 24.5 28.9125 26.8625 19.725 56-57 25.95 25.0125 28.012500000000003 21.025 58-59 23.3375 23.549999999999997 29.525000000000002 23.5875 60-61 24.825 24.25 28.95 21.975 62-63 22.8625 27.9375 30.825000000000003 18.375 64-65 25.324999999999996 26.0 28.000000000000004 20.674999999999997 66-67 26.2125 26.8375 26.9625 19.9875 68-69 24.1375 26.337500000000002 28.1625 21.3625 70-71 26.985119419782418 23.308740777791673 26.947605352007002 22.758534450418907 72-73 27.414095811387007 23.99046902432907 26.912465512916977 21.682969651366943 74-75 23.670647391577624 27.052168447517282 28.77435575109994 20.502828409805154 76-77 22.38467355684396 26.354928157297707 28.91353667759012 22.346861608268213 78-79 23.406940063091483 23.735015772870664 30.41009463722397 22.44794952681388 80-81 24.68682778691636 26.93913703656839 29.09021890421359 19.283816272301657 82-83 24.501966755487885 24.26088059890877 28.60043141733283 22.636721228270524 84-85 23.91912512716175 23.52492370295015 30.328077314343844 22.227873855544253 86-87 22.19679633867277 26.264937706585307 32.3162979913552 19.221967963386728 88-89 22.489193999491484 29.010933129926265 29.60844139333842 18.891431477243835 90-91 27.45995423340961 25.108060005085175 28.324434274091026 19.107551487414188 92-93 22.934146961606917 28.540554284261376 29.60844139333842 18.916857360793287 >>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.5 17 5.5 18 6.5 19 1.5 20 0.5 21 2.5 22 2.5 23 2.0 24 3.0 25 2.5 26 5.0 27 9.0 28 14.5 29 20.5 30 22.0 31 23.5 32 30.5 33 45.0 34 64.0 35 75.5 36 81.0 37 103.5 38 130.0 39 143.5 40 162.5 41 186.5 42 194.5 43 191.5 44 183.5 45 176.5 46 220.5 47 229.5 48 187.5 49 179.0 50 178.5 51 173.0 52 168.0 53 196.5 54 178.0 55 103.0 56 73.0 57 72.0 58 62.0 59 47.0 60 42.0 61 40.5 62 39.5 63 44.5 64 44.5 65 36.0 66 31.0 67 21.5 68 13.5 69 10.5 70 4.5 71 2.0 72 2.5 73 1.5 74 0.0 75 1.0 76 1.0 77 0.5 78 0.5 79 0.0 80 0.0 81 0.0 82 0.5 83 0.5 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.125 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.05083884087442806 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 6.0 72 8.0 73 3.0 74 5.0 75 6.0 76 4.0 77 2.0 78 1.0 79 8.0 80 5.0 81 6.0 82 5.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 3933.0 >>END_MODULE >>Sequence Duplication Levels pass #Total Deduplicated Percentage 79.05 #Duplication Level Percentage of deduplicated Percentage of total 1 92.40986717267552 73.05 2 4.427577482605946 7.000000000000001 3 1.1385199240986716 2.7 4 0.6008855154965212 1.9 5 0.31625553447185323 1.25 6 0.2213788741302973 1.05 7 0.06325110689437065 0.35000000000000003 8 0.15812776723592661 1.0 9 0.09487666034155598 0.675 >10 0.5376344086021506 9.0 >50 0.031625553447185324 2.025 >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 81 2.025 No Hit GGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGTG 50 1.25 No Hit GGGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTGG 50 1.25 No Hit GGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTT 35 0.8750000000000001 No Hit GGAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAA 32 0.8 No Hit GGGAGAGCAATACAAGCGTTGTGCTGCTAGGCGAAGCGGTTGAGTGCCGC 22 0.5499999999999999 No Hit GGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGTGTACAAGCTCGTAA 21 0.525 No Hit GGGGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTG 21 0.525 No Hit GGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTGGG 18 0.44999999999999996 No Hit GGTCGCGTTATTAATACTTGGGCTGATATCATCAACCGTGCTAATCTTGG 17 0.42500000000000004 No Hit GGGAGTATTATGAAAGGAGATTAATCATAGATTATCAAAAACCCTAGAAA 14 0.35000000000000003 No Hit GGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAATCG 13 0.325 No Hit GGCGTTCAACCTAAATGGATTCAATTTCAACCAATCTGTAGTTGATAGTC 13 0.325 No Hit GGGCGCGATCTTGCTCGTGAAGGTAATGAAATTATCCGAGCAGCTTGCAA 13 0.325 No Hit GGAGTCCCATATATATATGTATAAGATGCCAGCCCGGTTTCGTCAACCAT 11 0.27499999999999997 No Hit GGTGGAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGC 10 0.25 No Hit GGAGTATCCTTGATATATAAATATATAGATAAATAGATTTAAATGGAAAA 10 0.25 No Hit GGGGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTG 10 0.25 No Hit GGAAGAGCCCGAGCTGCTGCAGCAGGTTTTGAAAAGGGAATCGATCGTGA 9 0.22499999999999998 No Hit GGGCAAGGAGAAGTACAAGTGCGGATCCAACGTCTTCTGGAAATGGTGAA 9 0.22499999999999998 No Hit GGGAATCGATCGTGATTTAGAACCTGTTCTTTACATGAACCCTCTTAACT 9 0.22499999999999998 No Hit GGGATGGAGCGACAGAAGTATGGAAATAGGATAAGGTAGCGGCGAGACGA 8 0.2 No Hit GGATTCAATTTCAACCAATCTGTAGTTGATAGTCAAGGTCGCGTTATTAA 8 0.2 No Hit GGTAATGAAATTATCCGAGCAGCTTGCAAATGGAGTCCTGAACTAGCCGC 8 0.2 No Hit GGGCTTGGCTACAATTCCGAATACGCTATTACATGTTTGCGCTAGTTTTT 8 0.2 No Hit GGGGAACGCGAAATCACTTTAGGTTTTGTTGATTTATTGCGCGACGATTT 8 0.2 No Hit GGCGCGATCTTGCTCGTGAAGGTAATGAAATTATCCGAGCAGCTTGCAAA 7 0.17500000000000002 No Hit GGGCCGGGTATCTCTCTGCATGTACGTATGCCTGTAATGTACGTAGCTAC 7 0.17500000000000002 No Hit GGGGATGGAGCGACAGAAGTATGGAAATAGGATAAGGTAGCGGCGAGACG 6 0.15 No Hit GAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAAT 6 0.15 No Hit TACCTAGGCACCCAGAGACGAGGAAGGGCGTAGCAAGCGACGAAATGCTT 6 0.15 No Hit GAAGGTAATGAAATTATCCGAGCAGCTTGCAAATGGAGTCCTGAACTAGC 6 0.15 No Hit GGGCTCGAGGAGCATATGTACATTTGAACCCTGACTACACATATACACAC 6 0.15 No Hit GGGAACGCGAAATCACTTTAGGTTTTGTTGATTTATTGCGCGACGATTTT 6 0.15 No Hit GGTTTTGAAAAGGGAATCGATCGTGATTTAGAACCTGTTCTTTACATGAA 6 0.15 No Hit GGGCTGCGAGGAATCCGGCAAGGCATAAACAACCAAGGACAGCTCCGTAT 5 0.125 No Hit GGCAAGGAGAAGTACAAGTGCGGATCCAACGTCTTCTGGAAATGGTGAAG 5 0.125 No Hit GGAATAAGAATAAATCGCAACTCCTTTCCACTACACATAAAAATTGATTT 5 0.125 No Hit GGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGTGT 5 0.125 No Hit GGAAACAACGAGGTCATCATCGACATGATATATTGCTGCTATTTTCCACC 5 0.125 No Hit GGGCTGATATCATCAACCGTGCTAATCTTGGTATGGAAGTAATGCACGAA 5 0.125 No Hit GTGGAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCT 5 0.125 No Hit GTACAAGCTCGTAACGAAGGGCGCGATCTTGCTCGTGAAGGTAATGAAAT 5 0.125 No Hit GGGCCATTTGTGGCATGCAGGAAGAGCCCGAGCTGCTGCAGCAGGTTTTG 5 0.125 No Hit GGAACGCGAAATCACTTTAGGTTTTGTTGATTTATTGCGCGACGATTTTA 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.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 GCGATCT 25 0.0023441382 34.8 72-73 GCGCGAT 25 0.0023441382 34.8 70-71 GTGAAGG 25 0.0023441382 34.8 84-85 GGGCGCG 25 0.0023441382 34.8 68-69 TTGCTCG 25 0.0023441382 34.8 78-79 TCGTGAA 30 0.0057259556 29.0 82-83 ACCTGGT 30 0.0057259556 29.0 10-11 CTGGTGC 30 0.0057259556 29.0 12-13 GCTCGTG 30 0.0057259556 29.0 80-81 GATCTTG 30 0.0057259556 29.0 74-75 TCTTGCT 30 0.0057259556 29.0 76-77 GAAGGTA 30 0.0057259556 29.0 86-87 >>END_MODULE Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139240 READS because READLEN < 1 Read 139240 spots for ERR6133422.sra Written 139240 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra Rejected 139229 READS because READLEN < 1 Read 139229 spots for ERR6133422.sra Written 139229 spots for ERR6133422.sra SRR ids: ['ERR6133422.sra'] extra args: ['--split-files', '--defline-qual', '+'] tempdir: /tmp/pfd_6x0l25y3 ERR6133422.sra spots: 2784591 blocks: [[1, 139229], [139230, 278458], [278459, 417687], [417688, 556916], [556917, 696145], [696146, 835374], [835375, 974603], [974604, 1113832], [1113833, 1253061], [1253062, 1392290], [1392291, 1531519], [1531520, 1670748], [1670749, 1809977], [1809978, 1949206], [1949207, 2088435], [2088436, 2227664], [2227665, 2366893], [2366894, 2506122], [2506123, 2645351], [2645352, 2784591]] ERR6133422 file size 616121 ERR6133422 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 ERR6133422 ERR6133422_1.fastq Input file: ERR6133422_1.fastq trimmed: ERR6133422-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 04:48:53 2024 >> started Sat Dec 7 04:48:55 2024 >> done (1.913s) 2784591 reads processed; of these: 1449 ( 0.05%) short reads filtered out after trimming by size control 15 ( 0.00%) empty reads filtered out after trimming by size control 2783127 (99.95%) reads available; of these: 37510 ( 1.35%) trimmed reads available after processing 2745617 (98.65%) untrimmed reads available after processing Length distribution of reads after trimming: length count percentage 18 35 0.00% 19 36 0.00% 20 17 0.00% 21 18 0.00% 22 14 0.00% 23 12 0.00% 24 13 0.00% 25 11 0.00% 26 9 0.00% 27 9 0.00% 28 12 0.00% 29 46 0.00% 30 7 0.00% 31 5 0.00% 32 11 0.00% 33 8 0.00% 34 2 0.00% 35 9 0.00% 36 13 0.00% 37 10 0.00% 38 11 0.00% 39 23 0.00% 40 39 0.00% 41 18 0.00% 42 14 0.00% 43 21 0.00% 44 18 0.00% 45 20 0.00% 46 24 0.00% 47 9 0.00% 48 22 0.00% 49 17 0.00% 50 25 0.00% 51 35 0.00% 52 23 0.00% 53 24 0.00% 54 19 0.00% 55 22 0.00% 56 23 0.00% 57 8 0.00% 58 14 0.00% 59 15 0.00% 60 27 0.00% 61 11 0.00% 62 1 0.00% 63 0 0.00% 64 3 0.00% 65 0 0.00% 66 6 0.00% 67 4 0.00% 68 15 0.00% 69 38 0.00% 70 4080 0.15% 71 3243 0.12% 72 3421 0.12% 73 3085 0.11% 74 3547 0.13% 75 3521 0.13% 76 3120 0.11% 77 3151 0.11% 78 3611 0.13% 79 4246 0.15% 80 3746 0.13% 81 3601 0.13% 82 4266 0.15% 83 4666 0.17% 84 3470 0.12% 85 71 0.00% 86 140 0.01% 87 225 0.01% 88 464 0.02% 89 876 0.03% 90 2129 0.08% 91 6197 0.22% 92 25844 0.93% 93 2691561 96.71% 2783127 reads passed initial QC criterion=sequence-density sequence-density=0.22 sequence-density-rank=1 fanout-score=7.09 fanout-score-rank=16 prefix-density=0.55 prefix-fanout=2.8 sequence=GCCGCCGCCGCCAAGGAAGGCATGTTCGTCAAGAACTACAGCTACTGATCCTAATCGCATCAAGCTTCAACGCCTGTGAGTGAAAACCAGTGATGAGAGTGCTGCTGCTAGCTAGCGCCGGCATTGATGAGCTTGAGAGGGCACTGTAGCCAGTGTGTCAGTCGTTGTTAAATTACAGGTTGAGATCATCAGCGTACTCCGATGGGAGGTGGACATCAGAAAGTATACTGTGTTTTACCACCCTAATTAAGTAAACAACTTTTGGAATGGTGATCATCTTACTGTTTTAAATAAATCTGTTTC criterion=fanout-score sequence-density=0.02 sequence-density-rank=21 fanout-score=345.01 fanout-score-rank=1 prefix-density=0.71 prefix-fanout=11.1 sequence=AAGAAGAAGGTGCAGACCGAGTGCGCCTCCATGCCTTTCGATGACCAATGCGCCGTCTTGGAAAAGGAGGCCGTGAACGTGTCCCTCGAGAACCTCAAGACCTACCCGTTCGTCAAGGAAGGCGTCGCCAACGGAACCCTCAAGCTCGTGGGCGGCCACTACGACTTCGTCTCCGGCAAGTTCGACACATGGGAGCTCTAAGTCCTCTCATCCGGTTAACTCCTATACATACAACGTATACTTATACATACAGATATGGAGATGACCCTACAGATCGATCCATTGATGTGGATGCGATGCCATGGAGTATATGTACTCGCTATTTTCCAGTACTGCATGCCGGATGGCTCGATGTGAATTTGTAATAAGCAATAGAAGTTTCTACCATTTTCTGACGCGGGGTTGTACTTGTGATGCGTAATTTGGTCGTCTTGTGACCAAAAGACATCAACTATATAATTATAATACAATTTTCATCAAGACTCTGTATT Started job on | Dec 07 04:49:20 Started mapping on | Dec 07 04:49:21 Finished on | Dec 07 04:49:25 Mapping speed, Million of reads per hour | 2504.81 Number of input reads | 2783127 Average input read length | 92 UNIQUE READS: Uniquely mapped reads number | 2013711 Uniquely mapped reads % | 72.35% Average mapped length | 92.31 Number of splices: Total | 138892 Number of splices: Annotated (sjdb) | 117383 Number of splices: GT/AG | 133756 Number of splices: GC/AG | 3395 Number of splices: AT/AC | 38 Number of splices: Non-canonical | 1703 Mismatch rate per base, % | 0.47% Deletion rate per base | 0.05% Deletion average length | 1.81 Insertion rate per base | 0.02% Insertion average length | 1.87 MULTI-MAPPING READS: Number of reads mapped to multiple loci | 711544 % of reads mapped to multiple loci | 25.57% Number of reads mapped to too many loci | 6912 % of reads mapped to too many loci | 0.25% UNMAPPED READS: % of reads unmapped: too many mismatches | 0.00% % of reads unmapped: too short | 1.81% % of reads unmapped: other | 0.02% CHIMERIC READS: Number of chimeric reads | 0 % of chimeric reads | 0.00% N_unmapped 57872 57872 57872 N_multimapping 711544 711544 711544 N_noFeature 104899 122274 1920402 N_ambiguous 83840 8045 203 UnstrandedReadsAssigned:1824972 PositiveStrandReadsAssigned:1883392 NegativeStrandReadsAssigned:93106 Dataset is classified positive stranded MeadianReadLen=93 20thPercentileLength=93 echo kmer=89 ERR6133422 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: ERR6133422-trimmed.fastq [quant] finding pseudoalignments for the reads ... done [quant] processed 2,783,127 reads, 2,446,765 reads pseudoaligned [ em] quantifying the abundances ... done [ em] the Expectation-Maximization algorithm ran for 974 rounds 52973 ERR6133422.ke.tsv 35125 ERR6133422.se.tsv 88098 total ==> ERR6133422.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 69 27.2506 PNS24243 293 194 0 0 KQK14069 1603 1504 28 10.0877 KQK14071 474 375 0 0 ==> ERR6133422.se.tsv <== BRADI_1g14170v3 28 BRADI_1g53295v3 62 BRADI_1g59795v3 10 BRADI_1g07683v3 0 BRADI_1g00485v3 0 BRADI_1g20270v3 15 BRADI_1g74790v3 22 BRADI_1g09890v3 0 BRADI_1g77505v3 56 BRADI_1g48960v3 0 ERR6133422 completed mapping pipeline successfully