ShortStack version 3.8.2 Fri Apr 7 09:38:53 MSK 2017 hostname: ubuntu150 working directory: /hdd/sharing/bobrrep/ShortStack-3.8.2 Settings: RNAfold_version: 2 bowtie_cores: 32 bowtie_m: 50 dicermax: 24 dicermin: 20 error_logfile: ShortStack_1491547133/ErrorLogs.txt foldsize: 300 genomefile: ../ref/ivan_mirbase.fasta logfile: ShortStack_1491547133/Log.txt mincov: 0.5rpm mismatches: 0 mmap: u outdir: ShortStack_1491547133 pad: 75 ranmax: 3 readfile: ../files/first/1-tissue.fastq.gz.trimmed.fastq ../files/first/2-tissue.fastq.gz.trimmed.fastq ../files/first/3-tissue.fastq.gz.trimmed.fastq ../files/first/4-tissue.fastq.gz.trimmed.fastq ../files/first/5-tissue.fastq.gz.trimmed.fastq ../files/first/6-tissue.fastq.gz.trimmed.fastq ../files/first/7-tissue.fastq.gz.trimmed.fastq ../files/first/8-tissue.fastq.gz.trimmed.fastq ../files/first/9-tissue.fastq.gz.trimmed.fastq ../files/first/10-tissue.fastq.gz.trimmed.fastq ../files/first/11-tissue.fastq.gz.trimmed.fastq ../files/first/12-tissue.fastq.gz.trimmed.fastq ../files/first/13-tissue.fastq.gz.trimmed.fastq ../files/first/14-tissue.fastq.gz.trimmed.fastq ../files/first/15-tissue.fastq.gz.trimmed.fastq ../files/first/16-plasma.fastq.gz.trimmed.fastq ../files/first/17-plasma.fastq.gz.trimmed.fastq ../files/first/18-plasma.fastq.gz.trimmed.fastq ../files/first/19-plasma.fastq.gz.trimmed.fastq ../files/first/20-plasma.fastq.gz.trimmed.fastq ../files/first/21-plasma.fastq.gz.trimmed.fastq ../files/first/22-plasma.fastq.gz.trimmed.fastq ../files/first/23-plasma.fastq.gz.trimmed.fastq ../files/first/24-plasma.fastq.gz.trimmed.fastq sort_mem: 10G strand_cutoff: 0.8 Run Progress and Messages: Fri Apr 7 09:38:53 MSK 2017 Genome has more than 50 segments and more than two of them are < 1Mb. Stitching short references to improve performance ... Done Fri Apr 7 09:38:53 MSK 2017 Starting alignment of ../files/first/1-tissue.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:1-tissue.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/1-tissue.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/1-tissue.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 449607 / 3911763 (11.5 %) Multi mappers: 20539 / 3911763 (0.5 %) Non mappers: 3441617 / 3911763 (88.0 %) Fri Apr 7 09:39:45 MSK 2017 Starting alignment of ../files/first/2-tissue.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:2-tissue.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/2-tissue.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/2-tissue.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 2549198 / 12880889 (19.8 %) Multi mappers: 81287 / 12880889 (0.6 %) Non mappers: 10250404 / 12880889 (79.6 %) Fri Apr 7 09:43:06 MSK 2017 Starting alignment of ../files/first/3-tissue.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:3-tissue.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/3-tissue.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/3-tissue.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 467577 / 6113913 (7.6 %) Multi mappers: 13149 / 6113913 (0.2 %) Non mappers: 5633187 / 6113913 (92.1 %) Fri Apr 7 09:44:31 MSK 2017 Starting alignment of ../files/first/4-tissue.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:4-tissue.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/4-tissue.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/4-tissue.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 1441170 / 7617178 (18.9 %) Multi mappers: 41059 / 7617178 (0.5 %) Non mappers: 6134949 / 7617178 (80.5 %) Fri Apr 7 09:46:19 MSK 2017 Starting alignment of ../files/first/5-tissue.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:5-tissue.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/5-tissue.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/5-tissue.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 431242 / 3807196 (11.3 %) Multi mappers: 22078 / 3807196 (0.6 %) Non mappers: 3353876 / 3807196 (88.1 %) Fri Apr 7 09:47:13 MSK 2017 Starting alignment of ../files/first/6-tissue.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:6-tissue.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/6-tissue.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/6-tissue.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 997925 / 6779930 (14.7 %) Multi mappers: 25951 / 6779930 (0.4 %) Non mappers: 5756054 / 6779930 (84.9 %) Fri Apr 7 09:48:49 MSK 2017 Starting alignment of ../files/first/7-tissue.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:7-tissue.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/7-tissue.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/7-tissue.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 592736 / 6083336 (9.7 %) Multi mappers: 29233 / 6083336 (0.5 %) Non mappers: 5461367 / 6083336 (89.8 %) Fri Apr 7 09:50:17 MSK 2017 Starting alignment of ../files/first/8-tissue.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:8-tissue.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/8-tissue.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/8-tissue.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 1521933 / 10173673 (15.0 %) Multi mappers: 43644 / 10173673 (0.4 %) Non mappers: 8608096 / 10173673 (84.6 %) Fri Apr 7 09:52:43 MSK 2017 Starting alignment of ../files/first/9-tissue.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:9-tissue.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/9-tissue.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/9-tissue.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 444655 / 9342031 (4.8 %) Multi mappers: 20364 / 9342031 (0.2 %) Non mappers: 8877012 / 9342031 (95.0 %) Fri Apr 7 09:54:54 MSK 2017 Starting alignment of ../files/first/10-tissue.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:10-tissue.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/10-tissue.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/10-tissue.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 942856 / 8487850 (11.1 %) Multi mappers: 27952 / 8487850 (0.3 %) Non mappers: 7517042 / 8487850 (88.6 %) Fri Apr 7 09:56:49 MSK 2017 Starting alignment of ../files/first/11-tissue.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:11-tissue.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/11-tissue.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/11-tissue.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 725903 / 6938652 (10.5 %) Multi mappers: 35012 / 6938652 (0.5 %) Non mappers: 6177737 / 6938652 (89.0 %) Fri Apr 7 09:58:21 MSK 2017 Starting alignment of ../files/first/12-tissue.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:12-tissue.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/12-tissue.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/12-tissue.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 958614 / 8056108 (11.9 %) Multi mappers: 48782 / 8056108 (0.6 %) Non mappers: 7048712 / 8056108 (87.5 %) Fri Apr 7 10:00:22 MSK 2017 Starting alignment of ../files/first/13-tissue.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:13-tissue.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/13-tissue.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/13-tissue.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 653436 / 7374479 (8.9 %) Multi mappers: 38930 / 7374479 (0.5 %) Non mappers: 6682113 / 7374479 (90.6 %) Fri Apr 7 10:02:07 MSK 2017 Starting alignment of ../files/first/14-tissue.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:14-tissue.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/14-tissue.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/14-tissue.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 628905 / 6391365 (9.8 %) Multi mappers: 30326 / 6391365 (0.5 %) Non mappers: 5732134 / 6391365 (89.7 %) Fri Apr 7 10:03:40 MSK 2017 Starting alignment of ../files/first/15-tissue.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:15-tissue.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/15-tissue.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/15-tissue.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 308958 / 2775661 (11.1 %) Multi mappers: 14396 / 2775661 (0.5 %) Non mappers: 2452307 / 2775661 (88.4 %) Fri Apr 7 10:04:17 MSK 2017 Starting alignment of ../files/first/16-plasma.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:16-plasma.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/16-plasma.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/16-plasma.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 60685 / 6242556 (1.0 %) Multi mappers: 32733 / 6242556 (0.5 %) Non mappers: 6149138 / 6242556 (98.5 %) Fri Apr 7 10:05:47 MSK 2017 Starting alignment of ../files/first/17-plasma.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:17-plasma.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/17-plasma.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/17-plasma.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 21844 / 8939196 (0.2 %) Multi mappers: 15043 / 8939196 (0.2 %) Non mappers: 8902309 / 8939196 (99.6 %) Fri Apr 7 10:07:55 MSK 2017 Starting alignment of ../files/first/18-plasma.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:18-plasma.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/18-plasma.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/18-plasma.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 28289 / 5518851 (0.5 %) Multi mappers: 17939 / 5518851 (0.3 %) Non mappers: 5472623 / 5518851 (99.2 %) Fri Apr 7 10:09:16 MSK 2017 Starting alignment of ../files/first/19-plasma.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:19-plasma.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/19-plasma.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/19-plasma.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 26888 / 4147978 (0.6 %) Multi mappers: 16924 / 4147978 (0.4 %) Non mappers: 4104166 / 4147978 (98.9 %) Fri Apr 7 10:10:12 MSK 2017 Starting alignment of ../files/first/20-plasma.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:20-plasma.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/20-plasma.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/20-plasma.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 75313 / 8344022 (0.9 %) Multi mappers: 48509 / 8344022 (0.6 %) Non mappers: 8220200 / 8344022 (98.5 %) Fri Apr 7 10:12:10 MSK 2017 Starting alignment of ../files/first/21-plasma.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:21-plasma.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/21-plasma.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/21-plasma.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 166028 / 6166801 (2.7 %) Multi mappers: 123947 / 6166801 (2.0 %) Non mappers: 5876826 / 6166801 (95.3 %) Fri Apr 7 10:13:39 MSK 2017 Starting alignment of ../files/first/22-plasma.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:22-plasma.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/22-plasma.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/22-plasma.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 56545 / 5144615 (1.1 %) Multi mappers: 25282 / 5144615 (0.5 %) Non mappers: 5062788 / 5144615 (98.4 %) Fri Apr 7 10:14:50 MSK 2017 Starting alignment of ../files/first/23-plasma.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:23-plasma.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/23-plasma.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/23-plasma.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 156338 / 4133517 (3.8 %) Multi mappers: 89802 / 4133517 (2.2 %) Non mappers: 3887377 / 4133517 (94.0 %) Fri Apr 7 10:15:44 MSK 2017 Starting alignment of ../files/first/24-plasma.fastq.gz.trimmed.fastq with bowtie command: bowtie -q -v 0 -p 32 -S -a -m 50 --sam-RG ID:24-plasma.fastq.gz.trimmed ../ref/ivan_mirbase - < ../files/first/24-plasma.fastq.gz.trimmed.fastq Completed. Results are in temporary file ShortStack_1491547133/24-plasma.fastq.gz.trimmed_readsorted.sam.gz pending final processing Unique mappers: 233147 / 6038891 (3.9 %) Multi mappers: 175422 / 6038891 (2.9 %) Non mappers: 5630322 / 6038891 (93.2 %) Fri Apr 7 10:17:08 MSK 2017 Processing and sorting alignments Working on ShortStack_1491547133/1-tissue.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/2-tissue.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/3-tissue.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/4-tissue.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/5-tissue.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/6-tissue.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/7-tissue.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/8-tissue.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/9-tissue.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/10-tissue.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/11-tissue.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/12-tissue.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/13-tissue.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/14-tissue.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/15-tissue.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/16-plasma.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/17-plasma.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/18-plasma.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/19-plasma.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/20-plasma.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/21-plasma.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/22-plasma.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/23-plasma.fastq.gz.trimmed_readsorted.sam.gz Working on ShortStack_1491547133/24-plasma.fastq.gz.trimmed_readsorted.sam.gz Summary of primary alignments: XY:Z:N -- Unmapped because no valid alignments: 146432356 / 161410451 (90.7 %) XY:Z:M -- Unmapped because alignment number exceeded option bowtie_m 50: 0 / 161410451 (0.0 %) XY:Z:O -- Unmapped because alignment number exceeded option ranmax 3 and no guidance was possible: 6 / 161410451 (0.0 %) XY:Z:U -- Uniquely mapped: 13939792 / 161410451 (8.6 %) XY:Z:R -- Multi-mapped with primary alignment chosen randomly: 446353 / 161410451 (0.3 %) XY:Z:P -- Multi-mapped with primary alignment chosen based on u: 591944 / 161410451 (0.4 %) Fri Apr 7 10:53:38 MSK 2017 Merging files using the following samtools merge command: samtools merge ShortStack_1491547133/merged_alignments.bam ShortStack_1491547133/1-tissue.fastq.gz.trimmed.bam ShortStack_1491547133/2-tissue.fastq.gz.trimmed.bam ShortStack_1491547133/3-tissue.fastq.gz.trimmed.bam ShortStack_1491547133/4-tissue.fastq.gz.trimmed.bam ShortStack_1491547133/5-tissue.fastq.gz.trimmed.bam ShortStack_1491547133/6-tissue.fastq.gz.trimmed.bam ShortStack_1491547133/7-tissue.fastq.gz.trimmed.bam ShortStack_1491547133/8-tissue.fastq.gz.trimmed.bam ShortStack_1491547133/9-tissue.fastq.gz.trimmed.bam ShortStack_1491547133/10-tissue.fastq.gz.trimmed.bam ShortStack_1491547133/11-tissue.fastq.gz.trimmed.bam ShortStack_1491547133/12-tissue.fastq.gz.trimmed.bam ShortStack_1491547133/13-tissue.fastq.gz.trimmed.bam ShortStack_1491547133/14-tissue.fastq.gz.trimmed.bam ShortStack_1491547133/15-tissue.fastq.gz.trimmed.bam ShortStack_1491547133/16-plasma.fastq.gz.trimmed.bam ShortStack_1491547133/17-plasma.fastq.gz.trimmed.bam ShortStack_1491547133/18-plasma.fastq.gz.trimmed.bam ShortStack_1491547133/19-plasma.fastq.gz.trimmed.bam ShortStack_1491547133/20-plasma.fastq.gz.trimmed.bam ShortStack_1491547133/21-plasma.fastq.gz.trimmed.bam ShortStack_1491547133/22-plasma.fastq.gz.trimmed.bam ShortStack_1491547133/23-plasma.fastq.gz.trimmed.bam ShortStack_1491547133/24-plasma.fastq.gz.trimmed.bam Alignment completed. Final bam file is ShortStack_1491547133/merged_alignments.bam Fri Apr 7 11:01:22 MSK 2017 Performing de-novo cluster identification and analyses. At specified mincov of 0.5rpm with 161410451 mapped reads, mincov is 81 raw alignments Completed at Fri Apr 7 11:07:34 MSK 2017 Performing search for unplaced small RNAs. Completed at Fri Apr 7 11:43:35 MSK 2017 Fri Apr 7 11:43:36 MSK 2017 Tally of loci by predominant RNA size (DicerCall): DicerCall NotMIRNA MIRNA N or NA 61 0 20 25 0 21 98 0 22 187 0 23 53 0 24 5 0 Unplaced small RNAs Size NoAlignments <20 39892 20 8005 21 8661 22 10540 23 9026 24 6548 >24 18977