GT-SHREDDER(1)

๐Ÿ“ฆ
version 1.6.6 ยท 5 versions
, in ubuntu resolute, ubuntu stonking

๐Ÿ’– 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, \n etc) adds offset and length to ensure unique identifier (default: no)
  • -width [value] โ€” set output width for FASTA sequence printing (0 disables 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 GenomeTools 1.6.6 release, dated 11/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