Starting /dee2/code/volunteer_pipeline.sh ERR6133570 current disk space = 1544439767040 free memory = 1599413520 ERR6133570 SRAfilesize 10fd8ce4d3287edbeeef849fddef7cd7 ERR6133570.sra ERR6133570.sra file validated ERR6133570 is single end ERR6133570 is conventional basespace ERR6133570 read1 length is 70-93 nt ##FastQC 0.11.5 >>Basic Statistics pass #Measure Value Filename ERR6133570_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.29875 37.0 33.0 37.0 33.0 37.0 2 36.4865 37.0 37.0 37.0 37.0 37.0 3 35.98325 37.0 37.0 37.0 33.0 37.0 4 35.46875 37.0 37.0 37.0 33.0 37.0 5 35.32975 37.0 37.0 37.0 33.0 37.0 6 35.7485 37.0 37.0 37.0 33.0 37.0 7 37.5765 40.0 37.0 40.0 33.0 40.0 8 37.61075 40.0 37.0 40.0 33.0 40.0 9 37.6725 40.0 37.0 40.0 33.0 40.0 10-11 37.512375 40.0 37.0 40.0 33.0 40.0 12-13 37.455875 37.0 37.0 40.0 33.0 40.0 14-15 37.399625 37.0 37.0 40.0 33.0 40.0 16-17 37.248999999999995 37.0 37.0 40.0 33.0 40.0 18-19 37.287625 37.0 37.0 40.0 33.0 40.0 20-21 37.11375 37.0 37.0 40.0 33.0 40.0 22-23 37.119749999999996 37.0 37.0 40.0 33.0 40.0 24-25 37.266875 37.0 37.0 40.0 33.0 40.0 26-27 37.131 37.0 37.0 40.0 33.0 40.0 28-29 37.185249999999996 37.0 37.0 40.0 33.0 40.0 30-31 37.020624999999995 37.0 37.0 40.0 33.0 40.0 32-33 36.779125 37.0 37.0 40.0 33.0 40.0 34-35 36.748375 37.0 37.0 40.0 33.0 40.0 36-37 36.64975 37.0 37.0 40.0 33.0 40.0 38-39 36.332750000000004 37.0 37.0 40.0 33.0 40.0 40-41 36.135 37.0 37.0 40.0 33.0 40.0 42-43 36.0715 37.0 37.0 40.0 33.0 40.0 44-45 35.757000000000005 37.0 33.0 40.0 33.0 40.0 46-47 35.5055 37.0 33.0 37.0 33.0 40.0 48-49 35.461124999999996 37.0 33.0 37.0 33.0 40.0 50-51 35.25725 37.0 33.0 37.0 33.0 40.0 52-53 34.993875 37.0 33.0 37.0 30.0 40.0 54-55 34.858000000000004 37.0 33.0 37.0 27.0 40.0 56-57 34.62375 37.0 33.0 37.0 27.0 37.0 58-59 33.364875 35.0 33.0 37.0 27.0 37.0 60-61 34.06325 37.0 33.0 37.0 27.0 37.0 62-63 34.149249999999995 37.0 33.0 37.0 27.0 37.0 64-65 34.082750000000004 37.0 33.0 37.0 27.0 37.0 66-67 34.1775 37.0 33.0 37.0 27.0 37.0 68-69 33.5125 35.0 33.0 37.0 27.0 37.0 70-71 33.60411479719579 35.0 33.0 37.0 27.0 37.0 72-73 33.98901148793836 37.0 33.0 37.0 27.0 37.0 74-75 33.786870942381825 37.0 33.0 37.0 27.0 37.0 76-77 33.852130426137705 37.0 33.0 37.0 27.0 37.0 78-79 33.84619554179493 37.0 33.0 37.0 27.0 37.0 80-81 33.8764933435833 37.0 33.0 37.0 27.0 37.0 82-83 33.70062833345031 37.0 33.0 37.0 27.0 37.0 84-85 33.504258245225714 37.0 33.0 37.0 27.0 37.0 86-87 33.36813611755607 37.0 33.0 37.0 27.0 37.0 88-89 33.522170662541896 37.0 33.0 37.0 27.0 37.0 90-91 33.35202371745295 35.0 33.0 37.0 27.0 37.0 92-93 33.153905645784995 37.0 33.0 37.0 27.0 37.0 >>END_MODULE >>Per sequence quality scores pass #Quality Count 20 18.0 21 15.0 22 19.0 23 25.0 24 20.0 25 24.0 26 31.0 27 44.0 28 63.0 29 67.0 30 89.0 31 131.0 32 141.0 33 205.0 34 269.0 35 420.0 36 804.0 37 991.0 38 607.0 39 17.0 >>END_MODULE >>Per base sequence content fail #Base G A T C 1 87.0 2.475 4.125 6.4 2 68.45 18.175 8.5 4.875 3 33.175 38.15 16.625 12.049999999999999 4 35.15 25.45 18.475 20.925 5 23.275000000000002 32.725 25.75 18.25 6 18.925 39.375 26.174999999999997 15.525 7 35.775 29.25 19.625 15.35 8 31.2 29.125 23.275000000000002 16.400000000000002 9 25.374999999999996 26.900000000000002 29.075 18.65 10-11 24.887500000000003 27.200000000000003 28.487499999999997 19.425 12-13 26.887499999999996 25.025 27.800000000000004 20.2875 14-15 23.400000000000002 29.2375 28.7 18.6625 16-17 23.5125 31.4 25.2125 19.875 18-19 23.3625 26.900000000000002 27.537499999999998 22.2 20-21 25.5625 25.7375 27.187499999999996 21.512500000000003 22-23 27.212500000000002 22.4875 26.887499999999996 23.4125 24-25 26.4625 24.474999999999998 28.012500000000003 21.05 26-27 25.4375 25.8125 29.9625 18.787499999999998 28-29 26.237500000000004 26.174999999999997 27.3125 20.275000000000002 30-31 26.5 26.275 26.8125 20.4125 32-33 24.1375 26.75 28.9875 20.125 34-35 25.2 25.825 27.1375 21.837500000000002 36-37 24.637500000000003 24.6875 27.8125 22.8625 38-39 26.275 25.6 30.1875 17.9375 40-41 27.6 25.575 25.825 21.0 42-43 25.374999999999996 27.800000000000004 26.8125 20.0125 44-45 24.5625 26.400000000000002 27.650000000000002 21.3875 46-47 24.825 23.45 28.4125 23.3125 48-49 25.1875 24.6125 30.075000000000003 20.125 50-51 24.9125 25.974999999999998 29.175 19.9375 52-53 25.119077463023316 27.049385810980198 27.06192028077212 20.76961644522437 54-55 25.1 28.000000000000004 27.55 19.35 56-57 25.687500000000004 26.125 28.0625 20.125 58-59 25.324999999999996 24.375 29.2 21.099999999999998 60-61 24.525 25.275 30.0875 20.1125 62-63 22.8375 27.474999999999998 30.662499999999998 19.025 64-65 24.325 27.224999999999998 28.249999999999996 20.200000000000003 66-67 24.096536201075402 27.19769913717644 29.048393147430286 19.65737151431787 68-69 23.2375 26.987499999999997 28.050000000000004 21.725 70-71 25.756817613209908 26.044533400050035 27.833375031273455 20.365273955466602 72-73 25.77164366373902 24.51693851944793 28.406524466750316 21.304893350062734 74-75 24.27294473120987 27.546267153468463 27.67216416970918 20.508623945612488 76-77 22.716923855299772 26.56210473058437 29.496584872248928 21.224386541866934 78-79 24.289700659563675 24.84779299847793 30.2130898021309 20.6494165398275 80-81 23.63659531090724 27.675840978593275 29.43425076452599 19.253312945973498 82-83 23.485721603278268 26.853630426431042 29.56844666410552 20.09220130618517 84-85 24.144144144144146 23.462033462033464 31.299871299871302 21.093951093951095 86-87 22.892498066511987 26.488785769528228 31.412735241041506 19.20598092291828 88-89 21.745295179169887 28.95076050528487 29.22144882701727 20.08249548852797 90-91 24.851765919051303 25.754060324825982 29.64681618973962 19.74735756638309 92-93 22.815158546017013 29.698375870069604 28.254704820830113 19.23176076308327 >>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 0.0 17 8.0 18 9.0 19 1.0 20 0.5 21 2.5 22 3.0 23 4.0 24 4.5 25 5.0 26 6.0 27 9.0 28 21.0 29 30.0 30 34.5 31 39.5 32 45.0 33 60.0 34 68.5 35 83.0 36 115.0 37 139.0 38 169.5 39 174.5 40 170.0 41 187.0 42 204.0 43 215.5 44 206.5 45 190.0 46 200.5 47 192.5 48 165.0 49 157.5 50 164.5 51 166.0 52 152.0 53 159.0 54 128.0 55 78.0 56 62.0 57 55.5 58 54.0 59 45.5 60 46.0 61 45.0 62 40.0 63 41.0 64 40.5 65 36.5 66 20.5 67 13.5 68 13.5 69 14.5 70 13.5 71 9.0 72 5.0 73 3.0 74 3.0 75 1.0 76 0.5 77 0.0 78 0.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.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.27499999999999997 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.0375 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.012740476493820868 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 6.0 71 7.0 72 4.0 73 8.0 74 7.0 75 11.0 76 8.0 77 3.0 78 8.0 79 9.0 80 9.0 81 10.0 82 11.0 83 8.0 84 12.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 3879.0 >>END_MODULE >>Sequence Duplication Levels pass #Total Deduplicated Percentage 81.575 #Duplication Level Percentage of deduplicated Percentage of total 1 93.65614465216059 76.4 2 3.922770456634998 6.4 3 0.8887526815813668 2.175 4 0.4290530186944529 1.4000000000000001 5 0.18387986515476554 0.75 6 0.12258657676984369 0.6 7 0.061293288384921846 0.35000000000000003 8 0.061293288384921846 0.4 9 0.0 0.0 >10 0.6435795280416794 10.2 >50 0.030646644192460923 1.325 >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 53 1.325 No Hit GGTCGCGTTATTAATACTTGGGCTGATATCATCAACCGTGCTAATCTTGG 38 0.95 No Hit TACCTAGGCACCCAGAGACGAGGAAGGGCGTAGCAAGCGACGAAATGCTT 34 0.8500000000000001 No Hit GGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGTG 32 0.8 No Hit GGATTCAATTTCAACCAATCTGTAGTTGATAGTCAAGGTCGCGTTATTAA 31 0.775 No Hit GGGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTGG 31 0.775 No Hit GGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAATCGAGTGGCTTT 24 0.6 No Hit GTAGTAGGAATCTGGTTCACTGCTTTAGGTATTAGTACTATGGCGTTCAA 24 0.6 No Hit GGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTGGG 19 0.475 No Hit GGTGCAGCAGCTAATCGAGTGGCTTTAGAAGCCTGTGTACAAGCTCGTAA 18 0.44999999999999996 No Hit GGAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGCTAA 18 0.44999999999999996 No Hit GGAGTATCCTTGATATATAAATATATAGATAAATAGATTTAAATGGAAAA 17 0.42500000000000004 No Hit GGCGTTCAACCTAAATGGATTCAATTTCAACCAATCTGTAGTTGATAGTC 16 0.4 No Hit GGGGATGATTCTGTATTACAATTTGGTGGAGGAACTTTAGGACATCCTTG 14 0.35000000000000003 No Hit GAAGGTAATGAAATTATCCGAGCAGCTTGCAAATGGAGTCCTGAACTAGC 13 0.325 No Hit GGGAGAGCAATACAAGCGTTGTGCTGCTAGGCGAAGCGGTTGAGTGCCGC 13 0.325 No Hit GGGAGTATTATGAAAGGAGATTAATCATAGATTATCAAAAACCCTAGAAA 12 0.3 No Hit GTATTTAGCCTTGCAAGGTGGTCCTTGCTGATTCACACGGGATTCCACGT 12 0.3 No Hit GGCCTGTAGTAGGAATCTGGTTCACTGCTTTAGGTATTAGTACTATGGCG 11 0.27499999999999997 No Hit GGTAATGAAATTATCCGAGCAGCTTGCAAATGGAGTCCTGAACTAGCCGC 11 0.27499999999999997 No Hit GGTGGAGGAACTTTAGGACATCCTTGGGGAAATGCACCTGGTGCAGCAGC 10 0.25 No Hit GGAGTTGAAAATAAGCATAGATCCGGAGATTCCCAAATAGGTCAACCTTT 10 0.25 No Hit GGAGTCCTGAACTAGCCGCAGCTTGTGAAGTATGGAAGGCGATCAAATTC 8 0.2 No Hit GGAAAAGAGGGGTTACTTTTTTTCATTTTTCCCTTAAAAGATAGGCTTTG 8 0.2 No Hit GGGCGCGATCTTGCTCGTGAAGGTAATGAAATTATCCGAGCAGCTTGCAA 7 0.17500000000000002 No Hit GGGGAAGTACACCAGCGACGGCGAGGCCGCCGCCGCCAAGGAAGGCATGT 7 0.17500000000000002 No Hit CAAGGTCGCGTTATTAATACTTGGGCTGATATCATCAACCGTGCTAATCT 6 0.15 No Hit GGGCTGATATCATCAACCGTGCTAATCTTGGTATGGAAGTAATGCACGAA 6 0.15 No Hit GGGAATCGATCGTGATTTAGAACCTGTTCTTTACATGAACCCTCTTAACT 6 0.15 No Hit GGTTTTGAAAAGGGAATCGATCGTGATTTAGAACCTGTTCTTTACATGAA 6 0.15 No Hit GGGAAAAGAGGGGTTACTTTTTTTCATTTTTCCCTTAAAAGATAGGCTTT 5 0.125 No Hit CGAGGAAGGGCGTAGCAAGCGACGAAATGCTTCGGGGAGTTGAAAATAAG 5 0.125 No Hit GTAGGAATCTGGTTCACTGCTTTAGGTATTAGTACTATGGCGTTCAACCT 5 0.125 No Hit CAACCTAAATGGATTCAATTTCAACCAATCTGTAGTTGATAGTCAAGGTC 5 0.125 No Hit GGAGTGGGGGTGCCATACGCAAAAAGGAAGCGACTCATAGAATGGCAGAG 5 0.125 No Hit GTACAAGCTCGTAACGAAGGGCGCGATCTTGCTCGTGAAGGTAATGAAAT 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.0125 0.0 26-27 0.0 0.0 0.0 0.025 0.0 28-29 0.0 0.0 0.0 0.025 0.0 30-31 0.0 0.0 0.0 0.025 0.0 32-33 0.0 0.0 0.0 0.025 0.0 34-35 0.0 0.0 0.0 0.025 0.0 36-37 0.0 0.0 0.0 0.025 0.0 38-39 0.0 0.0 0.0 0.025 0.0 40-41 0.0 0.0 0.0 0.025 0.0 42-43 0.0 0.0 0.0 0.025 0.0 44-45 0.0 0.0 0.0 0.025 0.0 46-47 0.0 0.0 0.0 0.025 0.0 48-49 0.0 0.0 0.0 0.025 0.0 50-51 0.0 0.0 0.0 0.025 0.0 52-53 0.0 0.0 0.0 0.025 0.0 54-55 0.0 0.0 0.0 0.025 0.0 56-57 0.0 0.0 0.0 0.025 0.0 58-59 0.0 0.0 0.0 0.025 0.0 60-61 0.0 0.0 0.0 0.025 0.0 62-63 0.0 0.0 0.0 0.025 0.0 64-65 0.0 0.0 0.0 0.025 0.0 66-67 0.0 0.0 0.0 0.025 0.0 68-69 0.0 0.0 0.0 0.025 0.0 70-71 0.0 0.0 0.0 0.025 0.0 72-73 0.025 0.0 0.0 0.025 0.0 74-75 0.025 0.0 0.0 0.025 0.0 76-77 0.025 0.0 0.0 0.025 0.0 78-79 0.025 0.0 0.0 0.025 0.0 80-81 0.025 0.0 0.0 0.025 0.0 >>END_MODULE >>Kmer Content warn #Sequence Count PValue Obs/Exp Max Max Obs/Exp Position AATCTGG 20 0.0028106743 64.77187 9 TAGGAAT 20 0.0028106743 64.77187 5 AGTAGGA 20 0.0028106743 64.77187 3 AGGAATC 20 0.0028106743 64.77187 6 TAGTAGG 20 0.0028106743 64.77187 2 GTAGGAA 20 0.0028106743 64.77187 4 GGAATCT 20 0.0028106743 64.77187 7 GTAGTAG 20 0.0028106743 64.77187 1 GAATCTG 25 0.0068019433 51.817497 8 >>END_MODULE Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra Rejected 72002 READS because READLEN < 1 Read 72002 spots for ERR6133570.sra Written 72002 spots for ERR6133570.sra Rejected 71988 READS because READLEN < 1 Read 71988 spots for ERR6133570.sra Written 71988 spots for ERR6133570.sra SRR ids: ['ERR6133570.sra'] extra args: ['--split-files', '--defline-qual', '+'] tempdir: /tmp/pfd_8z_3czlj ERR6133570.sra spots: 1439774 blocks: [[1, 71988], [71989, 143976], [143977, 215964], [215965, 287952], [287953, 359940], [359941, 431928], [431929, 503916], [503917, 575904], [575905, 647892], [647893, 719880], [719881, 791868], [791869, 863856], [863857, 935844], [935845, 1007832], [1007833, 1079820], [1079821, 1151808], [1151809, 1223796], [1223797, 1295784], [1295785, 1367772], [1367773, 1439774]] ERR6133570 file size 316802 ERR6133570 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 ERR6133570 ERR6133570_1.fastq Input file: ERR6133570_1.fastq trimmed: ERR6133570-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:52:29 2024 >> started Sat Dec 7 08:52:30 2024 >> done (0.796s) 1439774 reads processed; of these: 95 ( 0.01%) short reads filtered out after trimming by size control 28 ( 0.00%) empty reads filtered out after trimming by size control 1439651 (99.99%) reads available; of these: 31555 ( 2.19%) trimmed reads available after processing 1408096 (97.81%) untrimmed reads available after processing Length distribution of reads after trimming: length count percentage 18 11 0.00% 19 34 0.00% 20 20 0.00% 21 24 0.00% 22 22 0.00% 23 10 0.00% 24 7 0.00% 25 3 0.00% 26 1 0.00% 27 8 0.00% 28 22 0.00% 29 13 0.00% 30 9 0.00% 31 12 0.00% 32 8 0.00% 33 15 0.00% 34 14 0.00% 35 75 0.01% 36 274 0.02% 37 7 0.00% 38 16 0.00% 39 55 0.00% 40 42 0.00% 41 19 0.00% 42 5 0.00% 43 6 0.00% 44 9 0.00% 45 3 0.00% 46 3 0.00% 47 1 0.00% 48 7 0.00% 49 3 0.00% 50 9 0.00% 51 20 0.00% 52 7 0.00% 53 4 0.00% 54 2 0.00% 55 5 0.00% 56 3 0.00% 57 5 0.00% 58 7 0.00% 59 4 0.00% 60 5 0.00% 61 4 0.00% 62 2 0.00% 63 1 0.00% 64 0 0.00% 65 3 0.00% 66 3 0.00% 67 4 0.00% 68 17 0.00% 69 36 0.00% 70 3642 0.25% 71 3341 0.23% 72 3762 0.26% 73 3233 0.22% 74 3283 0.23% 75 3304 0.23% 76 3002 0.21% 77 3035 0.21% 78 3285 0.23% 79 3569 0.25% 80 3196 0.22% 81 3442 0.24% 82 3974 0.28% 83 4074 0.28% 84 3531 0.25% 85 75 0.01% 86 138 0.01% 87 203 0.01% 88 357 0.02% 89 704 0.05% 90 1435 0.10% 91 4280 0.30% 92 22762 1.58% 93 1357125 94.27% 1439651 reads passed initial QC criterion=sequence-density sequence-density=0.40 sequence-density-rank=1 fanout-score=7.04 fanout-score-rank=15 prefix-density=0.62 prefix-fanout=4.5 sequence=ATGCATGCATGC criterion=fanout-score sequence-density=0.01 sequence-density-rank=32 fanout-score=104.68 fanout-score-rank=1 prefix-density=0.18 prefix-fanout=7.4 sequence=AGGAAGAGGAGGATGCAGTCAGGGCTCACAAACAACCTGCCACTGCCGCCATTGGCCTTCTAAAACAGGAGAGGAGGGGTTAGATAGTTTCGATCTGCAAGGGGGTGGGGCAATGGGCATTCCGTTGAGTTTGTATGGTGGGTGCTGTTTTGCAGAAGCTGTATGCTTATGAGCAGCTATGGTACTTATCGAGGACAGTAGCGTACTTGAGGTGTTGTATTTT Started job on | Dec 07 08:52:42 Started mapping on | Dec 07 08:52:42 Finished on | Dec 07 08:52:46 Mapping speed, Million of reads per hour | 1295.69 Number of input reads | 1439651 Average input read length | 92 UNIQUE READS: Uniquely mapped reads number | 1016239 Uniquely mapped reads % | 70.59% Average mapped length | 91.86 Number of splices: Total | 60085 Number of splices: Annotated (sjdb) | 50269 Number of splices: GT/AG | 56825 Number of splices: GC/AG | 1900 Number of splices: AT/AC | 23 Number of splices: Non-canonical | 1337 Mismatch rate per base, % | 0.51% Deletion rate per base | 0.05% Deletion average length | 1.86 Insertion rate per base | 0.02% Insertion average length | 1.77 MULTI-MAPPING READS: Number of reads mapped to multiple loci | 386213 % of reads mapped to multiple loci | 26.83% Number of reads mapped to too many loci | 9594 % of reads mapped to too many loci | 0.67% UNMAPPED READS: % of reads unmapped: too many mismatches | 0.00% % of reads unmapped: too short | 1.86% % of reads unmapped: other | 0.05% CHIMERIC READS: Number of chimeric reads | 0 % of chimeric reads | 0.00% N_unmapped 37199 37199 37199 N_multimapping 386213 386213 386213 N_noFeature 60699 70182 967280 N_ambiguous 43627 4178 123 UnstrandedReadsAssigned:911913 PositiveStrandReadsAssigned:941879 NegativeStrandReadsAssigned:48836 Dataset is classified positive stranded MeadianReadLen=93 20thPercentileLength=93 echo kmer=89 ERR6133570 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: ERR6133570-trimmed.fastq [quant] finding pseudoalignments for the reads ... done [quant] processed 1,439,651 reads, 1,210,335 reads pseudoaligned [ em] quantifying the abundances ... done [ em] the Expectation-Maximization algorithm ran for 1,004 rounds 52973 ERR6133570.ke.tsv 35125 ERR6133570.se.tsv 88098 total ==> ERR6133570.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 24 19.1101 PNS24243 293 194 0 0 KQK14069 1603 1504 6 4.35821 KQK14071 474 375 0 0 ==> ERR6133570.se.tsv <== BRADI_1g14170v3 6 BRADI_1g53295v3 37 BRADI_1g59795v3 7 BRADI_1g07683v3 0 BRADI_1g00485v3 0 BRADI_1g20270v3 33 BRADI_1g74790v3 5 BRADI_1g09890v3 0 BRADI_1g77505v3 22 BRADI_1g48960v3 0 ERR6133570 completed mapping pipeline successfully