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.
- checksums.yaml +4 -4
- data/README.md +60 -76
- data/ext/Rakefile +55 -0
- data/ext/cmappy/cmappy.c +129 -0
- data/ext/cmappy/cmappy.h +44 -0
- data/ext/minimap2/FAQ.md +46 -0
- data/ext/minimap2/LICENSE.txt +24 -0
- data/ext/minimap2/MANIFEST.in +10 -0
- data/ext/minimap2/Makefile +132 -0
- data/ext/minimap2/Makefile.simde +97 -0
- data/ext/minimap2/NEWS.md +821 -0
- data/ext/minimap2/README.md +403 -0
- data/ext/minimap2/align.c +1020 -0
- data/ext/minimap2/bseq.c +169 -0
- data/ext/minimap2/bseq.h +64 -0
- data/ext/minimap2/code_of_conduct.md +30 -0
- data/ext/minimap2/cookbook.md +243 -0
- data/ext/minimap2/esterr.c +64 -0
- data/ext/minimap2/example.c +63 -0
- data/ext/minimap2/format.c +559 -0
- data/ext/minimap2/hit.c +466 -0
- data/ext/minimap2/index.c +775 -0
- data/ext/minimap2/kalloc.c +205 -0
- data/ext/minimap2/kalloc.h +76 -0
- data/ext/minimap2/kdq.h +132 -0
- data/ext/minimap2/ketopt.h +120 -0
- data/ext/minimap2/khash.h +615 -0
- data/ext/minimap2/krmq.h +474 -0
- data/ext/minimap2/kseq.h +256 -0
- data/ext/minimap2/ksort.h +153 -0
- data/ext/minimap2/ksw2.h +184 -0
- data/ext/minimap2/ksw2_dispatch.c +96 -0
- data/ext/minimap2/ksw2_extd2_sse.c +402 -0
- data/ext/minimap2/ksw2_exts2_sse.c +416 -0
- data/ext/minimap2/ksw2_extz2_sse.c +313 -0
- data/ext/minimap2/ksw2_ll_sse.c +152 -0
- data/ext/minimap2/kthread.c +159 -0
- data/ext/minimap2/kthread.h +15 -0
- data/ext/minimap2/kvec.h +105 -0
- data/ext/minimap2/lchain.c +369 -0
- data/ext/minimap2/main.c +459 -0
- data/ext/minimap2/map.c +714 -0
- data/ext/minimap2/minimap.h +410 -0
- data/ext/minimap2/minimap2.1 +725 -0
- data/ext/minimap2/misc/README.md +179 -0
- data/ext/minimap2/misc/mmphase.js +335 -0
- data/ext/minimap2/misc/paftools.js +3149 -0
- data/ext/minimap2/misc.c +162 -0
- data/ext/minimap2/mmpriv.h +132 -0
- data/ext/minimap2/options.c +234 -0
- data/ext/minimap2/pe.c +177 -0
- data/ext/minimap2/python/README.rst +196 -0
- data/ext/minimap2/python/cmappy.h +152 -0
- data/ext/minimap2/python/cmappy.pxd +153 -0
- data/ext/minimap2/python/mappy.pyx +273 -0
- data/ext/minimap2/python/minimap2.py +39 -0
- data/ext/minimap2/sdust.c +213 -0
- data/ext/minimap2/sdust.h +25 -0
- data/ext/minimap2/seed.c +131 -0
- data/ext/minimap2/setup.py +55 -0
- data/ext/minimap2/sketch.c +143 -0
- data/ext/minimap2/splitidx.c +84 -0
- data/ext/minimap2/sse2neon/emmintrin.h +1689 -0
- data/ext/minimap2/test/MT-human.fa +278 -0
- data/ext/minimap2/test/MT-orang.fa +276 -0
- data/ext/minimap2/test/q-inv.fa +4 -0
- data/ext/minimap2/test/q2.fa +2 -0
- data/ext/minimap2/test/t-inv.fa +127 -0
- data/ext/minimap2/test/t2.fa +2 -0
- data/ext/minimap2/tex/Makefile +21 -0
- data/ext/minimap2/tex/bioinfo.cls +930 -0
- data/ext/minimap2/tex/blasr-mc.eval +17 -0
- data/ext/minimap2/tex/bowtie2-s3.sam.eval +28 -0
- data/ext/minimap2/tex/bwa-s3.sam.eval +52 -0
- data/ext/minimap2/tex/bwa.eval +55 -0
- data/ext/minimap2/tex/eval2roc.pl +33 -0
- data/ext/minimap2/tex/graphmap.eval +4 -0
- data/ext/minimap2/tex/hs38-simu.sh +10 -0
- data/ext/minimap2/tex/minialign.eval +49 -0
- data/ext/minimap2/tex/minimap2.bib +460 -0
- data/ext/minimap2/tex/minimap2.tex +724 -0
- data/ext/minimap2/tex/mm2-s3.sam.eval +62 -0
- data/ext/minimap2/tex/mm2-update.tex +240 -0
- data/ext/minimap2/tex/mm2.approx.eval +12 -0
- data/ext/minimap2/tex/mm2.eval +13 -0
- data/ext/minimap2/tex/natbib.bst +1288 -0
- data/ext/minimap2/tex/natbib.sty +803 -0
- data/ext/minimap2/tex/ngmlr.eval +38 -0
- data/ext/minimap2/tex/roc.gp +60 -0
- data/ext/minimap2/tex/snap-s3.sam.eval +62 -0
- data/ext/minimap2.patch +19 -0
- data/lib/minimap2/aligner.rb +4 -4
- data/lib/minimap2/alignment.rb +11 -11
- data/lib/minimap2/ffi/constants.rb +20 -16
- data/lib/minimap2/ffi/functions.rb +5 -0
- data/lib/minimap2/ffi.rb +4 -5
- data/lib/minimap2/version.rb +2 -2
- data/lib/minimap2.rb +51 -15
- metadata +97 -79
- data/lib/minimap2/ffi_helper.rb +0 -53
- 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,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
|