RAZERS(1)

πŸ“¦
version 2.5.2 Β· 2 versions
, in ubuntu stonking

πŸ’– razers(1) β€” Fast Read Mapping with Sensitivity Control πŸ’–

okay bestie, sit DOWN, because razers is about to become your new hyperfixation 🧚✨ (it lives in section 1 of the manual, hence the (1), iykyk)


🌸 NAME 🌸

razers β€” Fast Read Mapping with Sensitivity Control


πŸ¦‹ SYNOPSIS πŸ¦‹

razers [OPTIONS] <GENOME FILE> <READS FILE>
razers [OPTIONS] <GENOME FILE> <MP-READS FILE1> <MP-READS FILE2>

✨ DESCRIPTION ✨

RazerS is a versatile full-sensitive read mapper based on a k-mer counting filter. It supports single and paired-end mapping, and it optimally parametrizes the filter based on a user-defined minimal sensitivity. Wanna know more? Go look at https://www.seqan.de/apps/razers3P πŸ’…

Input to RazerS is a reference genome file and either one file with single-end reads or two files containing left or right mates of paired-end reads. And listen β€” use - to read single-end reads from stdin. slay.

(c) Copyright 2009 by David Weese. 🫢


πŸ‘‘ REQUIRED ARGUMENTS πŸ‘‘

  • ARGUMENT 0 INPUT_FILE β€” A reference genome file. Valid filetypes are: .sam[.*], .raw[.*], .gbk[.*], .frn[.*], .fq[.*], .fna[.*], .ffn[.*], .fastq[.*], .fasta[.*], .fas[.*], .faa[.*], .fa[.*], .embl[.*], and .bam, where * is any of the following extensions: gz, bz2, and bgzf for transparent (de)compression. πŸŽ€
  • READS List of INPUT_FILE's β€” Either one (single-end) or two (paired-end) read files. Valid filetypes are: .sam[.*], .raw[.*], .gbk[.*], .frn[.*], .fq[.*], .fna[.*], .ffn[.*], .fastq[.*], .fasta[.*], .fas[.*], .faa[.*], .fa[.*], .embl[.*], and .bam, where * is any of the following extensions: gz, bz2, and bgzf for transparent (de)compression. 🌈

🌟 OPTIONS 🌟

bestie there are like. SO many options. every single one is here, no skips, periodt πŸ’€
  • -h, --help β€” Display the help message. πŸ†˜
  • --version β€” Display version information.

