Starting /dee2/code/volunteer_pipeline.sh SRR1033812 current disk space = 1548083634176 free memory = 1596676920 SRR1033812 SRAfilesize 234906d822bf5675ade7e81f7f9b3b0f SRR1033812.sra SRR1033812.sra file validated SRR1033812 is single end SRR1033812 is conventional basespace SRR1033812 read1 length is 36 nt ##FastQC 0.11.5 >>Basic Statistics pass #Measure Value Filename SRR1033812_1.fastq File type Conventional base calls Encoding Sanger / Illumina 1.9 Total Sequences 4000 Sequences flagged as poor quality 0 Sequence length 36 %GC 50 >>END_MODULE >>Per base sequence quality fail #Base Mean Median Lower Quartile Upper Quartile 10th Percentile 90th Percentile 1 32.23675 33.0 33.0 33.0 32.0 33.0 2 31.74325 33.0 32.0 33.0 29.0 34.0 3 30.72875 32.0 30.0 33.0 26.0 33.0 4 31.1195 32.0 31.0 33.0 27.0 34.0 5 31.034 32.0 30.0 33.0 27.0 33.0 6 31.255 33.0 31.0 33.0 27.0 34.0 7 30.866 32.0 30.0 33.0 27.0 33.0 8 30.76675 32.0 30.0 33.0 26.0 33.0 9 30.14775 32.0 29.0 33.0 24.0 33.0 10 30.323 32.0 29.0 33.0 25.0 33.0 11 30.3495 32.0 30.0 33.0 24.0 33.0 12 30.564 32.0 30.0 33.0 26.0 33.0 13 30.691 32.0 30.0 33.0 26.0 33.0 14 29.2685 31.0 28.0 33.0 23.0 33.0 15 29.72525 32.0 29.0 33.0 24.0 33.0 16 29.05325 31.0 28.0 33.0 22.0 33.0 17 28.20375 30.0 26.0 32.0 21.0 33.0 18 28.70125 31.0 27.0 33.0 22.0 33.0 19 28.85125 31.0 28.0 33.0 22.0 33.0 20 28.18275 31.0 27.0 33.0 20.0 33.0 21 28.014 31.0 27.0 33.0 19.0 33.0 22 27.2515 30.0 26.0 32.0 16.0 33.0 23 25.232 28.0 22.0 31.0 8.0 33.0 24 23.8415 27.0 20.0 31.0 4.0 32.0 25 24.2025 27.0 21.0 31.0 4.0 33.0 26 23.05175 26.0 19.0 30.0 4.0 32.0 27 22.4605 25.0 18.0 29.0 4.0 32.0 28 21.18475 24.0 17.0 28.0 4.0 31.0 29 22.962 26.0 20.0 30.0 4.0 32.0 30 23.14925 27.0 20.0 30.0 4.0 32.0 31 21.9065 25.0 17.0 29.0 4.0 32.0 32 20.8075 24.0 15.0 28.0 4.0 31.0 33 20.75525 25.0 12.0 29.0 4.0 31.0 34 19.66325 24.0 4.0 29.0 4.0 32.0 35 15.2425 17.0 4.0 25.0 4.0 29.0 36 14.3065 4.0 4.0 25.0 4.0 31.0 >>END_MODULE >>Per tile sequence quality pass #Tile Base Mean 1 1 0.0 1 2 0.0 1 3 0.0 1 4 0.0 1 5 0.0 1 6 0.0 1 7 0.0 1 8 0.0 1 9 0.0 1 10 0.0 1 11 0.0 1 12 0.0 1 13 0.0 1 14 0.0 1 15 0.0 1 16 0.0 1 17 0.0 1 18 0.0 1 19 0.0 1 20 0.0 1 21 0.0 1 22 0.0 1 23 0.0 1 24 0.0 1 25 0.0 1 26 0.0 1 27 0.0 1 28 0.0 1 29 0.0 1 30 0.0 1 31 0.0 1 32 0.0 1 33 0.0 1 34 0.0 1 35 0.0 1 36 0.0 >>END_MODULE >>Per sequence quality scores pass #Quality Count 4 10.0 5 4.0 6 2.0 7 2.0 8 4.0 9 4.0 10 6.0 11 13.0 12 19.0 13 17.0 14 36.0 15 35.0 16 61.0 17 63.0 18 81.0 19 86.0 20 81.0 21 93.0 22 87.0 23 118.0 24 183.0 25 254.0 26 357.0 27 424.0 28 521.0 29 516.0 30 485.0 31 323.0 32 106.0 33 9.0 >>END_MODULE >>Per base sequence content fail #Base G A T C 1 19.747474747474747 39.77272727272727 29.873737373737374 10.606060606060606 2 37.02777082812109 25.769326995246434 23.167375531648737 14.035526644983737 3 32.59944958719039 29.522141606204656 21.015761821366024 16.86264698523893 4 27.950000000000003 35.199999999999996 16.5 20.349999999999998 5 27.575 31.6 18.575 22.25 6 27.525 26.674999999999997 25.724999999999998 20.075000000000003 7 28.14610958218664 29.29697272954716 27.34550913184889 15.211408556417313 8 21.935967983991997 33.16658329164582 20.51025512756378 24.387193596798397 9 33.033033033033036 26.926926926926924 24.174174174174173 15.865865865865866 10 32.3330832708177 32.808202050512634 15.753938484621155 19.10477619404851 11 29.132283070767688 27.7569392348087 22.155538884721178 20.955238809702426 12 31.874999999999996 30.2 17.825 20.1 13 37.26226226226226 28.653653653653656 14.814814814814813 19.26926926926927 14 27.900000000000002 31.1 24.9 16.1 15 36.459114778694676 25.656414103525883 22.55563890972743 15.328832208052013 16 34.8 30.599999999999998 17.1 17.5 17 34.15061295971979 31.923942957217914 18.388791593695274 15.536652489367025 18 34.33358339584896 26.70667666916729 16.27906976744186 22.680670167541887 19 27.85696424106027 28.657164291072768 20.7551887971993 22.73068267066767 20 22.275 28.249999999999996 27.0 22.475 21 22.775000000000002 23.775 29.975 23.474999999999998 22 21.875 17.7 34.625 25.8 23 26.700000000000003 21.825 18.8 32.675 24 27.956989247311824 16.27906976744186 36.00900225056264 19.754938734683673 25 9.00225056264066 11.177794448612154 66.5666416604151 13.253313328332084 26 11.15 18.6 22.875 47.375 27 56.38909727431858 7.226806701675419 29.957489372343087 6.426606651662915 28 22.7977977977978 10.46046046046046 54.87987987987988 11.861861861861863 29 9.632224168126093 47.43557668251188 15.836877658243683 27.09532149111834 30 19.50975487743872 5.52776388194097 50.40020010005003 24.562281140570285 31 64.95747873936969 2.4762381190595297 15.207603801900952 17.358679339669834 32 12.609457092819614 0.5253940455341506 21.1408556417313 65.72429321991493 33 12.887887887887889 0.47547547547547553 17.117117117117118 69.51951951951952 34 47.873936968484244 0.22511255627813906 36.84342171085543 15.057528764382191 35 6.426606651662915 0.27506876719179796 71.51787946986747 21.780445111277817 36 2.9007251812953236 0.3750937734433608 37.259314828707176 59.464866216554135 >>END_MODULE >>Per sequence GC content warn #GC Content Count 0 1.0 1 0.5 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 0.0 18 0.0 19 0.0 20 0.0 21 0.0 22 0.0 23 0.0 24 0.0 25 0.0 26 0.5 27 1.0 28 1.0 29 4.0 30 7.0 31 7.0 32 12.0 33 17.0 34 17.0 35 24.5 36 32.0 37 32.0 38 60.5 39 89.0 40 157.0 41 225.0 42 225.0 43 293.0 44 361.0 45 361.0 46 436.5 47 512.0 48 512.0 49 774.5 50 1037.0 51 853.0 52 669.0 53 669.0 54 575.0 55 481.0 56 481.0 57 385.5 58 290.0 59 290.0 60 222.0 61 154.0 62 154.0 63 107.5 64 61.0 65 51.5 66 42.0 67 42.0 68 27.0 69 12.0 70 12.0 71 8.5 72 5.0 73 5.0 74 3.5 75 2.0 76 1.5 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.5 91 1.0 92 1.0 93 0.5 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 1.0 2 0.075 3 0.075 4 0.0 5 0.0 6 0.0 7 0.075 8 0.05 9 0.1 10 0.025 11 0.025 12 0.0 13 0.1 14 0.0 15 0.025 16 0.0 17 0.075 18 0.025 19 0.025 20 0.0 21 0.0 22 0.0 23 0.0 24 0.025 25 0.025 26 0.0 27 0.025 28 0.1 29 0.075 30 0.05 31 0.05 32 0.075 33 0.1 34 0.05 35 0.025 36 0.025 >>END_MODULE >>Sequence Length Distribution pass #Length Count 36 4000.0 >>END_MODULE >>Sequence Duplication Levels pass #Total Deduplicated Percentage 76.1 #Duplication Level Percentage of deduplicated Percentage of total 1 95.2365308804205 72.475 2 2.726675427069645 4.15 3 0.821287779237845 1.875 4 0.22996057818659657 0.7000000000000001 5 0.19710906701708278 0.75 6 0.09855453350854139 0.44999999999999996 7 0.0328515111695138 0.17500000000000002 8 0.0328515111695138 0.2 9 0.09855453350854139 0.675 >10 0.42706964520367935 6.5 >50 0.0328515111695138 1.95 >100 0.0657030223390276 10.100000000000001 >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 TGACAGAAGAGAGTGAGCACTCGTATGCCGTCTTCT 232 5.800000000000001 No Hit TTGACAGAAGAGAGTGAGCACTCGTATGCCGTCTTC 172 4.3 No Hit GGTAGTTCGACCGCGGAATTTTCGTATGCCGTCTTC 78 1.95 No Hit CACGACTCTCGGCAACGGATATCTTCGTATGCCGTC 43 1.075 No Hit TCGCTTGGTGCAGATCGGGACTCGTATGCCGTCTTC 38 0.95 No Hit AGAATCTTGATGATGCTGCATTCGTATGCCGTCTTC 29 0.7250000000000001 No Hit TGAAGCTGCCAGCATGATCTGATCGTATGCCGTCTT 22 0.5499999999999999 No Hit GGTAGTTCGACCGCGGAATTTTCGTATGCCGTCTTT 20 0.5 No Hit GACACGACTCTCGGCAACGGATATTCGTATGCCGTC 16 0.4 No Hit CACGACTCTCGGCAACGGATATCTCTCGTATGCCGT 16 0.4 No Hit GGTAGTTCGACCGCGGAATTTATCGTATGCCGTCTT 15 0.375 No Hit GGTAGTTCGACCGCGGAATTTCGTATGCCGTCTTCT 14 0.35000000000000003 Illumina PCR Primer Index 8 (95% over 21bp) CACGACTCTCGGCAACGGATATCTCGTATGCCGTCT 13 0.325 No Hit GACACGACTCTCGGCAACGGATATCTTCGTATGCCG 13 0.325 No Hit GGTAGTTCGACCGCGGAATTTTTCGTATGCCGTCTT 11 0.27499999999999997 No Hit CTGACAGAAGAGAGTGAGCACTCGTATGCCGTCTTC 10 0.25 No Hit CACGACTCTCGGCAACGGATATCTCGTCGTATGCCG 9 0.22499999999999998 No Hit TTGACAGAAGAGAGTGAGCACTCGTATGCCGTCTTT 9 0.22499999999999998 No Hit CACGACTCTCGGCAACGGATATTCGTATGCCGTCTT 9 0.22499999999999998 No Hit GACACGACTCTCGGCAACGGATATCTCGTATGCCGT 8 0.2 No Hit GGGGATGTAGCTCAGATGGTAGATCGTATGCCGTCT 7 0.17500000000000002 No Hit ACGACTCTCGGCAACGGATATCTTCGTATGCCGTCT 6 0.15 No Hit GGGGATGTAGCTCAGATGGTAGAGTCGTATGCCGTC 6 0.15 No Hit GGTAGTTCGACCGCGGAATTATCGTATGCCGTCTTC 6 0.15 No Hit GGGGATGTAGCTCAGATGGTTCGTATGCCGTCTTCT 5 0.125 No Hit AGCTCTGATACCATGTGGATGAGATCGTATGCCGTC 5 0.125 No Hit GGGGATGTAGCTCAGATGGTATCGTATGCCGTCTTC 5 0.125 No Hit GACACGACTCTCGGCAACGGATTCGTATGCCGTCTT 5 0.125 No Hit TCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAAAAA 5 0.125 Illumina Single End Adapter 1 (100% over 21bp) TGGACCTGTCAGAAGAAGGCTTCGTATGCCGTCTTC 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 0.0 0.0 0.0 0.0 0.0 11 0.0 0.0 0.0 0.0 0.0 12 0.0 0.0 0.0 0.0 0.0 13 0.0 0.0 0.0 0.0 0.0 14 0.0 0.0 0.0 0.0 0.0 15 0.0 0.0 0.0 0.0 0.0 16 0.0 0.0 0.0 0.0 0.0 17 0.0 0.0 0.0 0.0 0.0 18 0.0 0.0 0.0 0.0 0.0 19 0.0 0.0 0.0 0.0 0.0 20 0.0 0.0 0.0 0.0 0.0 21 0.0 0.0 0.0 0.0 0.0 22 0.0 0.0 0.0 0.0 0.0 23 0.0 0.0 0.0 0.0 0.0 24 0.0 0.0 0.0 0.0 0.0 >>END_MODULE >>Kmer Content warn #Sequence Count PValue Obs/Exp Max Max Obs/Exp Position TAGTTCG 20 0.0043647434 30.265823 3 GGTAGTT 25 3.341556E-4 30.265823 1 GTAGTTC 20 0.0043647434 30.265823 2 TCGACCG 20 0.0046492554 29.8875 7 TTCGACC 20 0.0046492554 29.8875 6 ACCGCGG 25 3.6060356E-4 29.8875 10 GACCGCG 25 3.6060356E-4 29.8875 9 GGCTCGT 25 3.6060356E-4 29.8875 22 CGCGGAA 25 3.6060356E-4 29.8875 12 GGAATTT 25 3.6060356E-4 29.8875 15 AGTTCGA 20 0.0046492554 29.8875 4 CGACCGC 20 0.0046492554 29.8875 8 GTTCGAC 20 0.0046492554 29.8875 5 CCGCGGA 25 3.6060356E-4 29.8875 11 GCTCGTA 40 6.0214716E-6 26.151564 23 ATTCGTA 35 7.9584715E-5 25.617857 23 CGTCTTC 30 9.720369E-4 25.22152 30 GCGGAAT 30 0.001048598 24.90625 13 CGGAATT 30 0.001048598 24.90625 14 TTCGTAT 80 5.2241376E-9 20.547657 24 >>END_MODULE Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464235 spots for SRR1033812.sra Written 464235 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra Read 464231 spots for SRR1033812.sra Written 464231 spots for SRR1033812.sra SRR ids: ['SRR1033812.sra'] extra args: ['--split-files', '--defline-qual', '+'] tempdir: /tmp/pfd_yxexok_l SRR1033812.sra spots: 9284624 blocks: [[1, 464231], [464232, 928462], [928463, 1392693], [1392694, 1856924], [1856925, 2321155], [2321156, 2785386], [2785387, 3249617], [3249618, 3713848], [3713849, 4178079], [4178080, 4642310], [4642311, 5106541], [5106542, 5570772], [5570773, 6035003], [6035004, 6499234], [6499235, 6963465], [6963466, 7427696], [7427697, 7891927], [7891928, 8356158], [8356159, 8820389], [8820390, 9284624]] SRR1033812 file size 1202788 SRR1033812 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 SRR1033812 SRR1033812_1.fastq Input file: SRR1033812_1.fastq trimmed: SRR1033812-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 00:59:21 2024 >> started Sat Dec 7 00:59:26 2024 >> done (4.268s) 9284624 reads processed; of these: 205420 ( 2.21%) short reads filtered out after trimming by size control 34332 ( 0.37%) empty reads filtered out after trimming by size control 9044872 (97.42%) reads available; of these: 4828903 (53.39%) trimmed reads available after processing 4215969 (46.61%) untrimmed reads available after processing Length distribution of reads after trimming: length count percentage 18 75990 0.84% 19 99382 1.10% 20 134855 1.49% 21 154533 1.71% 22 112296 1.24% 23 106549 1.18% 24 106910 1.18% 25 101329 1.12% 26 81800 0.90% 27 110813 1.23% 28 88791 0.98% 29 107040 1.18% 30 200505 2.22% 31 312173 3.45% 32 336625 3.72% 33 551643 6.10% 34 1187618 13.13% 35 960051 10.61% 36 4215969 46.61% 9044872 reads passed initial QC criterion=sequence-density sequence-density=34.75 sequence-density-rank=1 fanout-score=44.17 fanout-score-rank=3 prefix-density=68.36 prefix-fanout=22.5 sequence=TCGTATGCCGTCTTCTGCTTGAAAA criterion=fanout-score sequence-density=1.61 sequence-density-rank=4 fanout-score=76.53 fanout-score-rank=1 prefix-density=2.62 prefix-fanout=47.1 sequence=GCTCGTATGCCG Potential 3prime adapter identified. Now checking if in reference sequence Warning: gzbuffer added in zlib v1.2.3.5. Unable to change buffer size from default of 8192. 1 reads; of these: 1 (100.00%) were unpaired; of these: 1 (100.00%) aligned 0 times 0 (0.00%) aligned exactly 1 time 0 (0.00%) aligned >1 times 0.00% overall alignment rate Adapter seq not found in reference. Now shuffling file before clipping skewer v0.2.2 [April 4, 2016] COMMAND LINE: skewer -f sanger -l 18 -t 20 -x TCGTATGCCGTCTTCTGCTTGAAAA -o SRR1033812 - Input file: STDIN trimmed: SRR1033812-trimmed.fastq Parameters used: -- 3' end adapter sequence (-x): TCGTATGCCGTCTTCTGCTTGAAAA -- maximum error ratio allowed (-r): 0.100 -- maximum indel error ratio allowed (-d): 0.030 -- minimum read length allowed after trimming (-l): 18 -- file format (-f): Sanger/Illumina 1.8+ FASTQ -- minimum overlap length for adapter detection (-k): 3 -- number of concurrent threads (-t): 20 Sat Dec 7 00:59:41 2024 >> started Sat Dec 7 00:59:51 2024 >> done (9.875s) 8528022 reads processed; of these: 12222 ( 0.14%) short reads filtered out after trimming by size control 5342 ( 0.06%) empty reads filtered out after trimming by size control 8510458 (99.79%) reads available; of these: 7614769 (89.48%) trimmed reads available after processing 895689 (10.52%) untrimmed reads available after processing Length distribution of reads after trimming: length count percentage 18 76741 0.90% 19 133832 1.57% 20 875711 10.29% 21 1593510 18.72% 22 680492 8.00% 23 838055 9.85% 24 3574712 42.00% 25 474416 5.57% 26 196562 2.31% 27 26456 0.31% 28 5219 0.06% 29 2156 0.03% 30 1720 0.02% 31 1916 0.02% 32 1613 0.02% 33 1877 0.02% 34 2595 0.03% 35 2803 0.03% 36 20072 0.24% criterion=sequence-density sequence-density=1.09 sequence-density-rank=1 fanout-score=0.00 fanout-score-rank=12 prefix-density=0.00 prefix-fanout=1.0 sequence=TCGCTTGGTGCAGATCGGGACTC criterion=fanout-score sequence-density=0.01 sequence-density-rank=33 fanout-score=5.97 fanout-score-rank=1 prefix-density=0.04 prefix-fanout=1.0 sequence=GGCTCGGGGCGTGGACTGTTGTCGGCCGTGCT Started job on | Dec 07 01:00:05 Started mapping on | Dec 07 01:00:05 Finished on | Dec 07 01:00:37 Mapping speed, Million of reads per hour | 1015.57 Number of input reads | 9027308 Average input read length | 23 UNIQUE READS: Uniquely mapped reads number | 3718179 Uniquely mapped reads % | 41.19% Average mapped length | 23.20 Number of splices: Total | 53134 Number of splices: Annotated (sjdb) | 3736 Number of splices: GT/AG | 51519 Number of splices: GC/AG | 987 Number of splices: AT/AC | 10 Number of splices: Non-canonical | 618 Mismatch rate per base, % | 0.37% Deletion rate per base | 0.00% Deletion average length | 1.08 Insertion rate per base | 0.00% Insertion average length | 1.02 MULTI-MAPPING READS: Number of reads mapped to multiple loci | 3172953 % of reads mapped to multiple loci | 35.15% Number of reads mapped to too many loci | 1721761 % of reads mapped to too many loci | 19.07% UNMAPPED READS: % of reads unmapped: too many mismatches | 0.00% % of reads unmapped: too short | 2.97% % of reads unmapped: other | 1.62% CHIMERIC READS: Number of chimeric reads | 0 % of chimeric reads | 0.00% N_unmapped 2136176 2136176 2136176 N_multimapping 3172953 3172953 3172953 N_noFeature 3101178 3298623 3512631 N_ambiguous 12265 3158 1211 UnstrandedReadsAssigned:604736 PositiveStrandReadsAssigned:416398 NegativeStrandReadsAssigned:204337 Dataset is classified unstranded MeadianReadLen=24 20thPercentileLength=21 echo kmer=19 SRR1033812 Starting Kallisto single end mapping to ensembl reference transcriptome. kmer=19 [quant] fragment length distribution is truncated gaussian with mean = 100, sd = 20 [index] k-mer length: 19 [index] number of targets: 52,972 [index] number of k-mers: 65,492,969 [index] number of equivalence classes: 320,172 [quant] running in single-end mode [quant] will process file 1: SRR1033812-trimmed.fastq [quant] finding pseudoalignments for the reads ... done [quant] processed 9,027,308 reads, 2,972,216 reads pseudoaligned [ em] quantifying the abundances ... done [ em] the Expectation-Maximization algorithm ran for 1,083 rounds 52973 SRR1033812.ke.tsv 35125 SRR1033812.se.tsv 88098 total ==> SRR1033812.ke.tsv <== target_id length eff_length est_counts tpm PNS24245 936 837 0 0 PNS24247 1044 945 0 0 PNS24249 1928 1829 3.99005 0.125292 PNS24246 1044 945 0 0 PNS24248 1044 945 0 0 PNS24244 1471 1372 0.00994766 0.000416413 PNS24243 293 194 0 0 KQK14069 1603 1504 2.3134 0.0883405 KQK14071 474 375 3.94681 0.604467 ==> SRR1033812.se.tsv <== BRADI_1g14170v3 15 BRADI_1g53295v3 2 BRADI_1g59795v3 6 BRADI_1g07683v3 1 BRADI_1g00485v3 1 BRADI_1g20270v3 2 BRADI_1g74790v3 111 BRADI_1g09890v3 1 BRADI_1g77505v3 1 BRADI_1g48960v3 1 SRR1033812 completed mapping pipeline successfully