minimap2 0.2.22.0 → 0.2.24.1

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (101) hide show
  1. checksums.yaml +4 -4
  2. data/README.md +60 -76
  3. data/ext/Rakefile +55 -0
  4. data/ext/cmappy/cmappy.c +129 -0
  5. data/ext/cmappy/cmappy.h +44 -0
  6. data/ext/minimap2/FAQ.md +46 -0
  7. data/ext/minimap2/LICENSE.txt +24 -0
  8. data/ext/minimap2/MANIFEST.in +10 -0
  9. data/ext/minimap2/Makefile +132 -0
  10. data/ext/minimap2/Makefile.simde +97 -0
  11. data/ext/minimap2/NEWS.md +821 -0
  12. data/ext/minimap2/README.md +403 -0
  13. data/ext/minimap2/align.c +1020 -0
  14. data/ext/minimap2/bseq.c +169 -0
  15. data/ext/minimap2/bseq.h +64 -0
  16. data/ext/minimap2/code_of_conduct.md +30 -0
  17. data/ext/minimap2/cookbook.md +243 -0
  18. data/ext/minimap2/esterr.c +64 -0
  19. data/ext/minimap2/example.c +63 -0
  20. data/ext/minimap2/format.c +559 -0
  21. data/ext/minimap2/hit.c +466 -0
  22. data/ext/minimap2/index.c +775 -0
  23. data/ext/minimap2/kalloc.c +205 -0
  24. data/ext/minimap2/kalloc.h +76 -0
  25. data/ext/minimap2/kdq.h +132 -0
  26. data/ext/minimap2/ketopt.h +120 -0
  27. data/ext/minimap2/khash.h +615 -0
  28. data/ext/minimap2/krmq.h +474 -0
  29. data/ext/minimap2/kseq.h +256 -0
  30. data/ext/minimap2/ksort.h +153 -0
  31. data/ext/minimap2/ksw2.h +184 -0
  32. data/ext/minimap2/ksw2_dispatch.c +96 -0
  33. data/ext/minimap2/ksw2_extd2_sse.c +402 -0
  34. data/ext/minimap2/ksw2_exts2_sse.c +416 -0
  35. data/ext/minimap2/ksw2_extz2_sse.c +313 -0
  36. data/ext/minimap2/ksw2_ll_sse.c +152 -0
  37. data/ext/minimap2/kthread.c +159 -0
  38. data/ext/minimap2/kthread.h +15 -0
  39. data/ext/minimap2/kvec.h +105 -0
  40. data/ext/minimap2/lchain.c +369 -0
  41. data/ext/minimap2/main.c +459 -0
  42. data/ext/minimap2/map.c +714 -0
  43. data/ext/minimap2/minimap.h +410 -0
  44. data/ext/minimap2/minimap2.1 +725 -0
  45. data/ext/minimap2/misc/README.md +179 -0
  46. data/ext/minimap2/misc/mmphase.js +335 -0
  47. data/ext/minimap2/misc/paftools.js +3149 -0
  48. data/ext/minimap2/misc.c +162 -0
  49. data/ext/minimap2/mmpriv.h +132 -0
  50. data/ext/minimap2/options.c +234 -0
  51. data/ext/minimap2/pe.c +177 -0
  52. data/ext/minimap2/python/README.rst +196 -0
  53. data/ext/minimap2/python/cmappy.h +152 -0
  54. data/ext/minimap2/python/cmappy.pxd +153 -0
  55. data/ext/minimap2/python/mappy.pyx +273 -0
  56. data/ext/minimap2/python/minimap2.py +39 -0
  57. data/ext/minimap2/sdust.c +213 -0
  58. data/ext/minimap2/sdust.h +25 -0
  59. data/ext/minimap2/seed.c +131 -0
  60. data/ext/minimap2/setup.py +55 -0
  61. data/ext/minimap2/sketch.c +143 -0
  62. data/ext/minimap2/splitidx.c +84 -0
  63. data/ext/minimap2/sse2neon/emmintrin.h +1689 -0
  64. data/ext/minimap2/test/MT-human.fa +278 -0
  65. data/ext/minimap2/test/MT-orang.fa +276 -0
  66. data/ext/minimap2/test/q-inv.fa +4 -0
  67. data/ext/minimap2/test/q2.fa +2 -0
  68. data/ext/minimap2/test/t-inv.fa +127 -0
  69. data/ext/minimap2/test/t2.fa +2 -0
  70. data/ext/minimap2/tex/Makefile +21 -0
  71. data/ext/minimap2/tex/bioinfo.cls +930 -0
  72. data/ext/minimap2/tex/blasr-mc.eval +17 -0
  73. data/ext/minimap2/tex/bowtie2-s3.sam.eval +28 -0
  74. data/ext/minimap2/tex/bwa-s3.sam.eval +52 -0
  75. data/ext/minimap2/tex/bwa.eval +55 -0
  76. data/ext/minimap2/tex/eval2roc.pl +33 -0
  77. data/ext/minimap2/tex/graphmap.eval +4 -0
  78. data/ext/minimap2/tex/hs38-simu.sh +10 -0
  79. data/ext/minimap2/tex/minialign.eval +49 -0
  80. data/ext/minimap2/tex/minimap2.bib +460 -0
  81. data/ext/minimap2/tex/minimap2.tex +724 -0
  82. data/ext/minimap2/tex/mm2-s3.sam.eval +62 -0
  83. data/ext/minimap2/tex/mm2-update.tex +240 -0
  84. data/ext/minimap2/tex/mm2.approx.eval +12 -0
  85. data/ext/minimap2/tex/mm2.eval +13 -0
  86. data/ext/minimap2/tex/natbib.bst +1288 -0
  87. data/ext/minimap2/tex/natbib.sty +803 -0
  88. data/ext/minimap2/tex/ngmlr.eval +38 -0
  89. data/ext/minimap2/tex/roc.gp +60 -0
  90. data/ext/minimap2/tex/snap-s3.sam.eval +62 -0
  91. data/ext/minimap2.patch +19 -0
  92. data/lib/minimap2/aligner.rb +4 -4
  93. data/lib/minimap2/alignment.rb +11 -11
  94. data/lib/minimap2/ffi/constants.rb +20 -16
  95. data/lib/minimap2/ffi/functions.rb +5 -0
  96. data/lib/minimap2/ffi.rb +4 -5
  97. data/lib/minimap2/version.rb +2 -2
  98. data/lib/minimap2.rb +51 -15
  99. metadata +97 -79
  100. data/lib/minimap2/ffi_helper.rb +0 -53
  101. data/vendor/libminimap2.so +0 -0
