/usr/share/genometools/gtdata/doc/shredder.lua is in genometools-common 1.5.1-2ubuntu1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 | print ([[
Each sequence given in sequence_file is shreddered into consecutive pieces of
random length (between -minlength and -maxlength) until it is consumed.
By this means the last shreddered fragment of a given sequence can be shorter
than the argument to option -minlength.
To get rid of such fragments use `gt seqfilter` (see example below).
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]])
|