Compare commits

..
4 Commits
Author SHA1 Message Date
Heng Li 6840370f3c Release minimap2-2.0 (r275) 2017-08-08 21:16:25 -04:00
Heng Li 8c2917b391 Added citation. 2017-08-06 21:30:47 -04:00
Heng Li 46de0fbdad Merge pull request #8 from rikuu/master
fix self-comparison in index parameter override check
2017-07-30 15:10:29 -04:00
Riku Walve 9e09c1ae72 fix self-comparison in index parameter override check 2017-07-30 21:46:25 +03:00
4 changed files with 47 additions and 16 deletions
+13
View File
@@ -1,3 +1,16 @@
Release 2.0-r275 (8 August 2017)
--------------------------------
This release is identical to version 2.0rc1, except the version number. It is
described and evaluated in the following technical report:
* Li, H. (2017). Minimap2: fast pairwise alignment for long DNA sequences.
[arXiv:1708.01492v1](https://arxiv.org/abs/1708.01492v1).
(2.0: 8 August 2017, r275)
Release 2.0rc1-r232 (30 July 2017) Release 2.0rc1-r232 (30 July 2017)
---------------------------------- ----------------------------------
+4
View File
@@ -31,6 +31,10 @@ and [longISLND][longislnd]), better chaining and the ability to produce CIGAR
with fast extension alignment (see also [libgaba][gaba] and [ksw2][ksw2]) and with fast extension alignment (see also [libgaba][gaba] and [ksw2][ksw2]) and
piece-wise affine gap cost. piece-wise affine gap cost.
If you use minimap2 in your work, please consider to cite:
> Li, H. (2017). Minimap2: fast pairwise alignment for long DNA sequences. [arXiv:1708.01492](https://arxiv.org/abs/1708.01492).
## Installation ## Installation
For modern x86-64 CPUs, just type `make` in the source code directory. This For modern x86-64 CPUs, just type `make` in the source code directory. This
+2 -2
View File
@@ -8,7 +8,7 @@
#include "minimap.h" #include "minimap.h"
#include "mmpriv.h" #include "mmpriv.h"
#define MM_VERSION "2.0rc1-r232" #define MM_VERSION "2.0-r275"
void liftrlimit() void liftrlimit()
{ {
@@ -194,7 +194,7 @@ int main(int argc, char *argv[])
mm_idx_t *mi = 0; mm_idx_t *mi = 0;
if (fpr) { if (fpr) {
mi = mm_idx_load(fpr); mi = mm_idx_load(fpr);
if (idx_par_set && mm_verbose >= 2 && (mi->k != k || mi->w != w || mi->is_hpc != mi->is_hpc)) if (idx_par_set && mm_verbose >= 2 && (mi->k != k || mi->w != w || mi->is_hpc != is_hpc))
fprintf(stderr, "[W::%s::%.3f*%.2f] Indexing parameters on the command line (-k/-w/-H) overridden by parameters in the prebuilt index.\n", fprintf(stderr, "[W::%s::%.3f*%.2f] Indexing parameters on the command line (-k/-w/-H) overridden by parameters in the prebuilt index.\n",
__func__, realtime() - mm_realtime0, cputime() / (realtime() - mm_realtime0)); __func__, realtime() - mm_realtime0, cputime() / (realtime() - mm_realtime0));
} else if (!mm_bseq_eof(fp)) { } else if (!mm_bseq_eof(fp)) {
+28 -14
View File
@@ -1,4 +1,4 @@
.TH minimap2 1 "30 July 2017" "minimap2-2.0rc1-r232" "Bioinformatics tools" .TH minimap2 1 "8 August 2017" "minimap2-2.0-r275" "Bioinformatics tools"
.SH NAME .SH NAME
.PP .PP
minimap2 - mapping and alignment between collections of DNA sequences minimap2 - mapping and alignment between collections of DNA sequences
@@ -247,35 +247,49 @@ are:
.RS .RS
.TP 8 .TP 8
.B map-pb .B map-pb
PacBio/Oxford Nanopore read to reference mapping (-Hk19) PacBio/Oxford Nanopore read to reference mapping
.RB ( -Hk19 )
.TP .TP
.B map10k .B map10k
The same as The same as
.B map-pb .B map-pb
(-Hk19) .RB ( -Hk19 )
.TP .TP
.B map-ont .B map-ont
Slightly more sensitive for Oxford Nanopore to reference mapping (-k15). For Slightly more sensitive for Oxford Nanopore to reference mapping
PacBio reads, HPC minimizers consistently leads to faster performance and more .RB ( -k15 ).
sensitive results in comparison to normal minimizers. For Oxford Nanopore data, For PacBio reads, HPC minimizers consistently leads to faster performance and
normal minimizers are better, though not much. The effectiveness of HPC is more sensitive results in comparison to normal minimizers. For Oxford Nanopore
determined by the sequencing error mode. data, normal minimizers are better, though not much. The effectiveness of HPC
is determined by the sequencing error mode.
.TP .TP
.B asm5 .B asm5
Long assembly to reference mapping (-k19 -w19 -A1 -B19 -O39,81 -E3,1 -s200 -z200). Long assembly to reference mapping
.RB ( -k19
.B -w19 -A1 -B19 -O39,81 -E3,1 -s200
.BR -z200 ).
Typically, the alignment will not extend to regions with 5% or higher sequence Typically, the alignment will not extend to regions with 5% or higher sequence
divergence. Only use this preset if the average divergence is far below 5%. divergence. Only use this preset if the average divergence is far below 5%.
.TP .TP
.B asm10 .B asm10
Long assembly to reference mapping (-k19 -w19 -A1 -B9 -O16,41 -E2,1 -s200 -z200). Up Long assembly to reference mapping
to 10% sequence divergence. .RB ( -k19
.B -w19 -A1 -B9 -O16,41 -E2,1 -s200
.BR -z200 ).
Up to 10% sequence divergence.
.TP 8 .TP 8
.B ava-pb .B ava-pb
PacBio all-vs-all overlap mapping (-Hk19 -w5 -Xp0 -m100 -K500m -g10000 --max-chain-skip 25) PacBio all-vs-all overlap mapping
.RB ( -Hk19
.B -w5 -Xp0 -m100 -K500m -g10000 --max-chain-skip
.BR 25 ).
.TP 8 .TP 8
.B ava-ont .B ava-ont
Oxford Nanopore all-vs-all overlap mapping (-k15 -w5 -Xp0 -m100 -K500m -g10000 Oxford Nanopore all-vs-all overlap mapping
--max-chain-skip 25). Similarly, the major difference from .RB ( -k15
.B -w5 -Xp0 -m100 -K500m -g10000 --max-chain-skip
.BR 25 ).
Similarly, the major difference from
.B ava-pb .B ava-pb
is that this preset is not using HPC minimizers. is that this preset is not using HPC minimizers.
.RE .RE