diff --git a/main.c b/main.c index beee9ff..46a4711 100644 --- a/main.c +++ b/main.c @@ -80,8 +80,7 @@ static ko_longopt_t long_options[] = { { "ds", ko_no_argument, 355 }, { "rmq-inner", ko_required_argument, 356 }, { "spsc", ko_required_argument, 357 }, - { "spsc0", ko_required_argument, 358 }, - { "junc-pen", ko_required_argument, 358 }, // an alias to --spsc0; this is not a typo + { "junc-pen", ko_required_argument, 358 }, { "dbg-seed-occ", ko_no_argument, 501 }, { "help", ko_no_argument, 'h' }, { "max-intron-len", ko_required_argument, 'G' }, diff --git a/minimap.h b/minimap.h index 6311379..9678d1f 100644 --- a/minimap.h +++ b/minimap.h @@ -5,7 +5,7 @@ #include #include -#define MM_VERSION "2.28-r1215-dirty" +#define MM_VERSION "2.28-r1216-dirty" #define MM_F_NO_DIAG (0x001LL) // no exact diagonal hit #define MM_F_NO_DUAL (0x002LL) // skip pairs where query name is lexicographically larger than target name diff --git a/minimap2.1 b/minimap2.1 index 9e15ccd..2a2a985 100644 --- a/minimap2.1 +++ b/minimap2.1 @@ -433,6 +433,21 @@ on SIRV data, please add .B --splice-flank=no to the command line. .TP +.BR --spsc \ FILE +Splice scores []. Each line consists of five fields: 1) contig, 2) offset, 3) `+' or `-', 4) `D' or `A', and 5) score, +where offset is the number of bases before a splice junction, `D' indicates the +line corresponds to a donor site and `A' for an acceptor site. +A positive score suggests the junction is preferred and a negative score +suggests the junction is not preferred. +.TP +.BR --junc-pen \ INT +Penalty for a position not in FILE specified by +.B --spsc +[5]. Effective with +.B --spsc +but not +.BR --junc-bed . +.TP .BR --junc-bed \ FILE Gene annotations in the BED12 format (aka 12-column BED), or intron positions in 5-column BED. With this option, minimap2 prefers splicing in annotations. @@ -440,9 +455,10 @@ BED12 file can be converted from GTF/GFF3 with `paftools.js gff2bed anno.gtf' []. .TP .BR --junc-bonus \ INT -Score bonus for a splice donor or acceptor found in annotation (effective with -.BR --junc-bed ) -[9]. +Score bonus for a splice donor or acceptor found in annotation [9]. Effective with +.B --junc-bed +but not +.BR --spsc . .TP .BI --end-seed-pen \ INT Drop a terminal anchor if