subjunc - Online in the Cloud

PROGRAM:

NAME

subjunc - an RNA-seq aligner suitable for all purposes of RNA-seq analyses

USAGE

subjunc [options] -i <index_name> -r <input> -o <output>

Required arguments:

-i <index>
Base name of the index.

-r <string>
Name of the input file. Input formats including gzipped fastq, fastq, and fasta can
be automatically detected. If paired-end, this should give the name of file
including first reads.

Optional arguments:

-o <string>
Name of the output file. By default, the output is in BAM format.

-n <int>
Number of selected subreads, 14 by default.

-m <int>
Consensus threshold for reporting a hit (minimal number of subreads that map in
consensus) . If paired-end, this gives the consensus threshold for the anchor read.
1 by default

-M <int>
Specify the maximum number of mis-matched bases allowed in the alignment. 3 by
default. Mis-matches found in softclipped bases are not counted.

-T <int>
Number of CPU threads used, 1 by default.

-I <int>
Maximum length (in bp) of indels that can be detected. 5 by default. The program
can detect indels of up to 200bp long.

-B <int>
Maximal number of equally-best mapping locations to be reported. 1 by default. Note
that -u option takes precedence over -B.

-P <3:6>
Format of Phred scores used in input files, '3' for phred+33 and '6' for phred+64.
'3' by default.

-u Report uniquely mapped reads only. Number of mis-matched bases is used to break the
tie.

-b Convert color-space read bases to base-space read bases in the mapping output. Note
that read mapping is performed at color-space.

--SAMinput
Input reads are in SAM format.

--BAMinput
Input reads are in BAM format.

--SAMoutput
Save mapping result in SAM format.

--trim5 <int>
Trim off <int> number of bases from 5' end of each read. 0 by default.

--trim3 <int>
Trim off <int> number of bases from 3' end of each read. 0 by default.

--rg-id <string>
Add read group ID to the output.

--rg <string>
Add <tag:value> to the read group (RG) header in the output.

--DPGapOpen <int> Penalty for gap opening in short indel detection. -1 by
default.

--DPGapExt <int>
Penalty for gap extension in short indel detection. 0 by default.

--DPMismatch <int> Penalty for mismatches in short indel detection. 0 by
default.

--DPMatch <int>
Score for matched bases in short indel detection. 2 by default.

--allJunctions
Detect exon-exon junctions (both canonical and non-canonical junctions) and
structural variants in RNA-seq data. Refer to Users Guide for reporting of
junctions and fusions.

--complexIndels
Detect multiple short indels that occur concurrently in a small genomic region
(these indels could be as close as 1bp apart).

-v Output version of the program.

Optional arguments for paired-end reads:

-R <string>
Name of the file including second reads.

-p <int>
Consensus threshold for the non-anchor read (receiving less votes than the anchor
read from the same pair). 1 by default.

-d <int>
Minimum fragment/insert length, 50bp by default.

-D <int>
Maximum fragment/insert length, 600bp by default.

-S <ff:fr:rf>
Orientation of first and second reads, 'fr' by default ( forward/reverse).

Refer to Users Manual for detailed description to the arguments.

Use subjunc online using onworks.net services