πŸŽ€ Main Options πŸŽ€

  • -f, --forward β€” Map reads only to forward strands. ➑️
  • -r, --reverse β€” Map reads only to reverse strands. ⬅️
  • -i, --percent-identity DOUBLE β€” Percent identity threshold. In range [50..100]. Default: 92.
  • -rr, --recognition-rate DOUBLE β€” Percent recognition rate. In range [80..100]. Default: 99.
  • -pd, --param-dir STRING β€” Read user-computed parameter files in the directory <DIR>.
  • -id, --indels β€” Allow indels. Default: mismatches only.
  • -ll, --library-length INTEGER β€” Paired-end library length. In range [1..inf]. Default: 220.
  • -le, --library-error INTEGER β€” Paired-end library length tolerance. In range [0..inf]. Default: 50.
  • -m, --max-hits INTEGER β€” Output only <NUM> of the best hits. In range [1..inf]. Default: 100.
  • --unique β€” Output only unique best matches β€” that's literally -m 1 -dr 0 -pa, all bundled up. πŸ’…
  • -tr, --trim-reads INTEGER β€” Trim reads to given length. Default: off. In range [14..inf].
  • -o, --output OUTPUT_FILE β€” Change output filename (use - to dump to stdout in razers format). Default: <READS FILE>.razers. Valid filetypes are: .razers, .gff, .fasta, .fa, and .eland. πŸ“
  • -v, --verbose β€” Verbose mode.
  • -vv, --vverbose β€” Very verbose mode. (yes, VERY. it's giving main character energy 🫢)

πŸ’œ Output Format Options πŸ’œ

  • -a, --alignment β€” Dump the alignment for each match (only razer or fasta format).
  • -pa, --purge-ambiguous β€” Purge reads with more than <max-hits> best matches. πŸ—‘οΈ
  • -dr, --distance-range INTEGER β€” Only consider matches with at most NUM more errors compared to the best. Default: output all.
  • -gn, --genome-naming INTEGER β€” Select how genomes are named (see Naming section below). In range [0..1]. Default: 0.
  • -rn, --read-naming INTEGER β€” Select how reads are named (see Naming section below). In range [0..2]. Default: 0.
  • -so, --sort-order INTEGER β€” Select how matches are sorted (see Sorting section below). In range [0..1]. Default: 0.
  • -pf, --position-format INTEGER β€” Select begin/end position numbering (see Coordinate section below). In range [0..1]. Default: 0.

πŸ“ Filtration Options πŸ“

  • -s, --shape STRING β€” Manually set k-mer shape. Default: 11111111111.
  • -t, --threshold INTEGER β€” Manually set minimum k-mer count threshold. In range [1..inf].
  • -oc, --overabundance-cut INTEGER β€” Set k-mer overabundance cut ratio. In range [0..1].
  • -rl, --repeat-length INTEGER β€” Skip simple-repeats of length <NUM>. In range [1..inf]. Default: 1000.
  • -tl, --taboo-length INTEGER β€” Set taboo length. In range [1..inf]. Default: 1.
  • -lm, --low-memory β€” Decrease memory usage at the expense of runtime. (trade-offs, bestie, we love a lil sacrifice πŸ™ˆ)

⭐ Verification Options ⭐

  • -mN, --match-N β€” N matches all other characters. Default: N matches nothing.
  • -ed, --error-distr STRING β€” Write error distribution to FILE.
  • -mcl, --min-clipped-len INTEGER β€” Set minimal read length for read clipping. In range [0..inf]. Default: 0.
  • -qih, --quality-in-header β€” Quality string in fasta header.

🌈 FORMATS, NAMING, SORTING, AND COORDINATE SCHEMES 🌈

RazerS supports various output formats. The output format is detected automatically from the file name suffix. no cap. πŸ”

  • .razers β€” Razer format
  • .fa, .fasta β€” Enhanced Fasta format
  • .eland β€” Eland format
  • .gff β€” GFF format

By default, reads and contigs are referred by their Fasta ids given in the input files. With the -gn and -rn options this behaviour can be changed:

  • 0 β€” Use Fasta id.
  • 1 β€” Enumerate beginning with 1.
  • 2 β€” Use the read sequence (only for short reads!). πŸ‘€

The way matches are sorted in the output file can be changed with the -so option for the following formats: razer, fasta, sam, and amos. Primary and secondary sort keys are:

  • 0 β€” 1. read number, 2. genome position
  • 1 β€” 1. genome position, 2. read number

The coordinate space used for begin and end positions can be changed with the -pf option for the razer and fasta formats:

  • 0 β€” Gap space. Gaps between characters are counted from 0.
  • 1 β€” Position space. Characters are counted from 1.

πŸ”₯ EXAMPLES πŸ”₯

  • razers example/genome.fa example/reads.fa -id -a -mN -v Map single-end reads with 4% error rate, indels, and output the alignments. Ns are considered to match everything. ✨
  • razers example/genome.fa example/reads.fa example/reads2.fa -id -mN Map paired-end reads with up to 4% errors, indels, and output concordantly mapped pairs within default library size. Ns are considered to match everything. ✨

and that's the whole thing, bestie, you ATE that documentation πŸ½οΈπŸ’… this gorgeousness is razers version 1.5.12 [tarball] β€” go forth and map, you absolute icon πŸ‘‘πŸ’–πŸ¦‹

πŸ“¦ where this page comes from