chromsize

annoyed to have to create an index and cut it?
have to look for that old script every time?
got you. just get your chrom sizes. very fast.

but first, how is this better than any other option? yeah, just check the image below.

googled 'get chromosome sizes from fasta', grab every command/tool I found and benchmarked it. surprisingly, you can lose 14 seconds of your life just waiting for those chrom sizes to be calculated. crazy.

What's new on v.0.0.32?

• now accepts .2bit files as input too!
• --fasta argument now is --input (or -i) [accounts for .2bit files]

Usage

Binary

Usage: chromsize --input <FASTA/FASTA.GZ/2BIT> --output <OUTPUT> [-t <THREADS>]

Arguments:
    -i, --input <FASTA>: FASTA file
    -o, --output <OUTPUT>: path to chrom.sizes

Options:
    -t, --threads <THREADS>: number of threads [default: your max ncpus]
    -a, --accession-only     only keep the accession id part of the header (stop after blank)
    --help: print help
    --version: print version

crate: https://crates.io/crates/chromsize

Installation

to install rust and use chromsize on your system follow this steps:

1. get installer: curl https://sh.rustup.rs -sSf | sh on unix, or go here for other options

2. run cargo install chromsize (make sure ~/.cargo/bin is in your $PATH before running it)

3. use chromsize with the required arguments

Library

use chromsize;

fn main() {
    let input = PathBuf::new("/path/to/fasta.fa"); // INFO: can be .2bit too
    let output = PathBuf::new("/path/to/chrom.sizes");

    let sizes: Vec<(String, u64)> = chromsize::chromsize(&input);
    chromsize::write(sizes, &output)
}

Python

build the port to install it as a pkg:

git clone https://github.com/alejandrogzi/chromsize.git && cd chromsize/py-chromsize
hatch shell
maturin develop --release

use it as a binary wrapper [for .fa and .2bit]:

import chromsize as cs

input = "/path/to/fasta.fa" # INFO: can be .2bit too
output = "/path/to/chrom.sizes"
cs.write_chromsizes(input, output)

or just get them directly

import chromsize as cs

input = "/path/to/fasta.fa" # INFO: can be .2bit too
sizes = cs.get_chromsizes(input)

>>> print(sizes)
[
    ('chr1', 123),
    ('chr2', 456),
    ...
]

Build

to build chromsize from this repo, do:

1. get rust
2. run git clone https://github.com/alejandrogzi/chromsize.git && cd chromsize
3. run cargo run --release -- -i <FASTA/FASTA.GZ/2BIT> -o <OUTPUT>

Container image

to build the development container image:

1. run git clone https://github.com/alejandrogzi/chromsize.git && cd chromsize
2. initialize docker with start docker or systemctl start docker
3. build the image docker image build --tag chromsize .
4. run docker run --rm -v "[dir_where_your_fa_is]:/dir" chromsize -f /dir/<INPUT> -o /dir/<OUTPUT>

Conda

to use chromsize through Conda just:

1. conda install chromsize -c bioconda or conda create -n chromsize -c bioconda chromsize

Nextflow (not available yet)

Benchmark

do not believe me? run the benchmark on your own:

1. get .fa from any species you want (or download the ones I used from UCSC/NCBI)
2. install hyperfine: https://github.com/sharkdp/hyperfine
3. go to chromsize/bench and modify the ASSEMBLIES const with the .fa you've download
4. run cargo run release --bin chromsize-benchmark -- -d /dir/where/my/fastas/are -a show-output ignore-failure

here is all the info and metadata from my experiment:

which tools I used?

detailed data?

how well performs with .gz?

CHM13-T2T.fa.gz