LUMPY(1)

πŸ’… LUMPY β€” ok bestie, let's find some structural variants

Heyyy so you're reading the manual page for lumpy, aka LUMPY, living in the User Commands section of the manual. Automated breakpoint detection for standard analyses β€” that's the whole vibe. Let's get into it. πŸ¦‹


✨ NAME ✨

lumpy β€” Automated breakpoint detection for standard analyses


🌸 SYNOPSIS 🌸

lumpy [OPTIONS]

That's it. That's the whole synopsis. She's a simple girl with a lot of flags. πŸ’…


🧚 DESCRIPTION 🧚

So LUMPY is a novel SV discovery framework (structural variants, bestie, not student version πŸ’€) and the way it works is that it naturally integrates multiple SV signals jointly across multiple samples. That's the whole main character move.

And the payoff? LUMPY yields improved sensitivity πŸ”₯ β€” especially when the SV signal is reduced, which happens either because you've got low coverage data or because your intra-sample variant allele frequency is low. Basically when the signal is giving nothing, LUMPY still shows up and eats. πŸ‘‘


✨ OPTIONS ✨

Every single one, no skips, no shortcuts. Let's go. 🌈

  • -g β€” Genome file (this is what defines chromosome order, so be nice to it)
  • -e β€” Show evidence for each call. Receipts! 🧾
  • -w β€” File read windows size (default 1000000)
  • -mw β€” minimum weight for a call
  • -msw β€” minimum per-sample weight for a call
  • -tt β€” trim threshold
  • -x β€” exclude file bed file
  • -t β€” temp file prefix, must be to a writeable directory ⚠️ (yes, it has to actually be writeable, do not play with her)
  • -P β€” output probability curve for each variant πŸ“ˆ
  • -b β€” output BEDPE instead of VCF
  • -sr β€” split read settings, and she takes a whole little list of them:
bam_file:<file name>,
id:<sample name>,
back_distance:<distance>,
min_mapping_threshold:<mapping quality>,
weight:<sample weight>,
min_clip:<minimum clip length>,
read_group:<string>
  • -pe β€” paired end settings, and she takes an even longer list:
bam_file:<file name>,
id:<sample name>,
histo_file:<file name>,
mean:<value>,
stdev:<value>,
read_length:<length>,
min_non_overlap:<length>,
discordant_z:<z value>,
back_distance:<distance>,
min_mapping_threshold:<mapping quality>,
weight:<sample weight>,
read_group:<string>
  • -bedpe β€” bedpe input, and here's her list:
bedpe_file:<bedpe file>,
id:<sample name>, weight:<sample weight>

πŸ’– AUTHOR πŸ’–

This manpage was written by Andreas Tille for the Debian distribution and can be used for any other usage of the program. Giving credit where credit is due, periodt. 🫢


ok that's the whole page, you're fully certified now πŸ’œ go call some variants

Signing off for real: this is lumpy 0.3.0, dated June 2020, and that's the LUMPY(1) energy. bye bestie, slay your breakpoints 🌟