RAPTOR-PREPARE(1)
π RAPTOR-PREPARE β ok bestie, we're pre-filtering nucleotides today π
This is the manual page for raptor-prepare, section 1. Pull up a chair, grab a strawberry π, and let's get into it.
β¨ NAME β¨
Raptor-prepare β a fast and space-efficient pre-filter for querying very large collections of nucleotide sequences. She's giving speed, she's giving efficiency, she's giving "your RAM could never" π
π¦ SYNOPSIS π¦
raptor prepare --input <file> --output <directory> [--threads <number>]
[--quiet] [--kmer <number>|--shape <01-pattern>] [--window <number>]
[--kmer-count-cutoff <number>|--use-filesize-dependent-cutoff]
Those brackets? Optional. That vertical bar? Pick one, not both β no cap. π₯
πΈ DESCRIPTION πΈ
So here's the tea: raptor prepare computes minimisers for use with raptor layout and raptor build. That's the whole job. She's the prep girlie so the other raptor commands can eat later. π
And the best part? If you crash π β or you just wanna run it in pieces β she can continue where she left off. In multiple runs, even. Resilient queen energy, honestly.
π OPTIONS π
π General options π
--input(std::filesystem::path) β The file that contains your file names. Real talk: this file must hold at least one file path per line, and if a line has multiple paths they're separated by a whitespace. Each line = one bin. Now, valid extensions for those paths:[minimiser]if you're using preprocessed input fromraptor prepare, or[embl,fasta,fa,fna,ffn,faa,frn,fas,fastq,fq,genbank,gb,gbk,sam], and that list may be followed by[bz2,gz,bgzf]. The input file has to exist and you must have read permissions. No permissions, no minimisers, periodt. π«--output(std::filesystem::path) β A valid path for the output directory. β¨
She will create a minimiser.list inside the output directory. What's in it? A list of the generated minimiser files, in the same order as the input. Keeping it tidy, keeping it ordered, we love her.
β οΈ MAJOR heads up bestie: when you manually delete a.in_progressfile, you also have to delete the matching.headerand.minimiserfile! Do not leave those orphans hanging, it's giving chaos π
Here's every output file she creates for each file:
*.header: Contains the shape, window size, cutoff and minimiser count.*.minimiser: Contains binary minimiser values, one minimiser per line.*.in_progress: Temporary file to track process. Deleted after finishing computation.
--threads(unsigned 8 bit integer) β How many threads you wanna throw at it. Default: 1. Value must be a positive integer. Go wild, but not too wild πββοΈ--quietβ Do not print time and memory usage. For the girlies who don't wanna see the stats, just the results π€«
π k-mer options π
--kmer(unsigned 8 bit integer) β The k-mer size. Default: 20. Value must be in range[1,32]. Stay in your lane, range-wise. π§--window(unsigned 32 bit integer) β The window size. Default: k-mer size. Value must be a positive integer.--shape(std::string) β The shape to use for k-mers. Mutually exclusive with--kmerβ pick a struggle, you can't have both. Parsed from right to left, which is lowkey iconic. Default:11111111111111111111(that's a k-mer of size 20). Value must match the pattern'[01]+'. Yes, literally just zeroes and ones. Binary baddie. π
π«Ά Processing options π«Ά
--kmer-count-cutoff(unsigned 8 bit integer) β Only store k-mers with at least (>=) x occurrences. Mutually exclusive with--use-filesize-dependent-cutoff. Default: 1. Value must be in range[1,254].--use-filesize-dependent-cutoffβ Apply cutoffs from Mantis(Pandey et al., 2018). Mutually exclusive with--kmer-count-cutoff. Again β one or the other, bestie π
β Common options β
-h,--helpβ Prints the help page.-hh,--advanced-helpβ Prints the help page including advanced options. The deluxe edition β¨--versionβ Prints the version information.--copyrightβ Prints the copyright/license information.--export-help(std::string) β Export the help page information. Value must be one of[html, man, ctd, cwl].
π₯ EXAMPLES π₯
k-mer and window, the classic combo:
raptor prepare --input bins.list --output some_directory --kmer 20 --window 24
Wanna cut off the low-count k-mers?
raptor prepare --input bins.list --output some_directory --kmer-count-cutoff 2
Or let the file size decide for you, main character style:
raptor prepare --input bins.list --output some_directory --use-filesize-dependent-cutoff
π VERSION π
- Last update: Unavailable
- Raptor-prepare version: 3.0.1 (commit unavailable)
- Sharg version: 1.1.1
- SeqAn version: 3.4.0
π URL π
https://github.com/seqan/raptor β go follow, go star, go say hi π«Ά
π LEGAL π
- Raptor-prepare Copyright: BSD 3-Clause License
- Author: Enrico Seiler
- Contact:
enrico.seiler@fu-berlin.de - SeqAn Copyright: 2006-2023 Knut Reinert, FU-Berlin; released under the 3-clause BSDL.
And listen β if you're using this in your academic works, please cite:
Raptor: A fast and space-efficient pre-filter for querying very large collections of nucleotide sequences; Enrico Seiler, Svenja Mehringer, Mitra Darvish, Etienne Turc, and Knut Reinert; iScience 2021 24 (7): 102782. doi: https://doi.org/10.1016/j.isci.2021.102782
For full copyright and/or warranty information see --copyright. π
And that's the whole page, bestie π this is raptor-prepare version 3.0.1 (commit u... Unavailable, from the RAPTOR-PREPARE(1) manual β last update info says Unavailable, so don't come for me about dates. Now go forth and minimise πβ¨
π¦ where this page comes from
- ubuntu resolute: seqan-raptor 3.0.1+ds-11 Β· BSD-3-Clause
- ubuntu stonking: seqan-raptor 3.0.1+ds-12build1 Β· BSD-3-Clause
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash on : a language model's rewrite of the original, under its terms