Starting /dee2/code/volunteer_pipeline.sh SRR921326 current disk space = 1547750924288 free memory = 1600617700 SRR921326 SRAfilesize e184feea48af1df8df09a8665896dc40 SRR921326.sra SRR921326.sra file validated SRR921326 is single end SRR921326 is conventional basespace SRR921326 read1 length is 35 nt ##FastQC 0.11.5 >>Basic Statistics pass #Measure Value Filename SRR921326_1.fastq File type Conventional base calls Encoding Sanger / Illumina 1.9 Total Sequences 4000 Sequences flagged as poor quality 0 Sequence length 35 %GC 49 >>END_MODULE >>Per base sequence quality pass #Base Mean Median Lower Quartile Upper Quartile 10th Percentile 90th Percentile 1 32.62975 39.0 36.0 40.0 2.0 40.0 2 33.95425 39.0 35.0 40.0 12.0 40.0 3 34.04425 39.0 36.0 40.0 13.0 40.0 4 33.914 39.0 35.0 40.0 12.0 40.0 5 33.9115 39.0 36.0 40.0 12.0 40.0 6 34.796 39.0 36.0 40.0 20.0 40.0 7 34.7165 39.0 36.0 40.0 19.0 40.0 8 34.6525 39.0 35.0 40.0 19.0 40.0 9 34.5505 39.0 35.0 40.0 18.0 40.0 10 34.55575 39.0 35.0 40.0 18.0 40.0 11 36.048 39.0 35.0 40.0 31.0 40.0 12 36.0545 39.0 35.0 40.0 31.0 40.0 13 36.1255 39.0 35.0 40.0 31.0 40.0 14 36.1345 39.0 35.0 40.0 31.0 40.0 15 36.213 39.0 35.0 40.0 31.0 40.0 16 36.157 39.0 35.0 40.0 31.0 40.0 17 36.22825 39.0 35.0 40.0 31.0 40.0 18 36.08425 39.0 35.0 40.0 31.0 40.0 19 35.97175 39.0 35.0 40.0 31.0 40.0 20 35.94275 39.0 34.0 40.0 31.0 40.0 21 35.78525 39.0 35.0 40.0 31.0 40.0 22 35.494 39.0 34.0 40.0 29.0 40.0 23 35.29325 39.0 34.0 40.0 28.0 40.0 24 35.36725 39.0 34.0 40.0 29.0 40.0 25 35.27425 39.0 34.0 40.0 29.0 40.0 26 34.81625 38.0 33.0 40.0 27.0 40.0 27 34.814 38.0 33.0 40.0 27.0 40.0 28 34.69775 38.0 33.0 40.0 27.0 40.0 29 34.4985 38.0 33.0 40.0 27.0 40.0 30 34.35175 38.0 33.0 40.0 26.0 40.0 31 34.2925 38.0 33.0 40.0 25.0 40.0 32 34.12475 38.0 32.0 40.0 25.0 40.0 33 34.0835 38.0 32.0 40.0 25.0 40.0 34 33.78625 38.0 31.0 40.0 25.0 40.0 35 33.6605 38.0 31.0 40.0 24.0 40.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 >>END_MODULE >>Per sequence quality scores pass #Quality Count 2 96.0 3 1.0 4 0.0 5 0.0 6 1.0 7 1.0 8 0.0 9 2.0 10 2.0 11 8.0 12 2.0 13 2.0 14 1.0 15 8.0 16 4.0 17 9.0 18 8.0 19 11.0 20 25.0 21 25.0 22 39.0 23 58.0 24 74.0 25 115.0 26 128.0 27 140.0 28 44.0 29 30.0 30 37.0 31 31.0 32 55.0 33 68.0 34 88.0 35 112.0 36 177.0 37 345.0 38 597.0 39 1633.0 40 23.0 >>END_MODULE >>Per base sequence content fail #Base G A T C 1 22.394282310899342 26.474091721262656 24.95533055390113 26.176295413936867 2 24.975 28.050000000000004 22.6 24.375 3 21.925 20.225 23.35 34.5 4 21.7 24.425 37.525 16.35 5 20.5 24.05 20.325 35.125 6 22.125 35.325 23.849999999999998 18.7 7 21.275 36.85 18.525 23.35 8 21.575 39.425 20.025000000000002 18.975 9 24.6 32.975 22.85 19.575 10 40.725 22.75 17.8 18.725 11 28.125 33.324999999999996 18.95 19.6 12 34.5 29.15 18.9 17.45 13 41.375 24.375 16.45 17.8 14 23.7 41.275 17.375 17.65 15 28.050000000000004 33.275 20.7 17.974999999999998 16 22.825 33.775 22.95 20.45 17 34.725 23.375 23.025000000000002 18.875 18 37.0 21.825 19.650000000000002 21.525 19 29.325000000000003 21.675 19.2 29.799999999999997 20 22.05 27.55 29.225 21.175 21 24.275 22.05 34.025 19.650000000000002 22 36.675000000000004 24.4 16.400000000000002 22.525000000000002 23 24.7 22.175 18.725 34.4 24 33.050000000000004 25.575 24.125 17.25 25 37.925 23.575 19.725 18.775 26 29.075 21.025 31.25 18.65 27 34.225 24.625 20.0 21.15 28 34.1 22.45 22.15 21.3 29 21.05 35.825 19.35 23.775 30 22.175 27.975 30.95 18.9 31 29.15 34.125 17.599999999999998 19.125 32 26.75 22.175 20.3 30.775000000000002 33 24.3 25.6 20.599999999999998 29.5 34 23.674999999999997 22.175 34.825 19.325 35 27.05 35.825 20.349999999999998 16.775000000000002 >>END_MODULE >>Per sequence GC content warn #GC Content Count 0 1.0 1 2.0 2 3.0 3 3.0 4 3.5 5 4.0 6 4.0 7 6.0 8 8.0 9 8.0 10 5.0 11 2.0 12 2.0 13 5.5 14 9.0 15 9.0 16 8.0 17 7.0 18 7.0 19 11.5 20 16.0 21 18.5 22 21.0 23 21.0 24 23.5 25 26.0 26 26.0 27 38.5 28 51.0 29 51.0 30 70.5 31 90.0 32 90.0 33 115.0 34 140.0 35 140.0 36 166.0 37 192.0 38 192.0 39 229.0 40 266.0 41 326.0 42 386.0 43 386.0 44 442.5 45 499.0 46 499.0 47 510.5 48 522.0 49 522.0 50 465.0 51 408.0 52 408.0 53 352.0 54 296.0 55 296.0 56 274.5 57 253.0 58 253.0 59 220.5 60 188.0 61 173.5 62 159.0 63 159.0 64 143.5 65 128.0 66 128.0 67 126.0 68 124.0 69 124.0 70 104.0 71 84.0 72 84.0 73 71.5 74 59.0 75 59.0 76 40.5 77 22.0 78 22.0 79 19.0 80 16.0 81 14.0 82 12.0 83 12.0 84 8.5 85 5.0 86 5.0 87 3.5 88 2.0 89 2.0 90 1.0 91 0.0 92 0.0 93 0.5 94 1.0 95 1.0 96 0.5 97 0.0 98 0.0 99 0.0 100 0.0 >>END_MODULE >>Per base N content warn #Base N-Count 1 16.05 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.0 27 0.0 28 0.0 29 0.0 30 0.0 31 0.0 32 0.0 33 0.0 34 0.0 35 0.0 >>END_MODULE >>Sequence Length Distribution pass #Length Count 35 4000.0 >>END_MODULE >>Sequence Duplication Levels pass #Total Deduplicated Percentage 86.05000000000001 #Duplication Level Percentage of deduplicated Percentage of total 1 95.2643811737362 81.975 2 1.6560139453805929 2.85 3 0.9877977919814062 2.55 4 0.5520046484601976 1.9 5 0.3486345148169669 1.5 6 0.31958163858221966 1.6500000000000001 7 0.17431725740848344 1.05 8 0.08715862870424172 0.6 9 0.2324230098779779 1.7999999999999998 >10 0.3776873910517141 4.125 >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 ATTTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 19 0.475 No Hit CTCTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 16 0.4 No Hit ATGTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 15 0.375 No Hit AAGTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 14 0.35000000000000003 No Hit CTTTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 13 0.325 No Hit CCGTCTCATGGAGAGTTCGATCCTGGCTCAGGATG 13 0.325 No Hit AATTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 13 0.325 No Hit ACCTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 11 0.27499999999999997 No Hit TCTTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 11 0.27499999999999997 No Hit CACTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 10 0.25 No Hit CCATCTCATGGAGAGTTCGATCCTGGCTCAGGATG 10 0.25 No Hit CTGTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 10 0.25 No Hit AAATCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 10 0.25 No Hit TACTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 9 0.22499999999999998 No Hit AGCTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 9 0.22499999999999998 No Hit CTCTCTCATGGAGAGTTCGATCCTGGCTCAGGATG 9 0.22499999999999998 No Hit TGCTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 9 0.22499999999999998 No Hit TTTTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 9 0.22499999999999998 No Hit NATTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 9 0.22499999999999998 No Hit TTCTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 9 0.22499999999999998 No Hit NACTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 9 0.22499999999999998 No Hit GCCTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 8 0.2 No Hit AACTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 8 0.2 No Hit NTGTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 8 0.2 No Hit NGCTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 7 0.17500000000000002 No Hit NTTTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 7 0.17500000000000002 No Hit GCTTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 7 0.17500000000000002 No Hit TTATCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 7 0.17500000000000002 No Hit NTCTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 7 0.17500000000000002 No Hit ATCTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 7 0.17500000000000002 No Hit NCTTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 6 0.15 No Hit ACTTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 6 0.15 No Hit ACATCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 6 0.15 No Hit CCCTCTCATGGAGAGTTCGATCCTGGCTCAGGATG 6 0.15 No Hit ACGTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 6 0.15 No Hit GACTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 6 0.15 No Hit TCGTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 6 0.15 No Hit TCCTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 6 0.15 No Hit TATTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 6 0.15 No Hit CGCTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 6 0.15 No Hit GCATCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 6 0.15 No Hit CATTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 5 0.125 No Hit AGCTACCTGGTTGATCCTGCCAGTAGTCATATGCT 5 0.125 No Hit CAATCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 5 0.125 No Hit AACTCTCATGGAGAGTTCGATCCTGGCTCAGGATG 5 0.125 No Hit TGGTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 5 0.125 No Hit TAATCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 5 0.125 No Hit AGGTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 5 0.125 No Hit CCATCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 5 0.125 No Hit GGCGCGACCCCAGGTCAGGCGGGACTACCCGCTGA 5 0.125 No Hit GGCTCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 5 0.125 No Hit CTATCAAAAGAGGAAAGGCTTGCGGTGGATACCTA 5 0.125 No Hit TGCTACCTGGTTGATCCTGCCAGTAGTCATATGCT 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.025 0.0 0.0 0.0 0.0 2 0.075 0.0 0.0 0.0 0.0 3 0.075 0.0 0.0 0.0 0.0 4 0.125 0.0 0.0 0.0 0.0 5 0.15 0.0 0.0 0.0 0.0 6 0.15 0.0 0.0 0.0 0.0 7 0.15 0.0 0.0 0.0 0.0 8 0.15 0.0 0.0 0.0 0.0 9 0.15 0.0 0.0 0.0 0.0 10 0.15 0.0 0.0 0.0 0.0 11 0.175 0.0 0.0 0.0 0.0 12 0.2 0.0 0.0 0.0 0.0 13 0.25 0.0 0.0 0.0 0.0 14 0.25 0.0 0.0 0.0 0.0 15 0.25 0.0 0.0 0.0 0.0 16 0.275 0.0 0.0 0.0 0.0 17 0.275 0.0 0.0 0.0 0.0 18 0.275 0.0 0.0 0.0 0.0 19 0.275 0.0 0.0 0.0 0.0 20 0.275 0.0 0.0 0.0 0.0 21 0.275 0.0 0.0 0.0 0.0 22 0.275 0.0 0.0 0.0 0.0 23 0.275 0.0 0.0 0.0 0.0 >>END_MODULE >>Kmer Content warn #Sequence Count PValue Obs/Exp Max Max Obs/Exp Position ATCCTGG 25 4.447609E-4 28.8375 20 ATCTCAT 20 0.005538505 28.8375 3 TCGATCC 25 4.447609E-4 28.8375 17 GAGGAAA 20 0.005538505 28.8375 10 GATACCT 20 0.005538505 28.8375 28 TCAGGAT 25 4.447609E-4 28.8375 28 AAGAGGA 20 0.005538505 28.8375 8 GCGGTGG 20 0.005538505 28.8375 22 AGGAAAG 20 0.005538505 28.8375 11 TGGAGAG 25 4.447609E-4 28.8375 9 ATGGAGA 25 4.447609E-4 28.8375 8 CTGGCTC 25 4.447609E-4 28.8375 23 GAGTTCG 25 4.447609E-4 28.8375 13 GGATACC 20 0.005538505 28.8375 27 AGAGTTC 25 4.447609E-4 28.8375 12 TCAAAAG 20 0.005538505 28.8375 4 AAAAGAG 20 0.005538505 28.8375 6 AGGCTTG 20 0.005538505 28.8375 16 GGCTTGC 20 0.005538505 28.8375 17 AAGGCTT 20 0.005538505 28.8375 15 >>END_MODULE Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686960 spots for SRR921326.sra Written 686960 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra Read 686952 spots for SRR921326.sra Written 686952 spots for SRR921326.sra SRR ids: ['SRR921326.sra'] extra args: ['--split-files', '--defline-qual', '+'] tempdir: /tmp/pfd_t3yenjmy SRR921326.sra spots: 13739048 blocks: [[1, 686952], [686953, 1373904], [1373905, 2060856], [2060857, 2747808], [2747809, 3434760], [3434761, 4121712], [4121713, 4808664], [4808665, 5495616], [5495617, 6182568], [6182569, 6869520], [6869521, 7556472], [7556473, 8243424], [8243425, 8930376], [8930377, 9617328], [9617329, 10304280], [10304281, 10991232], [10991233, 11678184], [11678185, 12365136], [12365137, 13052088], [13052089, 13739048]] SRR921326 file size 2003950 SRR921326 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 SRR921326 SRR921326_1.fastq Input file: SRR921326_1.fastq trimmed: SRR921326-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:06:41 2024 >> started Sat Dec 7 00:06:53 2024 >> done (11.857s) 13739048 reads processed; of these: 76000 ( 0.55%) short reads filtered out after trimming by size control 28578 ( 0.21%) empty reads filtered out after trimming by size control 13634470 (99.24%) reads available; of these: 540129 ( 3.96%) trimmed reads available after processing 13094341 (96.04%) untrimmed reads available after processing Length distribution of reads after trimming: length count percentage 18 9315 0.07% 19 16797 0.12% 20 37715 0.28% 21 8213 0.06% 22 14307 0.10% 23 21329 0.16% 24 39531 0.29% 25 96656 0.71% 26 14289 0.10% 27 22232 0.16% 28 36558 0.27% 29 51672 0.38% 30 101342 0.74% 31 12446 0.09% 32 14162 0.10% 33 19895 0.15% 34 23670 0.17% 35 13094341 96.04% 13634470 reads passed initial QC criterion=sequence-density sequence-density=11.65 sequence-density-rank=1 fanout-score=44.88 fanout-score-rank=1 prefix-density=11.44 prefix-fanout=44.9 sequence=TCAAAAGAGGAAAGGCTTGCGGTGGATACCTAGGCACCCAGAGACGAGGAAGGGCGTAGCAAGCGACGAAATGCTTCGGGGAGT criterion=fanout-score sequence-density=11.65 sequence-density-rank=1 fanout-score=44.88 fanout-score-rank=1 prefix-density=11.44 prefix-fanout=44.9 sequence=TCAAAAGAGGAAAGGCTTGCGGTGGATACCTAGGCACCCAGAGACGAGGAAGGGCGTAGCAAGCGACGAAATGCTTCGGGGAGT 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 TCAAAAGAGGAAAGGCTTGCGGTGGATACCTAGGCACCCAGAGACGAGGAAGGGCGTAGCAAGCGACGAAATGCTTCGGGGAGT -o SRR921326 - Input file: STDIN trimmed: SRR921326-trimmed.fastq Parameters used: -- 3' end adapter sequence (-x): TCAAAAGAGGAAAGGCTTGCGGTGGATACCTAGGCACCCAGAGACGAGGAAGGGCGTAGCAAGC -- 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:07:11 2024 >> started Sat Dec 7 00:07:22 2024 >> done (11.119s) 11362058 reads processed; of these: 1309957 (11.53%) short reads filtered out after trimming by size control 3147 ( 0.03%) empty reads filtered out after trimming by size control 10048954 (88.44%) reads available; of these: 203739 ( 2.03%) trimmed reads available after processing 9845215 (97.97%) untrimmed reads available after processing Length distribution of reads after trimming: length count percentage 18 8373 0.08% 19 13762 0.14% 20 29252 0.29% 21 7707 0.08% 22 12339 0.12% 23 16240 0.16% 24 30641 0.30% 25 69560 0.69% 26 11607 0.12% 27 17396 0.17% 28 24633 0.25% 29 40635 0.40% 30 85303 0.85% 31 54979 0.55% 32 142993 1.42% 33 11136 0.11% 34 18239 0.18% 35 9454159 94.08% criterion=sequence-density sequence-density=4.77 sequence-density-rank=1 fanout-score=38.37 fanout-score-rank=1 prefix-density=4.69 prefix-fanout=38.4 sequence=TCTCATGGAGAGTTCGATCCTGGCTCAGGATG criterion=fanout-score sequence-density=4.77 sequence-density-rank=1 fanout-score=38.37 fanout-score-rank=1 prefix-density=4.69 prefix-fanout=38.4 sequence=TCTCATGGAGAGTTCGATCCTGGCTCAGGATG Started job on | Dec 07 00:07:52 Started mapping on | Dec 07 00:07:52 Finished on | Dec 07 00:08:03 Mapping speed, Million of reads per hour | 4032.45 Number of input reads | 12321366 Average input read length | 30 UNIQUE READS: Uniquely mapped reads number | 8421982 Uniquely mapped reads % | 68.35% Average mapped length | 30.46 Number of splices: Total | 677612 Number of splices: Annotated (sjdb) | 643867 Number of splices: GT/AG | 665917 Number of splices: GC/AG | 8416 Number of splices: AT/AC | 387 Number of splices: Non-canonical | 2892 Mismatch rate per base, % | 0.32% Deletion rate per base | 0.00% Deletion average length | 1.36 Insertion rate per base | 0.00% Insertion average length | 1.17 MULTI-MAPPING READS: Number of reads mapped to multiple loci | 1883631 % of reads mapped to multiple loci | 15.29% Number of reads mapped to too many loci | 798347 % of reads mapped to too many loci | 6.48% UNMAPPED READS: % of reads unmapped: too many mismatches | 0.00% % of reads unmapped: too short | 9.48% % of reads unmapped: other | 0.40% CHIMERIC READS: Number of chimeric reads | 0 % of chimeric reads | 0.00% N_unmapped 2015753 2015753 2015753 N_multimapping 1883631 1883631 1883631 N_noFeature 470652 647466 8095629 N_ambiguous 164465 14617 1220 UnstrandedReadsAssigned:7786865 PositiveStrandReadsAssigned:7759899 NegativeStrandReadsAssigned:325133 Dataset is classified positive stranded MeadianReadLen=31 20thPercentileLength=31 echo kmer=27 SRR921326 Starting Kallisto single end mapping to ensembl reference transcriptome. kmer=27 [quant] fragment length distribution is truncated gaussian with mean = 100, sd = 20 [index] k-mer length: 27 [index] number of targets: 52,972 [index] number of k-mers: 66,546,128 [index] number of equivalence classes: 124,662 [quant] running in single-end mode [quant] will process file 1: SRR921326-trimmed.fastq [quant] finding pseudoalignments for the reads ... done [quant] processed 12,321,366 reads, 7,287,366 reads pseudoaligned [ em] quantifying the abundances ... done [ em] the Expectation-Maximization algorithm ran for 1,149 rounds 52973 SRR921326.ke.tsv 35125 SRR921326.se.tsv 88098 total ==> SRR921326.ke.tsv <== target_id length eff_length est_counts tpm PNS24245 936 837 0 0 PNS24247 1044 945 0 0 PNS24249 1928 1829 10.1241 1.15536 PNS24246 1044 945 0 0 PNS24248 1044 945 0 0 PNS24244 1471 1372 96.8759 14.7379 PNS24243 293 194 0 0 KQK14069 1603 1504 3202.86 444.493 KQK14071 474 375 1854.6 1032.27 ==> SRR921326.se.tsv <== BRADI_1g14170v3 5382 BRADI_1g53295v3 29 BRADI_1g59795v3 144 BRADI_1g07683v3 0 BRADI_1g00485v3 1 BRADI_1g20270v3 143 BRADI_1g74790v3 35 BRADI_1g09890v3 5 BRADI_1g77505v3 71 BRADI_1g48960v3 1 SRR921326 completed mapping pipeline successfully