EXPRESS(1)

πŸ’– EXPRESS(1) πŸ’–

okay bestie, welcome to EXPRESS(1) β€” that's the manual page name, don't @ me β€” aka berkeley-express, the streaming quantification icon for high-throughput sequencing. she's a User Command, she's lowkey legendary, and she's about to make your RNA-seq life so much easier. let's get into it 🧚✨


🌸 NAME 🌸

berkeley-express β€” Streaming quantification for high-throughput sequencing πŸ’…


πŸ¦‹ DESCRIPTION πŸ¦‹

so the way that eXpress works is, you can either hand her files directly or pipe stuff straight into her. both are valid, both are iconic. pick your fighter:

File Usage:

express [options] <target_seqs.fa> <hits.(sam/bam)>

Piped Usage:

bowtie [options] -S <index> <reads.fq> | express [options] <target_seqs.fa>

fr fr, she reads from stdin if you pipe it, no cap 🫢


πŸŽ€ Required arguments:

like, these ones are not optional, bestie. she needs them. periodt.

  • <target_seqs.fa> β€” target sequence file in fasta format 🌱
  • <hits.(sam/bam)> β€” read alignment file in SAM or BAM format πŸ“Ž

⭐ Standard Options:

  • -h, --help β€” produce help message. she's helpful like that πŸ’β€β™€οΈ
  • -o, --output-dir arg (=.) β€” write all output files to this directory πŸ“ default is literally just .
  • -D, --preprocess β€” run preprocess script for eXpressD ✨
  • -m, --frag-len-mean arg (=200) β€” prior estimate for average fragment length πŸ“
  • -s, --frag-len-stddev arg (=80) β€” prior estimate for fragment length std deviation πŸ“
  • -H, --haplotype-file arg β€” path to a file containing haplotype pairs πŸ‘―
  • -B, --additional-batch arg (=0) β€” number of additional batch EM rounds after initial online round πŸ”
  • -O, --additional-online arg (=0) β€” number of additional online EM rounds after initial online round πŸ”‚
  • -L, --max-read-len arg (=250) β€” maximum allowed length of a read πŸ“–
  • --output-align-prob β€” output alignments (sam/bam) with probabilistic assignments 🎲
  • --output-align-samp β€” output alignments (sam/bam) with sampled assignments 🎯
  • --fr-stranded β€” accept only forward->reverse alignments (second-stranded protocols) βž‘οΈβ¬…οΈ
  • --rf-stranded β€” accept only reverse->forward alignments (first-stranded protocols) β¬…οΈβž‘οΈ
  • --f-stranded β€” accept only forward single-end alignments (second-stranded protocols) ➑️
  • --r-stranded β€” accept only reverse single-end alignments (first-stranded protocols) ⬅️
  • --no-update-check β€” disables automatic check for update via web πŸ™ˆ
  • --logtostderr β€” prints all logging messages to stderr πŸ“’

πŸ‘‘ Advanced Options:

this is where the main character energy lives. buckle up πŸ“

  • -f, --forget-param arg (=0.84999999999999998) β€” sets the 'forgetting factor' parameter (0.5 < c <= 1) πŸ’­
  • --library-size arg β€” specifies library size for FPKM instead of calculating from alignments πŸ“Š
  • --max-indel-size arg (=10) β€” sets the maximum allowed indel size, affecting geometric indel prior 🧬
  • --calc-covar β€” calculate and output covariance matrix πŸ”’
  • --expr-alpha arg (=0.0050000000000000001) β€” sets the strength of the prior, per bp βš–οΈ
  • --stop-at arg (=0) β€” sets the number of fragments to process, disabled with 0 πŸ›‘
  • --burn-out arg (=5000000) β€” sets number of fragments after which to stop updating auxiliary parameters πŸ”₯
  • --no-bias-correct β€” disables bias correction 🚫
  • --no-error-model β€” disables error modelling 🚫
  • --aux-param-file arg β€” path to file containing auxiliary parameters to use instead of learning πŸ“‚

πŸ’– reminder bestie: every option here is exactly as written β€” no renaming, no merging, no inventing. copy them character for character or she will Not be happy with you. πŸ™ˆ

and that's the whole vibe check, bestie 🌈 go quantify something gorgeous. this page is for express version 1.5.1, dated February 2015 β€” yeah she's vintage but she ATE, and honestly? still slay. πŸ’…βœ¨πŸ¦‹πŸ‘‘