CNVKIT(1)

๐Ÿ“ฆ
version 0.9.9 ยท 2 versions
, in ubuntu jammy

๐Ÿ‘‘ CNVKIT(1) โ€” the girlies' guide to the copy number toolkit โœจ

okay so. bestie. this is cnvkit, and it lives in the User Commands section of the manual. she's a command-line toolkit for copy number analysis. that's the whole vibe. let's get into it ๐Ÿ’…


๐ŸŒธ NAME ๐ŸŒธ

cnvkit โ€” a command-line toolkit for copy number analysis.


๐Ÿ’– DESCRIPTION ๐Ÿ’–

okay so the usage is literally this, no cap, here's the exact shape of her:

usage: cnvkit [-h]

       {batch,target,access,antitarget,autobin,coverage,reference,fix,segment,call,diagram,scatter,heatmap,breaks,genemetrics,gainloss,sex,gender,metrics,segmetrics,bintest,import-picard,import-seg,import-theta,import-rna,export,version}
       ...

CNVkit, a command-line toolkit for copy number analysis. she contains multitudes ๐Ÿฆ‹


โœจ positional arguments โœจ

like, here's the whole squad. every single one of them, in order, periodt:

  • {batch,target,access,antitarget,autobin,coverage,reference,fix,segment,call,diagram,scatter,heatmap,breaks,genemetrics,gainloss,sex,gender,metrics,segmetrics,bintest,import-picard,import-seg,import-theta,import-rna,export,version} โ€” these are the Sub-commands. throw -h at any of them for more info, iykyk ๐Ÿ”

and here's the whole cast list with their little blurbs, because nobody gets left out:

  • batch โ€” Run the complete CNVkit pipeline on one or more BAM files. the whole glow-up, start to finish ๐ŸŒˆ
  • target โ€” Transform bait intervals into targets more suitable for CNVkit. glow-up for your intervals ๐Ÿ’ซ
  • access โ€” List the locations of accessible sequence regions in a FASTA file. where can we actually go ๐Ÿ—บ๏ธ
  • antitarget โ€” Derive off-target ("antitarget") bins from target regions. the off-screen characters, but make them bins ๐ŸŽญ
  • autobin โ€” Quickly calculate reasonable bin sizes from BAM read counts. she does the math so you don't have to ๐Ÿงฎ
  • coverage โ€” Calculate coverage in the given regions from BAM read depths. how deep does it go ๐Ÿ’ง
  • reference โ€” Compile a coverage reference from the given files (normal samples). building the baseline, the control group, the normal friend ๐Ÿงš
  • fix โ€” Combine target and antitarget coverages and correct for biases. Adjust raw coverage data according to the given reference, correct potential biases and recenter. she's the therapist of the pipeline, correcting all the biases ๐Ÿ˜Œ
  • segment โ€” Infer copy number segments from the given coverage table. connecting the dots, literally ๐Ÿ”—
  • call โ€” Call copy number variants from segmented log2 ratios. saying it out loud ๐Ÿ“ข
  • diagram โ€” Draw copy number (log2 coverages, segments) on chromosomes as a diagram. If both the raw probes and segments are given, show them side-by-side on each chromosome (segments on the left side, probes on the right side). art ๐ŸŽจ left side segments, right side probes, no notes.
  • scatter โ€” Plot probe log2 coverages and segmentation calls together. the two besties, side by side ๐Ÿซถ
  • heatmap โ€” Plot copy number for multiple samples as a heatmap. everyone on one vibe board ๐Ÿ”ฅ
  • breaks โ€” List the targeted genes in which a copy number breakpoint occurs. where it all went down ๐Ÿ’€
  • genemetrics โ€” Identify targeted genes with copy number gain or loss. who's up, who's down ๐Ÿ“ˆ๐Ÿ“‰
  • sex โ€” Guess samples' sex from the relative coverage of chromosomes X and Y. ๐Ÿงฌ
  • metrics โ€” Compute coverage deviations and other metrics for self-evaluation. the self-check, the mirror moment ๐Ÿชž
  • segmetrics โ€” Compute segment-level metrics from bin-level log2 ratios. zooming out, still counting ๐Ÿ”ญ
  • bintest โ€” Test for single-bin copy number alterations. one bin, one test, one diva ๐Ÿ’…
  • import-picard โ€” Convert Picard CalculateHsMetrics tabular output to CNVkit .cnn files. The input file is generated by the PER_TARGET_COVERAGE option in the CalculateHsMetrics script in Picard tools. If 'antitarget' is in the input filename, the generated output filename will have the suffix '.antitargetcoverage.cnn', otherwise '.targetcoverage.cnn'. filename detective work, she reads the name and picks the suffix ๐Ÿ”Ž
  • import-seg โ€” Convert a SEG file to CNVkit .cns files. translation era ๐ŸŒ
  • import-theta โ€” Convert THetA output to a BED-like, CNVkit-like tabular format. Equivalently, use the THetA results file to convert CNVkit .cns segments to integer copy number calls. two-for-one, she ate ๐Ÿฝ๏ธ
  • import-rna โ€” Convert a cohort of per-gene log2 ratios to CNVkit .cnr format. rna girlies, this one's for you ๐ŸŽ€
  • export โ€” Convert CNVkit output files to another format. leaving the group chat, going mainstream ๐Ÿ“ค
  • version โ€” Display this program's version. her birth chart ๐Ÿช
heads up bestie: gainloss and gender are both chilling in the sub-command roster above โ˜๏ธ but this page gives them no blurb of their own. so don't come for me, I'm not making anything up ๐Ÿ’€

๐ŸŒŸ optional arguments ๐ŸŒŸ

  • -h, --help โ€” show this help message and exit. the emergency bestie, always there ๐Ÿซ‚

See the online manual for details: https://cnvkit.readthedocs.io โ€” like, actually go read it, it's giving documentation โœจ


and that's the page, bestie!! this version is cnvkit 0.9.9, dated October 2021, and the whole thing is filed under CNVKIT(1). stay slaying, run your pipelines, love u bye ๐Ÿ’–๐Ÿฆ‹๐ŸŒธ

๐Ÿ“ฆ where this page comes from