RUN_SEPP.PY(1)

πŸ“¦ version 4.5.5, in ubuntu jammy, ubuntu questing

πŸ’– run_sepp.py β€” a phylogenetic placement tool πŸ’–

okay bestie, sit down, grab a boba, because we're about to talk about run_sepp.py, a User Commands girlie who lives in section 1 of your manual pages. she places stuff on trees. phylogenetically. it's giving ✨science✨ and honestly? she ate.


πŸŽ€ NAME πŸŽ€

run_sepp.py β€” a phylogenetic placement tool πŸ¦‹


🌸 DESCRIPTION 🌸

so here's the whole vibe. you run her like this:

usage: run_sepp.py [-h] [-v] [-A N] [-P N] [-F N] [-D DISTANCE] [-M
DIAMETER]
       [-S DECOMP] [-p DIR] [-rt] [-o OUTPUT] [-d OUTPUT_DIR] [-c
CONFIG] [-t TREE] [-r RAXML] [-a ALIGN] [-f FRAG] [-m MOLECULE]
       [--ignore-overlap] [-x N] [-cp CHCK_FILE] [-cpi N] [-seed N]

and what she does β€” no cap β€” is run the SEPP algorithm on an input tree, an alignment, a fragment file, and a RAxML info file. that's the whole job. she's booked and busy πŸ’…

optional arguments πŸ’«

  • -h, --help β€” shows this help message and dips. classic. πŸ™ˆ
  • -v, --version β€” shows the program's version number and exits. her little birth certificate. πŸ‘‘

🧬 DECOMPOSITION OPTIONS 🧬

these ones decide the alignment decomposition size and the taxon insertion size. if you pass None, the default is to align/place at 10% of total taxa. also super important: the alignment decomosition size must be less than the taxon insertion size. like, mandatory. it's the rules. 😀

  • -A N, --alignmentSize N β€” max alignment subset size of N. default is 10% of the total number of taxa or the placement subset size if given. ✨
  • -P N, --placementSize N β€” max placement subset size of N. default is 10% of the total number of taxa or the alignment length (whichever bigger). she said whichever bigger, periodt. πŸ’…
  • -F N, --fragmentChunkSize N β€” maximum fragment chunk size of N. this one helps you control memory, so if your machine is struggling, she's your girl. default: 20000. 🫢
  • -D DISTANCE, --distance DISTANCE β€” minimum p-distance before stopping the decomposition. default: 1. πŸ›‘
  • -M DIAMETER, --diameter DIAMETER β€” maximum tree diameter before stopping the decomposition. default: None. 🌳
  • -S DECOMP, --decomp_strategy DECOMP β€” the decomposition strategy. default: using tree branch length. iconic behaviour honestly. 🌈

πŸ“€ OUTPUT OPTIONS πŸ“€

these control the output. where the girlies go. πŸŽ€

  • -p DIR, --tempdir DIR β€” Tempfile files will be written to DIR. Full-path required. default: /tmp/sepp. don't be lazy with the path bestie. ⚠️
  • -rt, --remtemp β€” Remove tempfile directory. default: disabled. she cleans up after herself when you ask nicely. 🧹
  • -o OUTPUT, --output OUTPUT β€” output files with prefix OUTPUT. default: output. πŸ’–
  • -d OUTPUT_DIR, --outdir OUTPUT_DIR β€” output to the OUTPUT_DIR directory. full-path required. default: . πŸ—‚οΈ

πŸ“₯ INPUT OPTIONS πŸ“₯

these control the input. and listen, to actually run SEPP the following is required, so no skipping:

  • a backbone tree (newick format) 🌲
  • a RAxML_info file β€” this is the file RAxML generates while estimating the backbone tree, and Pplacer uses this info file to set model parameters 🧠
  • a backbone alignment file (fasta format) πŸ“„
  • a fasta file with the fragments 🧩

also heads up: the input sequences are assumed to be DNA unless you say otherwise. πŸ™ˆ

  • -c CONFIG, --config CONFIG β€” a config file with the options used to run SEPP. anything you pass as a command line argument overrides the config file value for that option. default: None. she respects your last word. πŸ‘‘
  • -t TREE, --tree TREE β€” input tree file (newick format). default: None. 🌳
  • -r RAXML, --raxml RAXML β€” RAxML_info file including model parameters, generated by RAxML. default: None. πŸ€–
  • -a ALIGN, --alignment ALIGN β€” aligned fasta file. default: None. πŸ“‘
  • -f FRAG, --fragment FRAG β€” fragment file. default: None. 🧩
  • -m MOLECULE, --molecule MOLECULE β€” molecule type of sequences. can be amino, dna, or rna. default: dna. choose your fighter. πŸ’…
  • --ignore-overlap β€” when a query sequence has the same name as a backbone sequence, ignore the query sequence and keep the backbone sequence. default: False. drama-free. ✌️

⭐ OTHER OPTIONS ⭐

these control how SEPP is run. the behind-the-scenes girlies. 🎬

  • -x N, --cpu N β€” use N cpus. default: number of cpus available on the machine. she'll use everything she can get. πŸ”₯
  • -cp CHCK_FILE, --checkpoint CHCK_FILE β€” checkpoint file. default: no checkpointing. 😭
  • -cpi N, --interval N β€” interval (in seconds) between checkpoint writes. has effect only with -cp provided. default: 3600. ⏰
  • -seed N, --randomseed N β€” random seed number. default: 297834. keep her consistent bestie, reproducibility is a love language. 🎲

πŸ‘€ SEE ALSO πŸ‘€

sepp.config(5) β€” go read her too, she's part of the friend group 🫢


that's the whole page bestie, go forth and place some taxa πŸ’–βœ¨

this was run_sepp.py, from the User Commands manual, dated Septbember 2020 β€” yes that's how it's spelled in the page and we love her for it, slay πŸ’…

πŸ“¦ where this page comes from