GT-SHREDDER(1)
๐ GT-SHREDDER(1) ๐
okay bestie, welcome to the GenomeTools Manual โจ pull up a chair, we're doing gt-shredder today ๐
โจ NAME โจ
gt-shredder โ shredder sequence file(s) into consecutive pieces of random length ๐ฆ
๐ธ SYNOPSIS ๐ธ
gt shredder [option ...] [sequence_file ...]
๐ DESCRIPTION ๐
okay so here's the whole vibe. every flag, every default, no cap:
-coverage[value] โ set the number of times the sequence_file is shreddered (default:1)-minlength[value] โ set the minimum length of the shreddered fragments (default:300)-maxlength[value] โ set the maximum length of the shreddered fragments (default:700)-overlap[value] โ set the overlap between consecutive pieces (default:0)-sample[value] โ take samples of the generated sequences pieces with the given probability (default:1.000000)-clipdesc[yes|no] โ clip descriptions after first space (fooled by\t,\netc) adds offset and length to ensure unique identifier (default:no)-width[value] โ set output width for FASTA sequence printing (0disables formatting) (default:0)-o[filename] โ redirect output to specified file (default:undefined)-gzip[yes|no] โ write gzip compressed output file (default:no)-bzip2[yes|no] โ write bzip2 compressed output file (default:no)-force[yes|no] โ force writing to output file (default:no)-helpโ display help and exit-versionโ display version information and exit
so the tea is this ๐ซ: each sequence you hand over in sequence_file gets shreddered into consecutive pieces of random length (somewhere between -minlength and -maxlength) until it is fully consumed. which means โ and listen, this is the important part โ the last shreddered fragment of a given sequence can be shorter than the argument to option -minlength ๐ญ it's just how the math shakes out. if those scraggly little leftovers are giving you the ick, use gt seqfilter to get rid of them (there's an example down below, i got you ๐).
๐ EXAMPLES: ๐
Shredder a given BAC:
$ gt shredder U89959_genomic.fas > fragments.fas
Shredder an EST collection into pieces between 50 and 100 bp and get rid of all (terminal) fragments shorter than 50 bp:
$ gt shredder -minlength 50 -maxlength 100 U89959_ests.fas \
| gt seqfilter -minlength 50 - > fragments.fas
# 130 out of 1260 sequences have been removed (10.317%)
Shredder an EST collection and show only random 10% of the resulting fragments:
$ gt shredder -sample 0.1 U89959_ests.fas
๐ซถ REPORTING BUGS ๐ซถ
found a bug? not the vibe, we gotta know. report bugs to https://github.com/genometools/genometools/issues ๐ฎ
friendly reminder bestie: this is the GenomeTools1.6.6release, dated11/09/2025๐ so if you're reading this in some distant future era, hi, hope the shredder still slays ๐
go forth and shred, main character ๐ฆโจ๐
๐ฆ where this page comes from
- ubuntu resolute: genometools 1.6.6+ds-2
- ubuntu stonking: genometools 1.6.6+ds-2
- ๐ girlified by deepseek-ai/DeepSeek-V4.1-Flash on : a language model's rewrite of the original, under its terms