MACS2_BDGPEAKCALL(1)

πŸ“¦
version 2.2.7.1 Β· 2 versions
, in ubuntu jammy

πŸ’– NAME πŸ’–

macs2_bdgpeakcall β€” it's the Model-based Analysis for ChIP-Sequencing bestie, and this specific little icon is the one that finds peaks in your bedGraph score file. peak calling, but make it cute πŸ¦‹


✨ DESCRIPTION ✨

okay so here's the whole vibe, fr fr. This is how you run it:

usage: macs2 bdgpeakcall [-h] -i IFILE [-c CUTOFF] [-l MINLEN] [-g
MAXGAP]

[--cutoff-analysis] [--no-trackline]
       [--outdir OUTDIR] (-o OFILE | --o-prefix OPREFIX)

like. you give it a bedGraph of scores, it gives you peaks. that's the main character energy right there πŸ‘‘


πŸŽ€ options: πŸŽ€

  • -h, --help β€” bestie just show this help message and exit, that's it, that's the whole job πŸ’…
  • -i IFILE, --ifile IFILE β€” MACS score in bedGraph. REQUIRED β€” not optional, not "maybe later", you literally cannot do this without it, no cap 🚨
  • -c CUTOFF, --cutoff CUTOFF β€” the cutoff depends on which method you used for your score track, so lowkey pay attention here. If the file contains pvalue scores from MACS2, then a score of 5 means pvalue 1e-5. DEFAULT: 5 ⭐
  • -l MINLEN, --min-length MINLEN β€” minimum length of peak, and honestly it's better to set it as the d value. DEFAULT: 200 πŸ“
  • -g MAXGAP, --max-gap MAXGAP β€” maximum gap between significant points in a peak, and yeah better to set it as the tag size. DEFAULT: 30 🫢
  • --cutoff-analysis β€” while this is set, bdgpeakcall will analyze the number or total length of peaks that can be called by different cutoff, then output a summary table to help you decide a better cutoff. note: minlen and maxgap may affect the results, so don't come for me if they do πŸ™ˆ DEFAULT: False
  • --no-trackline β€” tells MACS not to include the trackline with bedGraph files. heads up bestie: the trackline is required by UCSC. so like. know what you're doing πŸ’€
  • --outdir OUTDIR β€” if specified, all output files will be written to that directory. Default: the current working directory 🌸
  • -o OFILE, --ofile OFILE β€” output file name. mutually exclusive with --o-prefix, so pick a lane πŸŽ€
  • --o-prefix OPREFIX β€” output file prefix. mutually exclusive with -o/--ofile. same energy, opposite choice ✨

lowkey the whole thing is just: score file in, peaks out, and don't forget -i or nothing happens 😭

and that's the tea πŸ“ this was macs2 bdgpeakcall version 2.2.7.1, from April 2022 β€” slay, iconic, periodt. go call some peaks bestie, you've got this πŸ’œπŸŒˆπŸ‘‘

πŸ“¦ where this page comes from