CNVKIT_DIAGRAM(1)
β¨ NAME β¨
cnvkit_diagram β draw copy number (log2 coverages, segments) on chromosomes as a diagram π
okay so real quick, bestie: this is the cnvkit_diagram(1) page, and here's the tea π«
If both the raw probes and the segments are given, it shows them side-by-side on each chromosome β segments on the left side, probes on the right side. Iconic layout, no cap.
π DESCRIPTION π
Here's the whole vibe, in usage form:
usage: cnvkit.py diagram [-h] [-s SEGMENT] [-c CHROMOSOME] [-t THRESHOLD]
[-m MIN_PROBES] [-y]
[-x {m,y,male,Male,f,x,female,Female}] [--no-shift-xy] [-o FILENAME] [--title TITLE] [--no-gene-labels] [filename]
That's the entire shape of the thing. Let's break it down, babe πΈ
π positional arguments π
filenameβ your processed coverage data file (*.cnr), aka the output of thefixsub-command. That's the glow-up file right there π
π₯ options π₯
-h,--helpβ shows this help message and then dips. bye π-s SEGMENT,--segment SEGMENTβ the segmentation calls (.cns), which is the output of thesegmentcommand. feed it the segments and it will slay-c CHROMOSOME,--chromosome CHROMOSOMEβ the chromosome to display, e.g.'chr1'. and no, you can't do a chromosomal range here β no range allowed, periodt π«-t THRESHOLD,--threshold THRESHOLDβ the copy number change threshold used to label genes. [Default:0.5]-m MIN_PROBES,--min-probes MIN_PROBESβ the minimum number of covered probes needed before a gene gets labeled. [Default:3]-y,--male-reference,--haploid-x-referenceβ assume the inputs were normalized to a male reference. so basically: female samples will have+1log-CNR of chrX; otherwise male samples would have-1chrX. pick your fighter βοΈ-x {m,y,male,Male,f,x,female,Female},--sample-sex {m,y,male,Male,f,x,female,Female},-g {m,y,male,Male,f,x,female,Female},--gender {m,y,male,Male,f,x,female,Female}β tell it the sample's chromosomal sex, male or female. (Otherwise it just guesses from the X and Y coverage. lowkey impressive of it, honestly π)--no-shift-xyβ don't adjust the X and Y chromosomes according to sample sex. hands off the X and Y, thank you π-o FILENAME,--output FILENAMEβ the output PDF file name. that's your final look, bestie π
πΈ Plot aesthetics πΈ
this section is literally just about making it pretty, which is of course the most important part π
--title TITLEβ the plot title. [Default: sample ID, fromfilenameor-i]--no-gene-labelsβ disables thegene_namelabels on the plot. super useful when a lot of CNV were called and the labels are all fighting for the aux cord π
ngl bestie, if your plot is a chaotic mess of overlapping gene names, --no-gene-labels is your best friend. she's giving clean.
that's the whole page, you ate it, go draw some chromosomes π¦
signing off with love from cnvkit.py diagram 0.9.10, dated July 2023 πβ¨
π¦ where this page comes from
- debian trixie: cnvkit 0.9.12-1 Β· Apache-2.0
- ubuntu noble: cnvkit 0.9.10-1 Β· Apache-2.0
- ubuntu questing: cnvkit 0.9.12-1
- ubuntu resolute: cnvkit 0.9.12-1build1 Β· Apache-2.0
- ubuntu stonking: cnvkit 0.9.13-2 Β· Apache-2.0
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash on : a language model's rewrite of the original, under its terms