By downloading this code you agree to the
Source Code Use License (PDF). |
{version = 4.27; (* of hist.p 1994 sep 6}
(* begin module describe.hist *)
(*
name
hist: make a histogram of aligned sequences.
synopsis
hist(inst: in, book: in, hst: out, histp: in, output: out);
files
inst: the instructions which generated the book, used to
align the sequences; the format of the instructions must be
'get from # -x to # +y;' - the alignment is done on the base #.
if this file is empty the sequences are aligned by their 5
prime ends.
book: the sequences, not longer than 'dnamax' (see technical note);
hst: the histogram table, giving the position occurence
of all oligonucleotides from oligomin to oligomax in length
(see file histp). note - if the histogram is wider than can
be printed on a page, use the program split to print hst.
histp: to set the length of oligonucleotides searched for; contains two
integers, one per line, the first specifying oligomin
(the length of the shortest oligonucleotide which is searched for),
the second oligomax (the longest oligonucleotide searched for,
which cannot be greater than 4);
note: if oligomin is zero, then the bases are counted.
output: for error messages.
description
makes a histogram of the occurences of oligonucleotides at positions
relative to some aligned base. this is done for all oligonucleotides
with lengths from 'oligomin' to 'oligomax', set by file histp.
see also
split.p, alist.p, histan.p
author
gary stormo
bugs
none known
technical note
the constant 'dnamax' from the module book.const may be too
large for efficient use of this program. if you do not expect to do
histograms on aligned sequences of greater than, say, 120 nucleotides
you can go into the source program (hists) and change 'dnamax' before
compiling.
*)
(* end module describe.hist *)
{This manual page was created by makman 1.45}
{created by htmlink 1.62}