@@ -0,0 +1,17 @@
1
+ Q 60 32681 57 0.001744133
2
+ Q 39 3 1 0.001774569
3
+ Q 38 3 1 0.001804999
4
+ Q 35 5 1 0.001835311
5
+ Q 34 31 2 0.001894692
6
+ Q 20 11 2 0.001955154
7
+ Q 19 4 1 0.001985460
8
+ Q 15 29 5 0.002136296
9
+ Q 14 6 1 0.002166417
10
+ Q 10 11 1 0.002196193
11
+ Q 6 11 2 0.002256442
12
+ Q 5 1 1 0.002286864
13
+ Q 4 1 1 0.002317285
14
+ Q 3 36 15 0.002771602
15
+ Q 2 5 2 0.002832085
16
+ Q 1 12 9 0.003105023
17
+ Q 0 220 83 0.005594194
@@ -0,0 +1,28 @@
1
+ Q 42 16872292 669 0.000039651 16872292
2
+ Q 40 835329 636 0.000073697 17707621
3
+ Q 31 6544 2 0.000073783 17714165
4
+ Q 30 8882 6 0.000074084 17723047
5
+ Q 27 68499 9 0.000074305 17791546
6
+ Q 26 132041 81 0.000078277 17923587
7
+ Q 25 129378 96 0.000083033 18052965
8
+ Q 24 92056 382 0.000103665 18145021
9
+ Q 23 14341 402 0.000125720 18159362
10
+ Q 22 132838 146 0.000132789 18292200
11
+ Q 21 122274 124 0.000138641 18414474
12
+ Q 18 112183 103 0.000143361 18526657
13
+ Q 17 126981 213 0.000153804 18653638
14
+ Q 16 16356 208 0.000164810 18669994
15
+ Q 15 42804 782 0.000206223 18712798
16
+ Q 14 16026 318 0.000223025 18728824
17
+ Q 12 170250 814 0.000264087 18899074
18
+ Q 11 48351 1409 0.000337777 18947425
19
+ Q 8 1843 311 0.000354156 18949268
20
+ Q 7 62266 4435 0.000586276 19011534
21
+ Q 6 413997 50057 0.003150647 19425531
22
+ Q 5 404 58 0.003153568 19425935
23
+ Q 4 704 154 0.003161381 19426639
24
+ Q 3 1473 681 0.003196193 19428112
25
+ Q 2 17541 16462 0.004039875 19445653
26
+ Q 1 534344 354879 0.021693547 19979997
27
+ Q 0 11939 9917 0.022176642 19991936
28
+ U 8064
@@ -0,0 +1,52 @@
1
+ Q 60 18784147 3 0.000000160 18784147
2
+ Q 52 19002 1 0.000000213 18803149
3
+ Q 50 7152 2 0.000000319 18810301
4
+ Q 49 6797 1 0.000000372 18817098
5
+ Q 48 52188 2 0.000000477 18869286
6
+ Q 47 48775 3 0.000000634 18918061
7
+ Q 46 19447 2 0.000000739 18937508
8
+ Q 45 25983 3 0.000000896 18963491
9
+ Q 44 13455 1 0.000000949 18976946
10
+ Q 43 14573 2 0.000001053 18991519
11
+ Q 42 8697 4 0.000001263 19000216
12
+ Q 41 8645 2 0.000001368 19008861
13
+ Q 40 176603 75 0.000005264 19185464
14
+ Q 38 2503 2 0.000005368 19187967
15
+ Q 37 4117 3 0.000005523 19192084
16
+ Q 36 2924 16 0.000006356 19195008
17
+ Q 35 2323 8 0.000006772 19197331
18
+ Q 34 2344 10 0.000007292 19199675
19
+ Q 33 4279 6 0.000007603 19203954
20
+ Q 32 2092 4 0.000007810 19206046
21
+ Q 31 2625 11 0.000008382 19208671
22
+ Q 30 2828 13 0.000009057 19211499
23
+ Q 29 1581 1 0.000009108 19213080
24
+ Q 28 1543 6 0.000009420 19214623
25
+ Q 27 70916 223 0.000020948 19285539
26
+ Q 26 1288 16 0.000021777 19286827
27
+ Q 25 25551 122 0.000028065 19312378
28
+ Q 24 14345 84 0.000032390 19326723
29
+ Q 23 7308 87 0.000036878 19334031
30
+ Q 22 8358 125 0.000043325 19342389
31
+ Q 21 4836 71 0.000046983 19347225
32
+ Q 20 5888 123 0.000053325 19353113
33
+ Q 19 4656 83 0.000057600 19357769
34
+ Q 18 3948 87 0.000062081 19361717
35
+ Q 17 4418 114 0.000067954 19366135
36
+ Q 16 4226 131 0.000074702 19370361
37
+ Q 15 5760 164 0.000083144 19376121
38
+ Q 14 4697 257 0.000096384 19380818
39
+ Q 13 5246 313 0.000112503 19386064
40
+ Q 12 4170 241 0.000124908 19390234
41
+ Q 11 4095 304 0.000140557 19394329
42
+ Q 10 3857 360 0.000159087 19398186
43
+ Q 9 5300 438 0.000181617 19403486
44
+ Q 8 4206 572 0.000211050 19407692
45
+ Q 7 4676 787 0.000251541 19412368
46
+ Q 6 3923 688 0.000286924 19416291
47
+ Q 5 3294 708 0.000323333 19419585
48
+ Q 4 2936 693 0.000358965 19422521
49
+ Q 3 3928 816 0.000400897 19426449
50
+ Q 2 2613 810 0.000442533 19429062
51
+ Q 1 3515 1188 0.000503587 19432577
52
+ Q 0 567423 376636 0.019321100 20000000
@@ -0,0 +1,55 @@
1
+ Q 60 31721 27 0.000851171
2
+ Q 59 54 4 0.000975610
3
+ Q 58 29 5 0.001131933
4
+ Q 57 21 2 0.001194030
5
+ Q 56 14 4 0.001319137
6
+ Q 55 22 6 0.001506544
7
+ Q 54 12 4 0.001631475
8
+ Q 53 16 3 0.001724733
9
+ Q 51 10 1 0.001755541
10
+ Q 50 10 1 0.001786330
11
+ Q 49 11 3 0.001879699
12
+ Q 47 8 2 0.001941869
13
+ Q 46 17 1 0.001972140
14
+ Q 44 8 3 0.002065534
15
+ Q 43 10 1 0.002096174
16
+ Q 42 13 1 0.002126595
17
+ Q 41 14 3 0.002219444
18
+ Q 40 13 2 0.002281036
19
+ Q 38 17 4 0.002404747
20
+ Q 37 15 4 0.002528484
21
+ Q 36 12 1 0.002558742
22
+ Q 35 19 3 0.002650783
23
+ Q 34 12 3 0.002743313
24
+ Q 33 7 1 0.002773882
25
+ Q 32 21 3 0.002865508
26
+ Q 31 11 2 0.002926799
27
+ Q 30 14 3 0.003018891
28
+ Q 29 17 1 0.003048401
29
+ Q 28 11 2 0.003109549
30
+ Q 27 20 5 0.003262998
31
+ Q 26 11 1 0.003292948
32
+ Q 25 14 4 0.003415725
33
+ Q 24 16 5 0.003569212
34
+ Q 23 43 6 0.003750426
35
+ Q 21 15 1 0.003779664
36
+ Q 20 29 7 0.003992943
37
+ Q 19 22 2 0.004052089
38
+ Q 18 28 4 0.004172204
39
+ Q 16 25 5 0.004323390
40
+ Q 15 24 5 0.004474480
41
+ Q 14 25 5 0.004625204
42
+ Q 13 23 3 0.004714365
43
+ Q 12 22 1 0.004741963
44
+ Q 11 32 11 0.005075674
45
+ Q 10 35 7 0.005285315
46
+ Q 9 32 12 0.005648503
47
+ Q 8 33 8 0.005888126
48
+ Q 7 39 7 0.006095506
49
+ Q 6 42 14 0.006515953
50
+ Q 5 38 15 0.006966725
51
+ Q 4 37 12 0.007325113
52
+ Q 3 49 18 0.007862737
53
+ Q 2 63 21 0.008486434
54
+ Q 1 55 27 0.009292156
55
+ Q 0 153 77 0.011576593
@@ -0,0 +1,33 @@
1
+ #!/usr/bin/perl
2
+
3
+ use strict;
4
+ use warnings;
5
+ use Getopt::Std;
6
+
7
+ my %opts = (n=>33088, s=>100);
8
+ getopts('n:', \%opts);
9
+
10
+ my $pseudo = .5;
11
+ my $tot = $pseudo;
12
+ my $err = $pseudo;
13
+ my $tot_last_out = -$opts{s};
14
+ my $state = 0;
15
+ my $mapq = 0;
16
+ while (<>) {
17
+ chomp;
18
+ if (/^Q\t(\d+)\t(\d+)\t(\d+)/) {
19
+ $tot += $2;
20
+ $err += $3;
21
+ if ($tot - $tot_last_out >= $opts{s}) {
22
+ print join("\t", $1, $err/$tot, $tot / $opts{n}), "\n";
23
+ $tot_last_out = $tot;
24
+ $state = 0;
25
+ } else {
26
+ $state = 1;
27
+ $mapq = $1;
28
+ }
29
+ }
30
+ }
31
+ if ($state) {
32
+ print join("\t", $mapq, $err/$tot, $tot / $opts{n}), "\n";
33
+ }
@@ -0,0 +1,4 @@
1
+ Q 40 31897 63 0.001975107
2
+ Q 3 423 267 0.010210396
3
+ Q 2 162 120 0.013853827
4
+ Q 1 188 172 0.019038874
@@ -0,0 +1,10 @@
1
+ ./pbsim --prefix pb-1 --depth 0.1 --sample-fastq m131017_060208_42213_c100579642550000001823095604021496_s1_p0.1.subreads.fastq --length-min 1000 --length-max 30000 --seed 11 hs38.fa
2
+
3
+ bin/mason_variator -ir hs38.fa -s 1 -ov hs38-s1.vcf --snp-rate 1e-3 --small-indel-rate 2e-4 --sv-indel-rate 0 --sv-inversion-rate 0 --sv-translocation-rate 0 --sv-duplication-rate 0 --max-small-indel-size 10
4
+ bin/mason_simulator -ir hs38.fa -iv hs38-s1.vcf -n 1000000 --seed 1 -o s1_1.fq -or s1_2.fq -oa s1.sam --illumina-prob-mismatch-scale 2.5
5
+
6
+ bin/mason_variator -ir hs38.fa -s 2 -ov hs38-s2.vcf --snp-rate 1e-3 --small-indel-rate 2e-4 --sv-indel-rate 0 --sv-inversion-rate 0 --sv-translocation-rate 0 --sv-duplication-rate 0 --max-small-indel-size 10
7
+ bin/mason_simulator -ir hs38.fa -iv hs38-s2.vcf -n 1000000 --seed 2 -o mason-s2_1.fq -or mason-s2_2.fq -oa mason-s2.sam --illumina-prob-mismatch-scale 2.5 --illumina-read-length 150
8
+
9
+ bin/mason_variator -ir hs38.fa -s 3 -ov hs38-s3.vcf --snp-rate 1e-3 --small-indel-rate 2e-4 --sv-indel-rate 0 --sv-inversion-rate 0 --sv-translocation-rate 0 --sv-duplication-rate 0 --max-small-indel-size 10
10
+ bin/mason_simulator -ir hs38.fa -iv hs38-s3.vcf -n 10000000 --seed 3 -o mason-s3_1.fq -or mason-s3_2.fq -oa mason-s3.sam --illumina-prob-mismatch-scale 2.5 --illumina-read-length 150
@@ -0,0 +1,49 @@
1
+ Q 60 32070 190 0.005924540
2
+ Q 59 62 2 0.005975352
3
+ Q 58 37 5 0.006123908
4
+ Q 57 40 7 0.006333633
5
+ Q 56 39 6 0.006512032
6
+ Q 55 32 2 0.006567534
7
+ Q 54 54 2 0.006618420
8
+ Q 53 33 4 0.006735255
9
+ Q 52 39 2 0.006788866
10
+ Q 51 48 3 0.006871264
11
+ Q 50 34 2 0.006925634
12
+ Q 49 32 3 0.007011070
13
+ Q 48 35 2 0.007064967
14
+ Q 47 36 4 0.007179896
15
+ Q 46 23 1 0.007205495
16
+ Q 45 25 1 0.007230614
17
+ Q 44 17 3 0.007318716
18
+ Q 43 17 2 0.007376121
19
+ Q 42 31 5 0.007522016
20
+ Q 41 25 4 0.007638486
21
+ Q 40 26 4 0.007754541
22
+ Q 39 35 2 0.007807258
23
+ Q 37 18 4 0.007924896
24
+ Q 36 13 3 0.008013162
25
+ Q 35 15 2 0.008070411
26
+ Q 34 20 3 0.008156805
27
+ Q 33 11 1 0.008184501
28
+ Q 32 15 3 0.008272003
29
+ Q 31 25 1 0.008296107
30
+ Q 29 8 1 0.008324472
31
+ Q 28 7 2 0.008383452
32
+ Q 27 9 2 0.008441894
33
+ Q 26 30 2 0.008494888
34
+ Q 23 2 1 0.008524710
35
+ Q 22 11 3 0.008612846
36
+ Q 20 23 3 0.008697760
37
+ Q 19 6 1 0.008726479
38
+ Q 18 8 1 0.008754658
39
+ Q 16 6 1 0.008783354
40
+ Q 13 2 1 0.008813108
41
+ Q 12 4 2 0.008872604
42
+ Q 11 7 2 0.008931275
43
+ Q 10 4 3 0.009021009
44
+ Q 9 6 4 0.009140436
45
+ Q 8 6 3 0.009229559
46
+ Q 7 5 1 0.009258419
47
+ Q 6 8 3 0.009346925
48
+ Q 4 8 5 0.009495872
49
+ Q 3 17 8 0.009732801