document --spsc

This commit is contained in:
Heng Li
2024-10-12 23:27:46 -04:00
parent 68179ed195
commit 3e1ab2951d
3 changed files with 21 additions and 6 deletions

3
main.c
View File

@@ -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' },

View File

@@ -5,7 +5,7 @@
#include <stdio.h>
#include <sys/types.h>
#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

View File

@@ -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