Skip to content

CLI parameters

rustar-aligner accepts STAR’s --camelCase parameter names. Defaults match STAR. This page lists the currently-supported parameters; if a parameter is missing, the binary errors out at startup rather than silently ignoring it.

Run rustar-aligner --help for the full machine-generated listing.

ParameterDefaultDescription
--runModealignReadsalignReads or genomeGenerate.
--runThreadN1Number of threads.
--runRNGseed777RNG seed for tie-breaking among equal-scoring alignments.
ParameterDefaultDescription
--genomeDir./GenomeDirPath to the genome index directory.
--genomeFastaFilesOne or more FASTA files (required for genomeGenerate).
--genomeSAindexNbases14Length of the SA pre-indexing string (log2). Lower for small genomes.
--genomeChrBinNbits18Log2 of chromosome bin size.
--genomeSAsparseD1SA sparsity (higher = less RAM, slower mapping).
ParameterDefaultDescription
--readFilesInInput FASTQ file(s); second file is mate 2 for paired-end (required for alignReads).
--readFilesCommandDecompression command, e.g. zcat for .gz.
--readMapNumber-1Number of reads to map (-1 = all).
--clip5pNbases0Bases to clip from the 5’ end of each mate.
--clip3pNbases0Bases to clip from the 3’ end of each mate.
ParameterDefaultDescription
--outFileNamePrefix./Prefix (path + filename stem) for all output files.
--outSAMtypeSAMSAM, BAM Unsorted, BAM SortedByCoordinate, or None.
--outBAMcompression1BGZF level. -1/0 = uncompressed; 18 = flate2 levels; ≥9 = max.
--limitBAMsortRAM0Max RAM (bytes) for sorted BAM. 0 = unlimited.
--outStdNoneRoute primary output to stdout: None, SAM, BAM_Unsorted, BAM_SortedByCoordinate.
ParameterDefaultDescription
--outSAMstrandFieldNoneNone or intronMotif (sets XS tag from junction motifs).
--outSAMattributesStandardTags to include: Standard, All, None, or an explicit list (e.g. NH HI AS NM nM MD).
--outSAMattrRGline-Read group line(s). Multiple blocks separated by a literal ,.
--outSAMunmappedNoneUnmapped reads in SAM: None, Within, or Within KeepPairs.
--outSAMmapqUnique255MAPQ value for uniquely-mapping reads.
--outSAMmultNmax-1Max alignments per read in SAM (-1 = all up to outFilterMultimapNmax).
ParameterDefaultDescription
--outFilterTypeNormalNormal or BySJout (re-filter by discovered SJ pass).
--outFilterMultimapNmax10Max number of multi-mapping loci. Reads exceeding this are unmapped (MultiMapTooMany).
--outFilterMultimapScoreRange1Score range for keeping multi-mappers within best score.
--outFilterMismatchNmax10Max mismatches per pair.
--outFilterMismatchNoverLmax0.3Max ratio of mismatches to mapped length.
--outFilterScoreMin0Min absolute alignment score.
--outFilterScoreMinOverLread0.66Min alignment score normalized to read length.
--outFilterMatchNmin0Min absolute matched bases.
--outFilterMatchNminOverLread0.66Min matched bases normalized to read length.
--outFilterIntronMotifsNoneNone, RemoveNoncanonical, or RemoveNoncanonicalUnannotated.
--outFilterIntronStrandsRemoveInconsistentStrandsNone or RemoveInconsistentStrands.
ParameterDefaultDescription
--outReadsUnmappedNoneNone or Fastx (write Unmapped.out.mate1/mate2).
ParameterDefaultDescription
--outSJfilterOverhangMin30 12 12 12Min overhang per motif [noncan, GT/AG, GC/AG, AT/AC].
--outSJfilterCountUniqueMin3 1 1 1Min unique-mapping reads per motif.
--outSJfilterCountTotalMin3 1 1 1Min total reads per motif.
--outSJfilterDistToOtherSJmin10 0 5 10Min distance to other SJs per motif.
--outSJfilterIntronMaxVsReadN50000 100000 200000Max intron length per supporting-read count tier.
ParameterDefaultDescription
--alignIntronMin21Min intron size (smaller gaps are deletions).
--alignIntronMax0Max intron size; 0 = auto from genome / win params.
--alignMatesGapMax0Max genomic distance between PE mates; 0 = auto.
--alignSplicedMateMapLmin0Min mapped length for spliced PE mate (absolute).
--alignSplicedMateMapLminOverLmate0.66Min mapped length for spliced PE mate (fraction).
--alignSJoverhangMin5Min overhang for novel splice junctions.
--alignSJDBoverhangMin3Min overhang for annotated junctions.
--alignSJstitchMismatchNmax0 -1 0 0Max mismatches for SJ stitching [noncan, GC/AG, AT/AC, noncan].
ParameterDefaultDescription
--scoreGap0Canonical splice junction penalty.
--scoreGapNoncan-8Non-canonical junction penalty.
--scoreGapGCAG-4GC/AG junction penalty.
--scoreGapATAC-8AT/AC junction penalty.
--scoreDelOpen-2Deletion open penalty.
--scoreDelBase-2Deletion extension penalty per base.
--scoreInsOpen-2Insertion open penalty.
--scoreInsBase-2Insertion extension penalty per base.
--scoreStitchSJshift1Max score reduction for SJ stitching shift.
--scoreGenomicLengthLog2scale-0.25Log-scaled bonus per log2(genomicLength).
ParameterDefaultDescription
--winReadCoverageRelativeMin0.5Min read coverage for an alignment window (fraction).
--winBinNbits16Log2 of window bin size for seed clustering.
--winAnchorDistNbins9Max bins for seed anchor distance.
--winFlankNbins4Bins to extend each window by on each side.
--winAnchorMultimapNmax50Max loci an anchor can map to.
--seedMultimapNmax10000Max loci a seed can map to.
--seedPerReadNmax1000Max seeds per read.
--seedPerWindowNmax50Max seeds per window.
--seedSearchStartLmax50Max distance between seed search start positions.
--seedSearchStartLmaxOverLread1.0seedSearchStartLmax normalised by read length.
--seedSearchLmax0Max seed length; 0 = unlimited.
--seedMapMin5Min mappable length for seed search termination.
--alignWindowsPerReadNmax10000Max alignment windows per read.
--alignTranscriptsPerWindowNmax100Max transcripts per window.
ParameterDefaultDescription
--sjdbGTFfileGTF file with exon annotations.
--sjdbGTFchrPrefix""Prefix to add to chromosome names from GTF (e.g. chr).
--sjdbGTFfeatureExonexonGTF feature type to use as exon.
--sjdbGTFtagExonParentTranscripttranscript_idGTF attribute for transcript ID.
--sjdbGTFtagExonParentGenegene_idGTF attribute for gene ID.
--sjdbOverhang100Overhang length around junctions in the index. Set to read_length - 1.
--sjdbScore2Extra score for alignments crossing annotated junctions.
ParameterDefaultDescription
--quantModeGeneCounts and/or TranscriptomeSAM, space-separated.
--quantTranscriptomeSAMoutputBanSingleEnd_BanIndels_ExtendSoftclipVariant for transcriptome BAM: BanSingleEnd, BanSingleEnd_ExtendSoftclip, or the default RSEM-compatible form.
ParameterDefaultDescription
--twopassModeNoneNone or Basic.
--twopass1readsN-1Reads to use in pass 1 (-1 = all).
ParameterDefaultDescription
--chimSegmentMin0Min chimeric segment length. 0 disables chimeric detection.
--chimScoreMin0Min total chimeric alignment score.
--chimScoreDropMax20Max drop in chimeric score vs read length.
--chimScoreSeparation10Min score separation for unique chimeric.
--chimMainSegmentMultNmax10Max multimapping for main chimeric segment.
--chimSegmentReadGapMax0Max read-space gap between chimeric segments.
--chimJunctionOverhangMin20Min overhang at chimeric junction.
--chimScoreJunctionNonGTAG-1Score penalty for non-GT/AG chimeric junctions.
--chimOutTypeJunctionsJunctions, WithinBAM, or both space-separated.
ParameterDefaultDescription
--readNameFilter""If set, only emit detailed alignment logs for reads with this name.