table of contents
other versions
- jessie 1.5.3-2
- jessie-backports 1.5.9+ds-2~bpo8+1
- stretch 1.5.9+ds-4
- testing 1.5.10+ds-2
- unstable 1.5.10+ds-3
GT-SHREDDER(1) | GenomeTools Manual | GT-SHREDDER(1) |
NAME¶
gt-shredder - Shredder sequence file(s) into consecutive pieces of random length.SYNOPSIS¶
gt shredder [option ...] [sequence_file ...]DESCRIPTION¶
-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)
-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
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
$ 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%)
$ gt shredder -sample 0.1 U89959_ests.fas
REPORTING BUGS¶
Report bugs to <gt-users@genometools.org>.09/05/2014 | GenomeTools 1.5.3 |