Crate chromsize

Crate chromsize 

Source
Expand description

chromsize Alejandro Gonzales-Irribarren, 2024

chromsize is a utility designed to extract chromosome names and their corresponding lengths from input and 2bit files. It supports both plain and gzipped input formats [and .2bit] and offers an option to include only the accession ID from the input headers.

§Usage

To use chromsize, you typically provide the input input or 2bit file and specify the desired output file.

chromsize [OPTIONS] --input <input> --output <OUTPUT>

§Options

Here’s a breakdown of the available command-line options:

  • -i, --input <input>

    • Purpose: Specifies the path to the input input file. This is a required option.
    • Example: --input genome.input or --input sequences.input.gz

  • -o, --output <OUTPUT>

    • Purpose: Specifies the path where the output chromosome sizes will be written. This is a required option. The output will typically be a tab-separated file with chromosome names and their lengths.
    • Example: --output chrom_lengths.txt

  • -t, --threads <THREADS>

    • Purpose: Sets the number of threads to use for processing. This can speed up processing for large input files.
    • Default: 8
    • Example: --threads 4 (to use 4 threads)

  • -a, --accession-only

    • Purpose: A flag that, when present, instructs chromsize to only keep the accession ID part of the input header. This means it will stop reading the header at the first blank space. If omitted, the entire header line up to the first newline character will be used as the chromosome name.
    • Example: --accession-only

  • -h, --help

    • Purpose: Displays a help message with usage information and available options.
    • Example: chromsize --help

  • -V, --version

    • Purpose: Prints the version information of the chromsize tool.
    • Example: chromsize --version

§Example Usage Scenarios

  1. Get chromosome sizes from a plain input file, using full headers, with default threads:

    chromsize --input input.fa --output chrom_sizes.txt

  2. Get chromosome sizes from a gzipped input file, extracting only accession IDs, using 4 threads:

    chromsize --input input.input.gz --output accession_sizes.txt --accession-only --threads 4

Re-exports§

pub use size::*;

Modules§

size
chromsize Alejandro Gonzales-Irribarren, 2024