Wise2 API (version 2.1.19b)

Ewan Birney
Sanger Centre
Wellcome Trust Genome Campus
Hinxton, Cambridge CB10 1SA,
England.
Email: birney@sanger.ac.uk

1 Overview

This document describes the API of the Wise2 system. The API (application programming interface) allows other programmers use the functionality in the Wise2 package directly, rather than treating the executables as a black box through which you get ASCII output.

If you want to learn more about the Wise2 package itself, the algorithms in it or what it is used for, look for the Wise2 documentation (available as postscript), probably in the same place that you found this documentation!

The API is accessible in 3 different ways: As a C function calls made inside the Wise2 package namespace - this is the way the current executables (eg, genewise) access the API, as C function calls made from outside the Wise2 package namespace - this is for people writing C programs with their own set of functions who do not want name clashes of things like “Sequence” (in this API the name is exported as “Wise2_Sequence”), and finally as a Perl API, using the XS extension code where C function calls which are dynamic loaded into the Perl interpretter can be executed as if they were standard Perl commands.

Probably the most usable is the Perl API. Perl is a very forgiving language, and it is easier to learn for novice programmers - in particular memory management is handled for you. For people who want to use the Wise2 api from inside their own C program, I would use the external api. For people who want to extend wise2 programs to do other things, the internal api.

2 WARNING - still in alpha

After playing around with the API for a while, I have realised that a number of things are not clean enough in the interface. I am not currently considering the API to the 2.1.x series stable. An aim for the 2.2 series is to make a stable and useful API, that is well documented.

However, this API does work, and there is this documentation for it, so it maybe worth people who like this sort of thing to play around with it. Anyone who uses the API gets huge guru points from me...

3 API generation

The API is not manually generated but rather is generated by the Dynamite compiler. Dynamite is a language which I wrote specifically for the Wise2 project: it is a cranky but useful language based heavily on C (it converts its source code to C), with a portion dedicated to dynamic programming code (a common algorithm in bioinformatics). It also has a lightweight object model that supports scalars and lists of types. Because the API is generated through Dynamite, you can expect consistent documentation and memory handling of all the functions and objects.

4 Getting Started for the impatient

Here is 3 different ways of using the Wise2 API to reverse complement a sequence. Once in perl, once using the name space protected API, and once using the internal API.

These three programs all make the same output, using the same code. It is only how the programming is presented to the user (once in perl, twice in C) which changes.

4.1 Perl reverse complement

#!/usr/local/bin/perl


use Wise2; # loads in Wise2 api

$file = shift; # first argument

if( !defined $file ) {
    print "You must give a file to revcom for a reverse to work!";
    exit(1);
}

$seq = &Wise2::Sequence::read_fasta_file_Sequence($file);
$rev = $seq->revcomp(); 

print "Original sequence\n\n";
$seq->write_fasta(STDOUT);
print "Reversed sequence\n\n";
$rev->write_fasta(STDOUT);

4.2 Wise2 external API calls

#include "dyna_api.h"


int main(int argc,char ** argv)
{
  Wise2_Sequence * seq;
  Wise2_Sequence * rev;

  if( argc != 2 ) {
    fprintf(stderr,"have to give an argument for a file");
    exit(1);
  }

  seq = Wise2_read_fasta_file_Sequence(argv[1]);

  if( seq == NULL ) {
    fprintf(stderr,"Unable to read fasta file in %s\n",argv[1]);
    exit(1);
  }
  
  rev = Wise2_reverse_complement_Sequence(seq);

  printf("Original sequence\n\n");
  Wise2_write_fasta_Sequence(seq,stdout);
  printf("Revcomp sequence\n\n");
  Wise2_write_fasta_Sequence(rev,stdout);
 
  Wise2_free_Sequence(seq);
  Wise2_free_Sequence(rev);
}

4.3 Wise2 internal API calls

#include "dyna.h"


int main(int argc,char ** argv)
{
  Sequence * seq;
  Sequence * rev;

  if( argc != 2 ) {
    fprintf(stderr,"have to give an argument for a file");
    exit(1);
  }

  seq = read_fasta_file_Sequence(argv[1]);

  if( seq == NULL ) {
    fprintf(stderr,"Unable to read fasta file in %s\n",argv[1]);
    exit(1);
  }
  
  rev = reverse_complement_Sequence(seq);

  printf("Original sequence\n\n");
  write_fasta_Sequence(seq,stdout);
  printf("Revcomp sequence\n\n");
  write_fasta_Sequence(rev,stdout);
 
  free_Sequence(seq);
  free_Sequence(rev);
}

5 Navigating the source code

The Wise2 api has a bewildering number of objects and functions, and the biggest problem in using the API is knowing which objects can be made from what. This next section walks you through at an object level how to do some common tasks. This list is in no way complete, but it is better than just browsing around the index.

A very good place to start is to read the scripts in the perl/scripts area (halfwise.pl does not use the Wise2 API but all the others do).

5.1 Making a translation of a DNA sequence

Build a codon table object from a file (28.2)
Build a sequence object, from a file or strings (46.2)
Use the translate function on the Sequence object

5.2 Comparing two sequences using smith waterman

Build a Comparison matrix object from a file (31.1)
Build two Sequence objects, from a file or strings (46.2)
Optionally convert the Sequence objects into Protein objects (42.1). This ensures you have proteins
Read in the comparison matrix using CompMat (31.1)
Use one of the algorithm calls in sw_wrap module (20)
Show the alignment using a call in the seqaligndisaply module (18)

5.3 Running a smith waterman search of single protein sequence vs a db

Read in a sequence object and convert it to a protein object (42.1,46.2)
Make a protein database from the single protein object (43.2)
Make a protein database from a single fasta file (43.2)
Using one of the calls to the sw_wrap module, make a Hscore object (20)
Show the Hscore object using a show function (39.2)
Retrieve individual protein objects from the database by taking out the DataEntry objects (39.4) and passing them into the ProteinDB object (43.2), giving you a protein object
optionally align them as in the above section

5.4 Running a genewise on a single protein vs a single DNA sequence

See the script genewise.pl in the distribution

Make a Sequence object from a strings or a file (46.2)
Make that a protein object (42.1)
Make a Sequence object from a string or a file (46.2)
Make that a Genomic object (35.2)
Add any additional repeat areas from external information to the genomic object
Read in a gene frequency counts (??)
Read in a codon table (28.2)
Make a random DNA model (45.2)
Make an algorithm type (??)
Build an entire parameter set for genewise using Wise2::GeneParameter21_wrap (??)
Run the actual algorithm (??)
show the alignment using genedisplay (??)

6 Concepts and overview of the API

The API is organised in the following way. There are 4 main areas of source code in the wise2 package

wisebase - base memory, string and error handling libraries
dynlibsrc - generic bioinformatics objects
models - specific Wise2 objects
HMMer2 - HMMER 2 (Sean Eddy’s HMM package)

The API is mainly derived from the dynlibsrc and models directories. There is no distinction in the API of one directory from another

7 The reference section

The reference section is built automatically from the Dynamite source. This means that the function names, argument lists and in nearly all cases the documentation should be completely up to date with whatever version you got this documentation from.

The code is divided up into modules: each module having potentially a number of objects in it and a number of free standing functions (factory functions). The documentation lists each object and the fields in the object which are accessible by the Perl API and the external API (more fields maybe accessible by the internal API, but generally these are not fields that you are expected to use). Fields can either be scalar or list types. In either case the scalar or list can either be a basic type or another object type. The following access methods are available for scalar types

7.1 Accessing fields in objects

In both the external API and the Perl API you can access all the fields via function calls. In Perl these function calls have the correct names space system to be called using the OOP syntax of Perl.

7.1.1 Perl scalar accessors

$obj>fieldname() gets the value of this field
$obj>set_fieldname(new value) sets the value of this field

For example

   $name = $seq->name();       # get the name of a sequence
   $seq->set_name('NewName');  # set the name of a sequence

7.1.2 External C scalar accessors

Wise2_access_fieldname_ObjectName(obj) gets the value of this field
Wise2_replace_fieldname_ObjectName(obj,new value) sets the value of this field

For example

   char * name;
   Wise2_Sequence * seq;

   /* ... get a sequence object somehow ... */
 
   name = Wise2_access_name_Sequence(seq);
   
   Wise2_replace_name_Sequence(seq,"NewName");

7.1.3 Perl List accessors

$obj>each_fieldname() Gives a Perl array of all the items in a list
$obj>length_fieldname Length of the list
$obj>fieldname($i) The ith member of the list
$obj>add_fieldname($another_obj) Adds another object to the list
$obj>flush_fieldname() Destroys all the items in a list, sets list size to zero

7.1.4 External API List accessors

Wise2_access_fieldname_ObjectName(obj,i) access the ith position in the list
Wise2_flush_fieldname_ObjectName(obj) Flushes the list
Wise2_add_fieldname_ObjectName(obj,added_object) Adds an object onto the end of the list

7.2 Object Construction and handling

The good news is that in the Perl API all the memory handling is managed between the Perl memory handling method and the Wise2 handling method. Bascially you can completely forget about these things and code normally in Perl and all the memory is handled for you.

In the C external API, as in any C program, the programmer is responsible for the memory, and you need to read the documentation as to whether the objects you recieve from function calls need explict frees or not.

7.2.1 Low level Object Constructing

In both Perl and in C you have the possibility of making a new object from scratch.

In Perl it is $obj = new Wise2::ObjectName;
in C it is obj = ObjectName_alloc for objects with no lists, and obj = ObjectName_alloc_std for objects with lists (this is a mistake in the API I know).

However I would read carefully the documentation for an object first, as in some cases the objects have to made through specific functions. These are likely to be things like new_ObjectName or such like. They are likely to be “factory” functions, that is functions not attached to any object.

7.2.2 Object deconstructors

In Perl you don’t have to worry about this (heaven).

In the C API you have two functions to handle the memory of objects. The objects have a reference counted memory: when the free function is called it decrements the object reference count and if this count hits 0 then the object itself is free’d. To up the reference count you call the hard_link_ObjectName function.

free_ObjectName(obj) Releases this pointer on this object
obj = hard_link_ObjectName() Adds this pointer to object, increasing the reference count

8 Wise2 Specific Modules

There are a number of modules which are specific to Wise2 algorithms. These should be the starting point for how to use Wise2 algorithms: try to find a function in these modules which provide the functionality that you want. Then figure out how to make the appropiate objects to use this functionality.

gwrap: ?? The gwrap module has the main entry points for the genewise algorithm and how to build parameters for it
estwrap: 12 The estwrap module has the main entry points for the estwise algorithm
sw_wrap: 20 The sw_wrap has the main entry points for the smith waterman algorithm
genedisplay: 13 The pretty ascii output used for genewise and estwise output
seqaligndisplay: 18 The pretty ascii output used for smith waterman alignments
threestatemodel: 22 profile-HMM support
threestatedb: 21 profile-HMM database support
genefrequency: 14 Raw counts for the genewise model
geneparameter: 15 probabilities for the genewise model
cdparser: 10 probabilities for the estwise model

9 Dynamite library modules

9.1 Sequence modules

sequence: 46 Basic sequences
sequencedb: 48 Basic sequence database
protein: 42 Protein specific type
proteindb: 43 Protein database
genomic: 35 Genomic specific type
genomicdb: 36 Genomic database
cdna: 26 Cdna specific type
cdnadb: 27 Cdna database

9.2 Generic probabilistic modelling support

probability: 41 Probability to log space conversions
codon: 28 Codon Table support
compmat: 31 Protein Comparison matrix support
codonmat: ?? Codon Matrix comparison matrix
codonmapper: 29 Codon bias/substitution errors support

9.3 Generic Database Searching

hscore: 39 High Score list
histogram: 38 Extreme Value distribution fitting
dbimpl: ?? Database Implementation

9.4 Generic Dynamite algorithm support

aln: 23 Label alignments
packaln: 40 Raw (low level) alignments
basematrix: 25 Memory management for DP implementations

10 cdparser

This module contains the following objects

10.2 cDNAParser
This module also contains some factory methods

10.1 cdparser factory methods

10.1.1 flat_cDNAParser

External C: Wise2_flat_cDNAParser (p)
Perl: &Wise2::flat_cDNAParser (p)

Arguments

p: [READ ] probability of an indel [Probability]
returns: [UNKN ] Undocumented return value [cDNAParser *]

Makes a flat (ie, indels of 1 or 2 == p) cDNA parser. This means that insertions and deletions of both 1 or 2 bases are all parameterised at the same probability

10.2 Object cDNAParser

The cDNAParser object has the following fields. To see how to access them refer to 7.1

: trans[PCD_PARSER_TRANS_LEN] Type [Probability : Scalar] No documentation

This object holds the (very few) extra transition information needed for the estwise algorithm. It is sort of like the ǵene modelṕart of sequencing error (but very very simple)

Member functions of cDNAParser

11 dnaalign

This module contains the following objects

11.2 DnaStartEnd
This module also contains some factory methods

11.1 dnaalign factory methods

11.1.1 make_align_dnaalign

External C: Wise2_make_align_dnaalign (one,two,mat,se,qgap,qext,tgap,text,dpri)
Perl: &Wise2::make_align_dnaalign (one,two,mat,se,qgap,qext,tgap,text,dpri)

Arguments

one: [READ ] first sequence to align [Sequence *]
two: [READ ] second sequence to align [Sequence *]
mat: [READ ] DnaMatrix for the matching [DnaMatrix *]
se: [READ ] DnaStartEnd policy [DnaStartEnd *]
qgap: [READ ] gap open penalty in query (one) coordinate [int]
qext: [READ ] gap extension penalty in query (one) coordinate [int]
tgap: [READ ] gap open penalty in target (two) coordinate [int]
text: [READ ] gap extension penalty in target (two) coordinate [int]
dpri: [READ ] DPRunImpl structure [DPRunImpl *]
returns: [UNKN ] an alb structure of the alignment [AlnBlock *]

Makes an alignment out of two DNA sequences

11.1.2 DnaStartEnd_from_policy

External C: Wise2_DnaStartEnd_from_policy (policy)
Perl: &Wise2::DnaStartEnd_from_policy (policy)

Arguments

policy: [UNKN ] Undocumented argument [char *]
returns: [UNKN ] Undocumented return value [DnaStartEnd *]

Makes a DnaStartEnd from a particular string. Possible strings are:

local - fully local global - fully global edge - aligns only to edges

11.2 Object DnaStartEnd

The DnaStartEnd object has the following fields. To see how to access them refer to 7.1

: trans[DSE_NUMBER] Type [int : Scalar] start/end points possibilities

No documentation for DnaStartEnd

Member functions of DnaStartEnd

12 estwrap

This module only contains factory methods

12.1 estwrap factory methods

12.1.1 Hscore_from_TSM_estwise

External C: Wise2_Hscore_from_TSM_estwise (tdb,cdb,cp,cm,rmd,use_syn,alg,bits_cutoff,allN,flat_insert,report_level,die_on_error,dbsi)
Perl: &Wise2::Hscore_from_TSM_estwise (tdb,cdb,cp,cm,rmd,use_syn,alg,bits_cutoff,allN,flat_insert,report_level,die_on_error,dbsi)

Arguments

tdb: [READ ] a three state model database [ThreeStateDB *]
cdb: [READ ] a dna sequence database [cDNADB *]
cp: [READ ] the codon parser for this comparison [cDNAParser *]
cm: [READ ] the codon mapper for this comparison [CodonMapper *]
rmd: [READ ] random model used for the dna sequence comparison [RandomModelDNA *]
use_syn: [UNKN ] whether a synchronous coding model should be used or not [boolean]
alg: [UNKN ] algorithm to use [int]
bits_cutoff: [UNKN ] Undocumented argument [double]
allN: [UNKN ] Undocumented argument [Probability]
flat_insert: [UNKN ] Undocumented argument [boolean]
report_level: [UNKN ] Undocumented argument [int]
die_on_error: [UNKN ] if true, dies if there is an error [boolean]
dbsi: [UNKN ] Undocumented argument [DBSearchImpl *]
returns: [OWNER] a newly allocated Hscore structure of the search [Hscore *]

Runs a database search for the estwise set of algorithms

12.1.2 AlnBlock_from_Protein_estwise_wrap

External C: Wise2_AlnBlock_from_Protein_estwise_wrap (pro,cdna,cp,cm,ct,comp,gap,ext,is_global,rmd,alg,rm,use_syn,allN,dpri,palpoi)
Perl: &Wise2::AlnBlock_from_Protein_estwise_wrap (pro,cdna,cp,cm,ct,comp,gap,ext,is_global,rmd,alg,rm,use_syn,allN,dpri)

Arguments

pro: [READ ] protein to be used in the comparison [Protein *]
cdna: [READ ] cdna to be compared to [cDNA *]
cp: [READ ] cdna parser indicating insertion deletion probabilities [cDNAParser *]
cm: [READ ] codon mapper indicating substitution errors etc [CodonMapper *]
ct: [READ ] codon table for the codon->amino acid mappings [CodonTable *]
comp: [READ ] comparison matrix to use [CompMat *]
gap: [UNKN ] gap penalty [int]
ext: [UNKN ] extension penalty [int]
is_global: [UNKN ] if true, global start-end in protein is used [boolean]
rmd: [UNKN ] random model of dna to use [RandomModelDNA *]
alg: [UNKN ] est algorithm type to use [int]
rm: [UNKN ] random protein model for use with compmat [RandomModel *]
use_syn: [UNKN ] if true, uses a synchronous coding model [boolean]
allN: [UNKN ] Undocumented argument [Probability]
dpri: [UNKN ] Undocumented argument [DPRunImpl *]
palpoi: only for C api [WRITE] the raw packed alignment output if wanted [PackAln **]
returns: [UNKN ] Undocumented return value [AlnBlock *]

This function is the guts for the est single alignment mode.

It uses /AlnBlock_from_TSM_estwise_wrap for the heavy part of the call

12.1.3 AlnBlock_from_TSM_estwise_wrap

External C: Wise2_AlnBlock_from_TSM_estwise_wrap (tsm,cdna,cp,cm,ct,rmd,alg,use_syn,force_flat_insert,allN,dpri,palpoi)
Perl: &Wise2::AlnBlock_from_TSM_estwise_wrap (tsm,cdna,cp,cm,ct,rmd,alg,use_syn,force_flat_insert,allN,dpri)

Arguments

tsm: [READ ] threestatemodel to be compared to the dna [ThreeStateModel *]
cdna: [READ ] cdna to be compared to [cDNA *]
cp: [READ ] cdna parser indicating insertion deletion probabilities [cDNAParser *]
cm: [READ ] codon mapper indicating substitution errors etc [CodonMapper *]
ct: [READ ] codon table for the codon->amino acid mappings [CodonTable *]
rmd: [UNKN ] random model of dna to use [RandomModelDNA *]
alg: [UNKN ] est algorithm type to use [int]
use_syn: [UNKN ] if true, uses a synchronous coding model [boolean]
force_flat_insert: [UNKN ] Undocumented argument [boolean]
allN: [UNKN ] Undocumented argument [Probability]
dpri: [UNKN ] Undocumented argument [DPRunImpl *]
palpoi: only for C api [WRITE] the raw packed alignment output if wanted [PackAln **]
returns: [UNKN ] Undocumented return value [AlnBlock *]

This function is the basic wrap for Protein models vs cDNA sequences.

12.1.4 alg_estwrap_from_string

External C: Wise2_alg_estwrap_from_string (str)
Perl: &Wise2::alg_estwrap_from_string (str)

Arguments

str: [UNKN ] Undocumented argument [char *]
returns: [UNKN ] Undocumented return value [int]

This function returns the algorithm type for an est search from the string

13 genedisplay

This module only contains factory methods

13.1 genedisplay factory methods

13.1.1 protein2genomic_ascii_display

External C: Wise2_protein2genomic_ascii_display (alb,p,gen,ct,name,main,ofp)
Perl: &Wise2::protein2genomic_ascii_display (alb,p,gen,ct,name,main,ofp)

Arguments

alb: [UNKN ] logical alignment [AlnBlock *]
p: [UNKN ] protein sequence [Protein *]
gen: [UNKN ] genomic dna to do the comparison [Genomic *]
ct: [UNKN ] codon table for translation [CodonTable *]
name: [UNKN ] length of name block [int]
main: [UNKN ] length of main block [int]
ofp: [UNKN ] output file [FILE *]
returns: [UNKN ] Undocumented return value [boolean]

shows the alignment in alb between protsequence and protname with genomic into ofp with pretty formatting

13.1.2 protcdna_ascii_display

External C: Wise2_protcdna_ascii_display (alb,protsequence,protname,protoff,cdna,ct,name,main,mult,ofp)
Perl: &Wise2::protcdna_ascii_display (alb,protsequence,protname,protoff,cdna,ct,name,main,mult,ofp)

Arguments

alb: [UNKN ] logical alignment [AlnBlock *]
protsequence: [UNKN ] protein sequence - either real or an artifical consensus [char *]
protname: [UNKN ] name of the protein [char *]
protoff: [UNKN ] offset of the alb from the protein [int]
cdna: [UNKN ] cdna of the match [cDNA *]
ct: [UNKN ] codon table for translation [CodonTable *]
name: [UNKN ] length of name block [int]
main: [UNKN ] length of main block [int]
mult: [UNKN ] is multi-match [boolean]
ofp: [UNKN ] output file [FILE *]
returns: [UNKN ] Undocumented return value [boolean]

shows the alignment in alb between protsequence and protname with cdna into ofp with pretty formatting

14 genefrequency

This module contains the following objects

14.2 GeneFrequency21
14.3 GeneConsensus
14.4 GeneSingleCons
This module also contains some factory methods

14.1 genefrequency factory methods

14.1.1 read_GeneFrequency21_file

External C: Wise2_read_GeneFrequency21_file (filename)
Perl: &Wise2::read_GeneFrequency21_file (filename)

Arguments

filename: [UNKN ] will open from WISECONFIGDIR etc via openfile [char *]
returns: [UNKN ] a newly allocated structure [GeneFrequency21 *]

Opens the file with /openfile

Reads in a GeneFrequency (Mor-Ewan style)

14.1.2 read_GeneFrequency21

External C: Wise2_read_GeneFrequency21 (ifp)
Perl: &Wise2::read_GeneFrequency21 (ifp)

Arguments

ifp: [UNKN ] file pointer [FILE *]
returns: [UNKN ] a newly allocated structure [GeneFrequency21 *]

Reads in a GeneFrequency (Mor-Ewan style) file from ifp

14.2 Object GeneFrequency21

The GeneFrequency21 object has the following fields. To see how to access them refer to 7.1

: ss5 Type [GeneConsensus * : Scalar] No documentation
: ss3 Type [GeneConsensus * : Scalar] No documentation
: codon[64] Type [double : Scalar] No documentation
: central[4] Type [double : Scalar] No documentation
: py[4] Type [double : Scalar] No documentation
: spacer[4] Type [double : Scalar] No documentation
: transition[GENEFREQUENCY21_TRANSITION_LEN] Type [double : Scalar] No documentation
: cds_triplet[64] Type [double : Scalar] phase 0

No documentation for GeneFrequency21

Member functions of GeneFrequency21

14.3 Object GeneConsensus

The GeneConsensus object has the following fields. To see how to access them refer to 7.1

: center Type [int : Scalar] No documentation
: gsc Type [GeneSingleCons ** : List] No documentation

No documentation for GeneConsensus

Member functions of GeneConsensus

14.4 Object GeneSingleCons

The GeneSingleCons object has the following fields. To see how to access them refer to 7.1

: string Type [char * : Scalar] No documentation
: number Type [double : Scalar] No documentation

No documentation for GeneSingleCons

Member functions of GeneSingleCons

15 geneparameter

This module contains the following objects

15.1 GeneParameter21

15.1 Object GeneParameter21

The GeneParameter21 object has the following fields. To see how to access them refer to 7.1

: gp Type [GeneParser21 * : Scalar] No documentation
: cm Type [CodonMapper * : Scalar] No documentation
: cses Type [ComplexSequenceEvalSet * : Scalar] No documentation
: ss Type [SpliceSiteModel ** : List] held only to be freed when GeneParser21Set is freed
: rc Type [RandomCodon * : Scalar] needed to soak up the odd-and-sods of genes
: gwcm Type [GeneWiseCodonModel * : Scalar] No documentation
: ct Type [CodonTable * : Scalar] No documentation
: modelled_splice Type [boolean : Scalar] so we can alter balance scores.
: gms Type [GeneModel * : Scalar] No documentation

GeneParameter21 keeps all the parameters for genewise algorithms in one tidy unit.

This is also the switch between the old (compugen handled) and new statistics. This object can be made from either the old or the new statistics

I have made the object complete opaque to scripting apis because the contents have to be coordinated quite well

Member functions of GeneParameter21

16 matchsum

This module contains the following objects

16.1 MatchSummarySet
16.2 MatchSummary

16.1 Object MatchSummarySet

The MatchSummarySet object has the following fields. To see how to access them refer to 7.1

: ms Type [MatchSummary ** : List] No documentation

This holds a set of MatchSummaries,

Member functions of MatchSummarySet

16.1.1 MatchSummarySet_from_AlnBlock_estwise

External C: Wise2_MatchSummarySet_from_AlnBlock_estwise (alb,qname,offset,target)
Perl: &Wise2::MatchSummarySet::MatchSummarySet_from_AlnBlock_estwise (alb,qname,offset,target)
Perl-OOP call: $obj->MatchSummarySet_from_AlnBlock_estwise(qname,offset,target)

Arguments

alb: [UNKN ] Undocumented argument [AlnBlock *]
qname: [UNKN ] Undocumented argument [char *]
offset: [UNKN ] Undocumented argument [int]
target: [UNKN ] Undocumented argument [Sequence *]
returns: [UNKN ] Undocumented return value [MatchSummarySet *]

Builds a MatchSummarySet from a EstWise alignment. this makes alot of assumptions about the labels setc in alb, so make sure it was a estwise alignment - however as you can notice this is exactly the same labels as found in genewise set

16.1.2 MatchSummarySet_from_AlnBlock_genewise

External C: Wise2_MatchSummarySet_from_AlnBlock_genewise (alb,qname,protoff,target)
Perl: &Wise2::MatchSummarySet::MatchSummarySet_from_AlnBlock_genewise (alb,qname,protoff,target)
Perl-OOP call: $obj->MatchSummarySet_from_AlnBlock_genewise(qname,protoff,target)

Arguments

alb: [UNKN ] Undocumented argument [AlnBlock *]
qname: [UNKN ] Undocumented argument [char *]
protoff: [UNKN ] Undocumented argument [int]
target: [UNKN ] Undocumented argument [Sequence *]
returns: [UNKN ] Undocumented return value [MatchSummarySet *]

Builds a MatchSummarySet from a GeneWise alignment. this makes alot of assumptions about the labels setc in alb, so make sure it was a genewise alignment

16.2 Object MatchSummary

The MatchSummary object has the following fields. To see how to access them refer to 7.1

: bits Type [double : Scalar] No documentation
: qname Type [char * : Scalar] No documentation
: tname Type [char * : Scalar] No documentation
: qstart Type [int : Scalar] No documentation
: qend Type [int : Scalar] No documentation
: tstart Type [int : Scalar] No documentation
: tend Type [int : Scalar] No documentation
: qintron Type [int : Scalar] No documentation
: qframeshift Type [int : Scalar] No documentation
: tintron Type [int : Scalar] No documentation
: tframeshift Type [int : Scalar] No documentation

A Match Summary has summary statistics for a single alignment, with the two start/end ranges and the number of introns and frameshifts for each sequence (obviously, if one is a protein then neither are valid!)

Member functions of MatchSummary

17 pfamhmmer1db

This module contains the following objects

17.1 PfamHmmer1DB
17.2 PfamHmmer1Entry

17.1 Object PfamHmmer1DB

The PfamHmmer1DB object has the following fields. To see how to access them refer to 7.1

: en Type [PfamHmmer1Entry ** : List] No documentation
: dirname Type [char * : Scalar] directory name with the models
: cur Type [int : Scalar] No documentation
: def Type [RandomModel * : Scalar] default random model

Pfam Hmmer1db is a wrapper around a Pfam Hmmer database. This is file called HMM.s in a directory which has appropiate .HMM and .random files

Although this DB will be phased out, it is still around for a while.

This wont be used directly, but rather in a threestatedb model. It does not implement a full dynamite style db stream. rather it expects threestatedb to prod it in the correct manner.

Member functions of PfamHmmer1DB

17.2 Object PfamHmmer1Entry

The PfamHmmer1Entry object has the following fields. To see how to access them refer to 7.1

: entryname Type [char * : Scalar] No documentation
: is_random Type [boolean : Scalar] No documentation
: is_hmmls Type [boolean : Scalar] No documentation
: bits_cutoff Type [double : Scalar] No documentation

No documentation for PfamHmmer1Entry

Member functions of PfamHmmer1Entry

18 seqaligndisplay

This module only contains factory methods

18.1 seqaligndisplay factory methods

18.1.1 write_pretty_str_align

External C: Wise2_write_pretty_str_align (alb,qname,query,tname,target,name,main,ofp)
Perl: &Wise2::write_pretty_str_align (alb,qname,query,tname,target,name,main,ofp)

Arguments

alb: [UNKN ] alignment structure [AlnBlock *]
qname: [UNKN ] name of first sequence [char *]
query: [UNKN ] first sequence [char *]
tname: [UNKN ] name of second sequence [char *]
target: [UNKN ] second sequence [char *]
name: [UNKN ] length of the name block [int]
main: [UNKN ] length of the main block [int]
ofp: [UNKN ] output file [FILE *]
returns: [UNKN ] Undocumented return value [boolean]

This gives an interface into the alignment display using strings and files.

18.1.2 write_pretty_seq_align

External C: Wise2_write_pretty_seq_align (alb,q,t,name,main,ofp)
Perl: &Wise2::write_pretty_seq_align (alb,q,t,name,main,ofp)

Arguments

alb: [UNKN ] alignment structure [AlnBlock *]
q: [UNKN ] first sequence [Sequence *]
t: [UNKN ] second sequence [Sequence *]
name: [UNKN ] length of the name block [int]
main: [UNKN ] length of the main block [int]
ofp: [UNKN ] output file [FILE *]
returns: [UNKN ] Undocumented return value [boolean]

This gives an interface into the alignment display using sequences and files. A more generic function is write_pretty_str_align

18.1.3 write_pretty_Protein_align

External C: Wise2_write_pretty_Protein_align (alb,q,t,name,main,ofp)
Perl: &Wise2::write_pretty_Protein_align (alb,q,t,name,main,ofp)

Arguments

alb: [UNKN ] alignment structure [AlnBlock *]
q: [UNKN ] first sequence [Protein *]
t: [UNKN ] second sequence [Protein *]
name: [UNKN ] length of the name block [int]
main: [UNKN ] length of the main block [int]
ofp: [UNKN ] output file [FILE *]
returns: [UNKN ] Undocumented return value [boolean]

This gives an interface into the alignment display using Protein objects

19 seqhit

This module contains the following objects

19.1 DnaSequenceHitList
19.2 SegmentHitList
19.3 SegmentHit

19.1 Object DnaSequenceHitList

The DnaSequenceHitList object has the following fields. To see how to access them refer to 7.1

: forward Type [SegmentHitList * : Scalar] No documentation
: backward Type [SegmentHitList * : Scalar] No documentation

No documentation for DnaSequenceHitList

Member functions of DnaSequenceHitList

19.1.1 show_DnaSequenceHitList

External C: Wise2_show_DnaSequenceHitList (dsl,ofp)
Perl: &Wise2::DnaSequenceHitList::show_DnaSequenceHitList (dsl,ofp)
Perl-OOP call: $obj->show_DnaSequenceHitList(ofp)

Arguments

dsl: [UNKN ] Undocumented argument [DnaSequenceHitList *]
ofp: [UNKN ] Undocumented argument [FILE *]
returns: Nothing - no return value

shows a DnaSequenceHitsList -

only really useful for debugging

19.1.2 read_MSPcrunch_DnaSequenceHitList

External C: Wise2_read_MSPcrunch_DnaSequenceHitList (ifp)
Perl: &Wise2::DnaSequenceHitList::read_MSPcrunch_DnaSequenceHitList (ifp)
Perl-OOP call: $obj->read_MSPcrunch_DnaSequenceHitList()

Arguments

ifp: [UNKN ] input file to read [FILE *]
returns: [UNKN ] newly allocated structure [DnaSequenceHitList *]

Reads a MSPcrunch -x output file

19.2 Object SegmentHitList

The SegmentHitList object has the following fields. To see how to access them refer to 7.1

: seghit Type [SegmentHit ** : List] No documentation

No documentation for SegmentHitList

Member functions of SegmentHitList

19.3 Object SegmentHit

The SegmentHit object has the following fields. To see how to access them refer to 7.1

: name Type [char * : Scalar] No documentation
: qstart Type [int : Scalar] No documentation
: qend Type [int : Scalar] No documentation
: tstart Type [int : Scalar] No documentation
: tend Type [int : Scalar] No documentation
: score Type [double : Scalar] No documentation
: next_hit Type [SegmentHit * : Scalar] for collecting hits together ;)

No documentation for SegmentHit

Member functions of SegmentHit

20 sw_wrap

This module only contains factory methods

20.1 sw_wrap factory methods

20.1.1 Align_strings_ProteinSmithWaterman

External C: Wise2_Align_strings_ProteinSmithWaterman (one,two,comp,gap,ext,dpenv,dpri)
Perl: &Wise2::Align_strings_ProteinSmithWaterman (one,two,comp,gap,ext,dpenv,dpri)

Arguments

one: [UNKN ] string of the first sequence [char *]
two: [UNKN ] string of the second sequence [char *]
comp: [UNKN ] Comparison Matrix [CompMat *]
gap: [UNKN ] gap penalty [int]
ext: [UNKN ] extension penalty [int]
dpenv: [UNKN ] Undocumented argument [DPEnvelope *]
dpri: [UNKN ] Undocumented argument [DPRunImpl *]
returns: [UNKN ] Undocumented return value [AlnBlock *]

This is the most *stupidly* abstracted view of two sequences getting aligned, being two strings.

It would be much better if you used Sequence objects or Protein objects to carry the proteins.

20.1.2 Align_Sequences_ProteinSmithWaterman

External C: Wise2_Align_Sequences_ProteinSmithWaterman (one,two,comp,gap,ext,dpenv,dpri)
Perl: &Wise2::Align_Sequences_ProteinSmithWaterman (one,two,comp,gap,ext,dpenv,dpri)

Arguments

one: [READ ] First sequence to compare [Sequence *]
two: [READ ] Second sequecne to compare [Sequence *]
comp: [READ ] Comparison matrix to use [CompMat *]
gap: [UNKN ] gap penalty. Must be negative or 0 [int]
ext: [UNKN ] ext penalty. Must be negative or 0 [int]
dpenv: [UNKN ] Undocumented argument [DPEnvelope *]
dpri: [UNKN ] Undocumented argument [DPRunImpl *]
returns: [OWNER] new AlnBlock structure representing the alignment [AlnBlock *]

This function is a mid-level abstraction of comparing two sequences, which could be generic types (eg DNA!). This is tested for and warnings are given but the alignment is still calculated. To prevent this test warning either make sure the Sequence types are set to PROTEIN or, better still, use the high level abstraction Align_Proteins_SmithWaterman

Otherwise this performs a standard smith waterman protein alignment...

To display the alignment use write_pretty_seq_align

20.1.3 Align_Proteins_SmithWaterman

External C: Wise2_Align_Proteins_SmithWaterman (one,two,comp,gap,ext,dpenv,dpri)
Perl: &Wise2::Align_Proteins_SmithWaterman (one,two,comp,gap,ext,dpenv,dpri)

Arguments

one: [UNKN ] Protein to align [Protein *]
two: [UNKN ] Protein to align [Protein *]
comp: [UNKN ] Comparison Matrix [CompMat *]
gap: [UNKN ] gap penalty [int]
ext: [UNKN ] extension penalty [int]
dpenv: [UNKN ] Undocumented argument [DPEnvelope *]
dpri: [UNKN ] Undocumented argument [DPRunImpl *]
returns: [UNKN ] Undocumented return value [AlnBlock *]

This is the most correct way of aligning two Proteins, using Protein objects, which can be assummed to be proteins with no objections

To display the alignment use write_pretty_Protein_align

20.1.4 Align_Proteins_ABC

External C: Wise2_Align_Proteins_ABC (one,two,comp,a,b,c,dpenv,dpri)
Perl: &Wise2::Align_Proteins_ABC (one,two,comp,a,b,c,dpenv,dpri)

Arguments

one: [UNKN ] protein to align [Protein *]
two: [UNKN ] protein to align [Protein *]
comp: [UNKN ] comparison matrix [CompMat *]
a: [UNKN ] generalized affine gap cost a [int]
b: [UNKN ] generalized affine gap cost b [int]
c: [UNKN ] generalized affine gap cost c [int]
dpenv: [UNKN ] Undocumented argument [DPEnvelope *]
dpri: [UNKN ] Undocumented argument [DPRunImpl *]
returns: [UNKN ] Undocumented return value [AlnBlock *]

Analogous to Align_Proteins_SmithWaterman for ABC model

20.1.5 Align_Sequences_ProteinABC

External C: Wise2_Align_Sequences_ProteinABC (one,two,comp,a,b,c,dpenv,dpri)
Perl: &Wise2::Align_Sequences_ProteinABC (one,two,comp,a,b,c,dpenv,dpri)

Arguments

one: [UNKN ] Sequence to align [Sequence *]
two: [UNKN ] Sequence to align [Sequence *]
comp: [UNKN ] Comparison Matrix [CompMat *]
a: [UNKN ] genearlized affine gap cost [int]
b: [UNKN ] genearlized affine gap cost [int]
c: [UNKN ] genearlized affine gap cost [int]
dpenv: [UNKN ] Undocumented argument [DPEnvelope *]
dpri: [UNKN ] Undocumented argument [DPRunImpl *]
returns: [UNKN ] Undocumented return value [AlnBlock *]

Align_Sequences_ProteinABC this function is analogous to Align_Sequences_ProteinSmithWaterman but using the abc model

20.1.6 Hscore_from_ProteinSW

External C: Wise2_Hscore_from_ProteinSW (querydb,targetdb,comp,gap,ext,bits_cutoff,report_level,die_on_error,dbsi)
Perl: &Wise2::Hscore_from_ProteinSW (querydb,targetdb,comp,gap,ext,bits_cutoff,report_level,die_on_error,dbsi)

Arguments

querydb: [UNKN ] query database [ProteinDB*]
targetdb: [UNKN ] target database [ProteinDB*]
comp: [UNKN ] comparison matrix [CompMat*]
gap: [UNKN ] gap penalty [int]
ext: [UNKN ] extension penalty [int]
bits_cutoff: [UNKN ] [double]
report_level: [UNKN ] [int]
die_on_error: [UNKN ] [boolean]
dbsi: [UNKN ] [DBSearchImpl*]
returns: [UNKN ] Undocumented return value [Hscore *]

Runs a database psw search

20.1.7 Hscore_from_ProteinABC

External C: Wise2_Hscore_from_ProteinABC (querydb,targetdb,comp,a,b,c,bits_cutoff,report_level,die_on_error,dbsi)
Perl: &Wise2::Hscore_from_ProteinABC (querydb,targetdb,comp,a,b,c,bits_cutoff,report_level,die_on_error,dbsi)

Arguments

querydb: [UNKN ] query database [ProteinDB*]
targetdb: [UNKN ] target database [ProteinDB*]
comp: [UNKN ] comparison matrix [CompMat*]
a: [UNKN ] generalized affine gap cost a [int]
b: [UNKN ] generalized affine gap cost b [int]
c: [UNKN ] generalized affine gap cost c [int]
bits_cutoff: [UNKN ] [double]
report_level: [UNKN ] [int]
die_on_error: [UNKN ] [boolean]
dbsi: [UNKN ] [DBSearchImpl*]
returns: [UNKN ] Undocumented return value [Hscore *]

Runs a database abc search

20.1.8 Hscore_from_ProteinBA

External C: Wise2_Hscore_from_ProteinBA (querydb,targetdb,comp,bentry,bexit,bfor_trans,b_self_trans,b3exit,bits_cutoff,report_level,dbsi)
Perl: &Wise2::Hscore_from_ProteinBA (querydb,targetdb,comp,bentry,bexit,bfor_trans,b_self_trans,b3exit,bits_cutoff,report_level,dbsi)

Arguments

querydb: [UNKN ] query database [ProteinDB*]
targetdb: [UNKN ] target database [ProteinDB*]
comp: [UNKN ] comparison matrix [CompMat*]
bentry: [UNKN ] [Score]
bexit: [UNKN ] [Score]
bfor_trans: [UNKN ] [Score]
b_self_trans: [UNKN ] [Score]
b3exit: [UNKN ] [Score]
bits_cutoff: [UNKN ] [double]
report_level: [UNKN ] [int]
dbsi: [UNKN ] [DBSearchImpl*]
returns: [UNKN ] Undocumented return value [Hscore *]

Runs a database pba search

21 threestatedb

This module contains the following objects

21.1 ThreeStateDB

21.1 Object ThreeStateDB

The ThreeStateDB object has the following fields. To see how to access them refer to 7.1

: dbtype Type [int : Scalar] No documentation
: filename Type [char * : Scalar] No documentation
: type Type [int : Scalar] No documentation
: current_file Type [FILE * : Scalar] No documentation
: rm Type [RandomModel * : Scalar] NB, this is hard-linked...
: byte_position Type [long : Scalar] this is the file position for the current model
: single Type [ThreeStateModel * : Scalar] for single db cases
: phdb Type [PfamHmmer1DB * : Scalar] for Pfam Hmmer1 style databases.
: sdb Type [SequenceDB * : Scalar] for protein databases
: comp Type [CompMat * : Scalar] for protein databases
: gap Type [int : Scalar] for protein databases
: ext Type [int : Scalar] for protein databases
: seq_cache Type [Sequence * : Scalar] needed for a bit of inter-function communication
: reload_generic Type [ThreeStateModel * (*reload_generic)(ThreeStateDB * tdb,int * return_status) : Scalar] for generic parsing applications...
: open_generic Type [boolean (*open_generic)(ThreeStateDB * tdb) : Scalar] No documentation
: close_generic Type [boolean (*close_generic)(ThreeStateDB * tdb) : Scalar] No documentation
: dataentry_add Type [boolean (*dataentry_add)(ThreeStateDB * tdb,DataEntry * en) : Scalar] No documentation
: open_index_generic Type [boolean (*open_index_generic)(ThreeStateDB *tdb) : Scalar] No documentation
: index_generic Type [ThreeStateModel * (*index_generic)(ThreeStateDB *tdb,DataEntry *de) : Scalar] No documentation
: close_index_generic Type [boolean (*close_index_generic)(ThreeStateDB *tdb) : Scalar] No documentation
: data Type [void * : Scalar] whatever else the damn system wants to carry around with it!
: hmm_model_start Type [int : Scalar] No documentation
: hmm_model_end Type [int : Scalar] No documentation
: current_no Type [int : Scalar] No documentation

ThreeStateDB is the object that represents a database of profile-HMMs.

The object hold a variety of fields on some of which are occupied depending on the type.

Realistically we need a more abstract class idea, which is implemented here anyway via the generic stuff, in hacky C-style pointers to function plus a void pointer. This object therefore houses a switch system around the different types including the generic system... but as the generic function stuff was bolted on later, some things are handled with explicit datastructures. It is quite messy ;). Apologies. To be cleaned up.

The generic stuff was principly added in to allow a decoupling of this module from the HMMer2.o interface code which is held in wise2xhmmer.dy

The old static datastructure code can be made via protein sequences which are then converted or a Pfam 2.0 style directory + HMMs file.

Member functions of ThreeStateDB

21.1.1 indexed_ThreeStateModel_ThreeStateDB

External C: Wise2_indexed_ThreeStateModel_ThreeStateDB (mdb,en)
Perl: &Wise2::ThreeStateDB::indexed_model (mdb,en)
Perl-OOP call: $obj->indexed_model(en)

Arguments

mdb: [UNKN ] database where this is indexed [ThreeStateDB *]
en: [UNKN ] dataentry to pull the model from [DataEntry *]
returns: [UNKN ] Undocumented return value [ThreeStateModel *]

Retrieves a model from a database which has been opened for indexing by /open_for_indexing_ThreeStateDB

The index information comes from the dataentry which should have been from a search of the ThreeStateDB.

21.1.2 new_proteindb_ThreeStateDB

External C: Wise2_new_proteindb_ThreeStateDB (sdb,comp,gap,ext)
Perl: &Wise2::ThreeStateDB::new_proteindb_ThreeStateDB (sdb,comp,gap,ext)
Perl-OOP call: $obj->new_proteindb_ThreeStateDB(comp,gap,ext)

Arguments

sdb: [READ ] sequence database to use [SequenceDB *]
comp: [READ ] comparison matrix to use [CompMat *]
gap: [READ ] gap open penalty [int]
ext: [READ ] gap extensions penalty [int]
returns: [UNKN ] Undocumented return value [ThreeStateDB *]

makes a new ThreeStateDB from a sequencedb (better be protein!)

21.1.3 new_PfamHmmer1DB_ThreeStateDB

External C: Wise2_new_PfamHmmer1DB_ThreeStateDB (dirname)
Perl: &Wise2::ThreeStateDB::new_PfamHmmer1DB_ThreeStateDB (dirname)
Perl-OOP call: $obj->new_PfamHmmer1DB_ThreeStateDB()

Arguments

dirname: [UNKN ] Undocumented argument [char *]
returns: [UNKN ] Undocumented return value [ThreeStateDB *]

Makes a new PfamHmmer1DB from a filename indicating the directory

21.1.4 new_single_ThreeStateDB

External C: Wise2_new_single_ThreeStateDB (tsm,rm)
Perl: &Wise2::ThreeStateDB::new_single_ThreeStateDB (tsm,rm)
Perl-OOP call: $obj->new_single_ThreeStateDB(rm)

Arguments

tsm: [READ ] a single ThreeStateModel [ThreeStateModel *]
rm: [READ ] random model to be used in comparisons.. [RandomModel *]
returns: [UNKN ] Undocumented return value [ThreeStateDB *]

Making a new ThreeStateDB from a single model

22 threestatemodel

This module contains the following objects

22.2 ThreeStateModel
22.3 ThreeStateUnit
This module also contains some factory methods

22.1 threestatemodel factory methods

22.1.1 read_HMMer_1_7_ascii_file

External C: Wise2_read_HMMer_1_7_ascii_file (filename)
Perl: &Wise2::read_HMMer_1_7_ascii_file (filename)

Arguments

filename: [UNKN ] the name fo the hmmer file [char *]
returns: [UNKN ] Undocumented return value [ThreeStateModel *]

reads a HMMer ascii version 1.7 (1.8) file from filename.

22.1.2 read_HMMer_1_7_ascii

External C: Wise2_read_HMMer_1_7_ascii (ifp)
Perl: &Wise2::read_HMMer_1_7_ascii (ifp)

Arguments

ifp: [UNKN ] Undocumented argument [FILE *]
returns: [UNKN ] Undocumented return value [ThreeStateModel *]

Basic function to read HMMer version 1.7(1.8) files.

22.2 Object ThreeStateModel

The ThreeStateModel object has the following fields. To see how to access them refer to 7.1

: name Type [char * : Scalar] name of the model
: unit Type [ThreeStateUnit ** : List] the actuall three state probs and emissions
: alphabet Type [char * : Scalar] alphabet used
: accession Type [char * : Scalar] accession number
: threshold Type [double : Scalar] bits threshold (if sensible)
: rm Type [RandomModel * : Scalar] Random model for the model: maybe NULL!

This is profile-HMM object, similar to the SAM and HMMer plan9 architecture.

Member functions of ThreeStateModel

22.2.1 force_global_model

External C: Wise2_force_global_model (tsm,prob_into_model)
Perl: &Wise2::ThreeStateModel::force_global_model (tsm,prob_into_model)
Perl-OOP call: $obj->force_global_model(prob_into_model)

Arguments

tsm: [UNKN ] ThreeStateModel to be forced[ThreeStateModel *]
prob_into_model: [UNKN ] Probability to start the model: for true global will be 1.0 [double]
returns: Nothing - no return value

Makes start at position 0 and end at position end, no other positions being valid

22.2.2 force_weighted_local_model

External C: Wise2_force_weighted_local_model (tsm,prob_into_model,ratio_start,ratio_end)
Perl: &Wise2::ThreeStateModel::force_weighted_local_model (tsm,prob_into_model,ratio_start,ratio_end)
Perl-OOP call: $obj->force_weighted_local_model(prob_into_model,ratio_start,ratio_end)

Arguments

tsm: [UNKN ] ThreeStateModel to be forced[ThreeStateModel *]
prob_into_model: [UNKN ] Probability to start the model: for true global will be 1.0 [double]
ratio_start: [UNKN ] ratio of prob to unit 0 to the rest (1.0 means all goes to start) [double]
ratio_end: [UNKN ] ratio of prob to unit (last) to the rest (1.0 means all goes to the end) [double]
returns: Nothing - no return value

places the ratio of probability to start/end, and then distributes the rest over the start/end

22.2.3 ThreeStateModel_from_half_bit_Sequence

External C: Wise2_ThreeStateModel_from_half_bit_Sequence (pro,mat,rm,gap,ext)
Perl: &Wise2::ThreeStateModel::ThreeStateModel_from_half_bit_Sequence (pro,mat,rm,gap,ext)
Perl-OOP call: $obj->ThreeStateModel_from_half_bit_Sequence(mat,rm,gap,ext)

Arguments

pro: [READ ] protein sequence [Protein *]
mat: [READ ] comparison matrix to use [CompMat *]
rm: [READ ] random model which you assumme the matrix was built with [RandomModel *]
gap: [READ ] gap open penalty [int]
ext: [READ ] gap ext penalty [int]
returns: [UNKN ] Undocumented return value [ThreeStateModel *]

Makes a local three-state-model from a sequence. this is scary hackery, assumming that the matrix is half-bits and normalising in a *very* wrong way to get "probabilities" out.

Works though

22.2.4 write_HMMer_1_7_ascii_ThreeStateModel

External C: Wise2_write_HMMer_1_7_ascii_ThreeStateModel (tsm,ofp)
Perl: &Wise2::ThreeStateModel::write_HMMer_1_7_ascii_ThreeStateModel (tsm,ofp)
Perl-OOP call: $obj->write_HMMer_1_7_ascii_ThreeStateModel(ofp)

Arguments

tsm: [UNKN ] Undocumented argument [ThreeStateModel *]
ofp: [UNKN ] Undocumented argument [FILE *]
returns: Nothing - no return value

writes a HMMer version 1.7 (also ok with 1.8) file

22.3 Object ThreeStateUnit

The ThreeStateUnit object has the following fields. To see how to access them refer to 7.1

: match_emission[ALPHABET_SIZE] Type [Probability : Scalar] No documentation
: insert_emission[ALPHABET_SIZE] Type [Probability : Scalar] No documentation
: transition[TRANSITION_LEN] Type [Probability : Scalar] No documentation
: display_char Type [char : Scalar] No documentation

This object is the probability version of hte common unit to profile HMMs, ie the match,insert,delete triple

Member functions of ThreeStateUnit

23 aln

This module contains the following objects

23.1 AlnBlock
23.2 AlnColumn
23.3 AlnUnit
23.4 AlnSequence

23.1 Object AlnBlock

The AlnBlock object has the following fields. To see how to access them refer to 7.1

: start Type [AlnColumn * : Scalar] the first AlnColumn in the alignment
: seq Type [AlnSequence ** : List] a list of AlnSequences in the alignment
: length Type [int : Scalar] not used
: score Type [int : Scalar] not used

AlnBlock is the main representation of alignments from Dynamite. Each AlnBlock represents any number of śequences, of any type, which share things in common. The alignment is represented by a series of /AlnColumns (linked list) in which each AlnColumn has a series of AlnUnits, each Unit being a start/end/text_label triple. Alternatively, one can see each sequence in isolation, and not ask what it is aligned to, but rather what labels it has on it.

Member functions of AlnBlock

23.1.1 dump_ascii_AlnBlock

External C: Wise2_dump_ascii_AlnBlock (alb,ofp)
Perl: &Wise2::AlnBlock::dump_ascii_AlnBlock (alb,ofp)
Perl-OOP call: $obj->dump_ascii_AlnBlock(ofp)

Arguments

alb: [UNKN ] AlnBlock to dump [AlnBlock *]
ofp: [UNKN ] File stream to dump to [FILE *]
returns: Nothing - no return value

Dumps the alignment in rereadable ascii form.

Not really for human consumption

23.2 Object AlnColumn

The AlnColumn object has the following fields. To see how to access them refer to 7.1

: alu Type [AlnUnit ** : List] list of the AlnUnits in this column
: next Type [AlnColumn * : Scalar] the next AlnColumn in this block

This is a coupling of AlnUnits from different sequences. Each AlnUnit is meant to represent *the equivalent* piece of biological information in some sense (ie, they are alignmed), even though quite possibly they are very different types of information

Member functions of AlnColumn

23.2.1 at_end_AlnColumn

External C: Wise2_at_end_AlnColumn (alc)
Perl: &Wise2::AlnColumn::at_end (alc)
Perl-OOP call: $obj->at_end()

Arguments

alc: [READ ] AlnColumn [AlnColumn *]
returns: [UNKN ] Undocumented return value [boolean]

This tells you whether the AlnColumn is at the end without passing NULLś around

23.3 Object AlnUnit

The AlnUnit object has the following fields. To see how to access them refer to 7.1

: start Type [int : Scalar] start position in the sequence
: end Type [int : Scalar] end position in the sequence
: label Type [int : Scalar] not used
: text_label Type [char * : Scalar] text label of this position
: next Type [AlnUnit * : Scalar] next AlnUnit in this sequence
: score[AlnUnitSCORENUMBER] Type [int : Scalar] a series of scores for this position.
: in_column Type [boolean : Scalar] not used
: seq Type [AlnSequence * : Scalar] No documentation

This is the basic unit of the label alignment. It describes a single mark-up over one sequence: being a start, an end and a text_label.

Member functions of AlnUnit

23.3.1 bio_start_AlnUnit

External C: Wise2_bio_start_AlnUnit (alu)
Perl: &Wise2::AlnUnit::bio_start (alu)
Perl-OOP call: $obj->bio_start()

Arguments

alu: [UNKN ] Undocumented argument [AlnUnit *]
returns: [UNKN ] Undocumented return value [int]

Tells the bio-coordinate of the start point of this alnunit

23.3.2 bio_end_AlnUnit

External C: Wise2_bio_end_AlnUnit (alu)
Perl: &Wise2::AlnUnit::bio_end (alu)
Perl-OOP call: $obj->bio_end()

Arguments

alu: [UNKN ] Undocumented argument [AlnUnit *]
returns: [UNKN ] Undocumented return value [int]

Tells the bio-coordinate of the end point of this alnunit

23.4 Object AlnSequence

The AlnSequence object has the following fields. To see how to access them refer to 7.1

: start Type [AlnUnit * : Scalar] the first AlnUnit of this sequence
: data_type Type [int : Scalar] not used
: data Type [void * : Scalar] not used - dont use!
: bio_start Type [int : Scalar] start of this sequence in its bioćoordinates
: bio_end Type [int : Scalar] end of this sequence in its bioćoordinates

Each Sequence in an AlnBlock is represented by one of these, and in many ways this is an orthogonal way of looking at the alignment than the AlnColumns. If you look at the alignment from one AlnSequence you will just see the individual labels on this sequence

Member functions of AlnSequence

24 alnrange

This module contains the following objects

24.1 AlnRange
24.2 AlnRangeSet

24.1 Object AlnRange

The AlnRange object has the following fields. To see how to access them refer to 7.1

: starti Type [int : Scalar] No documentation
: startj Type [int : Scalar] No documentation
: startstate Type [int : Scalar] No documentation
: stopi Type [int : Scalar] No documentation
: stopj Type [int : Scalar] No documentation
: stopstate Type [int : Scalar] No documentation
: startscore Type [int : Scalar] No documentation
: stopscore Type [int : Scalar] No documentation

No documentation for AlnRange

Member functions of AlnRange

24.2 Object AlnRangeSet

The AlnRangeSet object has the following fields. To see how to access them refer to 7.1

: score Type [int : Scalar] over complete alignment
: alr Type [AlnRange ** : List] No documentation

No documentation for AlnRangeSet

Member functions of AlnRangeSet

25 basematrix

This module only contains factory methods

25.1 basematrix factory methods

25.1.1 change_max_BaseMatrix_kbytes

External C: Wise2_change_max_BaseMatrix_kbytes (new_kilo_number)
Perl: &Wise2::change_max_BaseMatrix_kbytes (new_kilo_number)

Arguments

new_kilo_number: [UNKN ] max kilobytes allowed [int]
returns: Nothing - no return value

This is to change, at run-time the maximum level of bytes basematrix *thinks* it can use. This number is *not* used for any actual calls to basematrix allocation: it is only used with /get_max_BaseMatrix_kbytes

25.1.2 get_max_BaseMatrix_kbytes

External C: Wise2_get_max_BaseMatrix_kbytes (void)
Perl: &Wise2::get_max_BaseMatrix_kbytes ()

Arguments

returns: [UNKN ] Undocumented return value [int]

returns the max. number of kilobytes suggested as a limited to BaseMatrix.

26 cdna

This module contains the following objects

26.1 cDNA

26.1 Object cDNA

The cDNA object has the following fields. To see how to access them refer to 7.1

: baseseq Type [Sequence * : Scalar] No documentation

No documentation for cDNA

Member functions of cDNA

26.1.1 truncate_cDNA

External C: Wise2_truncate_cDNA (cdna,start,stop)
Perl: &Wise2::cDNA::truncate_cDNA (cdna,start,stop)
Perl-OOP call: $obj->truncate_cDNA(start,stop)

Arguments

cdna: [READ ] cDNA that is truncated [cDNA *]
start: [UNKN ] Undocumented argument [int]
stop: [UNKN ] Undocumented argument [int]
returns: [UNKN ] Undocumented return value [cDNA *]

Truncates a cDNA sequence. Basically uses the /magic_trunc_Sequence function (of course!)

It does not alter cdna, rather it returns a new sequence with that truncation

26.1.2 read_fasta_file_cDNA

External C: Wise2_read_fasta_file_cDNA (filename)
Perl: &Wise2::cDNA::read_fasta_file_cDNA (filename)
Perl-OOP call: $obj->read_fasta_file_cDNA()

Arguments

filename: [UNKN ] filename to be opened and read [char *]
returns: [UNKN ] Undocumented return value [cDNA *]

Reads a fasta file assumming that it is cDNA. Will complain if it is not, and return NULL.

26.1.3 cDNA_name

External C: Wise2_cDNA_name (cdna)
Perl: &Wise2::cDNA::cDNA_name (cdna)
Perl-OOP call: $obj->cDNA_name()

Arguments

cdna: [UNKN ] Undocumented argument [cDNA *]
returns: [UNKN ] Undocumented return value [char *]

Returns the name of the cDNA

26.1.4 cDNA_length

External C: Wise2_cDNA_length (cdna)
Perl: &Wise2::cDNA::cDNA_length (cdna)
Perl-OOP call: $obj->cDNA_length()

Arguments

cdna: [UNKN ] Undocumented argument [cDNA *]
returns: [UNKN ] Undocumented return value [int]

Returns the length of the cDNA

26.1.5 cDNA_seqchar

External C: Wise2_cDNA_seqchar (cdna,pos)
Perl: &Wise2::cDNA::cDNA_seqchar (cdna,pos)
Perl-OOP call: $obj->cDNA_seqchar(pos)

Arguments

cdna: [UNKN ] cDNA [cDNA *]
pos: [UNKN ] position in cDNA to get char [int]
returns: [UNKN ] Undocumented return value [char]

Returns sequence character at this position.

26.1.6 cDNA_from_Sequence

External C: Wise2_cDNA_from_Sequence (seq)
Perl: &Wise2::cDNA::cDNA_from_Sequence (seq)
Perl-OOP call: $obj->cDNA_from_Sequence()

Arguments

seq: [OWNER] Sequence to make cDNA from [Sequence *]
returns: [UNKN ] Undocumented return value [cDNA *]

makes a new cDNA from a Sequence. It owns the Sequence memory, ie will attempt a /free_Sequence on the structure when /free_cDNA is called

If you want to give this cDNA this Sequence and forget about it, then just hand it this sequence and set seq to NULL (no need to free it). If you intend to use the sequence object elsewhere outside of the cDNA datastructure then use cDNA_from_Sequence(/hard_link_Sequence(seq))

27 cdnadb

This module contains the following objects

27.2 cDNADB
This module also contains some factory methods

27.1 cdnadb factory methods

27.1.1 new_cDNADB_from_single_seq

External C: Wise2_new_cDNADB_from_single_seq (seq)
Perl: &Wise2::new_cDNADB_from_single_seq (seq)

Arguments

seq: [UNKN ] sequence which as placed into cDNADB structure. [cDNA *]
returns: [UNKN ] Undocumented return value [cDNADB *]

To make a new cDNA database from a single cDNA Sequence with a eval system

27.1.2 new_cDNADB

External C: Wise2_new_cDNADB (seqdb)
Perl: &Wise2::new_cDNADB (seqdb)

Arguments

seqdb: [READ ] sequence database [SequenceDB *]
returns: [UNKN ] Undocumented return value [cDNADB *]

To make a new cDNA database

27.2 Object cDNADB

The cDNADB object has the following fields. To see how to access them refer to 7.1

: is_single_seq Type [boolean : Scalar] No documentation
: done_forward Type [boolean : Scalar] No documentation
: forward_only Type [boolean : Scalar] No documentation
: forw Type [ComplexSequence * : Scalar] No documentation
: rev Type [ComplexSequence * : Scalar] No documentation
: sdb Type [SequenceDB * : Scalar] No documentation
: current Type [Sequence * : Scalar] No documentation
: cses Type [ComplexSequenceEvalSet * : Scalar] No documentation
: error_handling Type [CdnaDBErrorType : Scalar] No documentation
: error_tol Type [double : Scalar] No documentation

This object hold a database of cDNA sequences.

You will probably use it in one of two ways

1 A sequence formatted database, which is provided by a /SequenceDB object is used to provide the raw sequences

2 A single cDNA sequence is used.

In each case this database provides both the forward and reverse strands into the system.

Notice that what is exported are /ComplexSequence objects, not cDNA dna, as this is what is generally needed. These are things with splice sites calculated etc. This is why for initialisation this needs a /ComplexSequenceEvalSet of the correct type.

Member functions of cDNADB

27.2.1 get_cDNA_from_cDNADB

External C: Wise2_get_cDNA_from_cDNADB (cdnadb,de)
Perl: &Wise2::cDNADB::get_entry (cdnadb,de)
Perl-OOP call: $obj->get_entry(de)

Arguments

cdnadb: [READ ] cDNA database [cDNADB *]
de: [READ ] DataEntry information [DataEntry *]
returns: [UNKN ] Undocumented return value [cDNA *]

Gets cDNA sequence out from the cDNADB using the information stored in dataentry

28 codon

This module contains the following objects

28.2 CodonTable
This module also contains some factory methods

28.1 codon factory methods

28.1.1 is_non_ambiguous_codon_seq

External C: Wise2_is_non_ambiguous_codon_seq (seq)
Perl: &Wise2::is_non_ambiguous_codon_seq (seq)

Arguments

seq: [READ ] pointer to DNA sequence [char *]
returns: [UNKN ] TRUE if real codon, FALSE if contains Nś [boolean]

Tells you if this codon is a real codon

28.1.2 codon_from_base4_codon

External C: Wise2_codon_from_base4_codon (c)
Perl: &Wise2::codon_from_base4_codon (c)

Arguments

c: [UNKN ] Undocumented argument [int]
returns: [UNKN ] Undocumented return value [codon]

maps a 0-63 codon to a 0-123 codon. Suprisingly useful.

28.1.3 base4_codon_from_codon

External C: Wise2_base4_codon_from_codon (c)
Perl: &Wise2::base4_codon_from_codon (c)

Arguments

c: [READ ] codon 0-125 [codon]
returns: [UNKN ] base 4 codon (0-63) [int]

maps a 0-125 codon to a 0-63 codon.

If ambiguous then returns 64 having issued a warning.

28.1.4 has_random_bases

External C: Wise2_has_random_bases (c)
Perl: &Wise2::has_random_bases (c)

Arguments

c: [READ ] codon number 0-124 [codon]
returns: [UNKN ] TRUE if has Nś , FALSE otherwise [boolean]

Tests to see if this codon number has any Nś in it

28.1.5 permute_possible_random_bases

External C: Wise2_permute_possible_random_bases (c,one,two,three)
Perl: &Wise2::permute_possible_random_bases (c,one,two,three)

Arguments

c: [READ ] codon number [codon]
one: [READ ] base to replace first position if N [base]
two: [READ ] base to replace second position if N [base]
three: [READ ] base to replace third position if N [base]
returns: [UNKN ] codon number [codon]

Bizarely useful function for calculating ambiguity scores.

This takes the codon c, and for each possible base, if it is N, replaces it with one, two or three.

If the base is not N, it remains the same

28.1.6 base_from_codon

External C: Wise2_base_from_codon (c,pos)
Perl: &Wise2::base_from_codon (c,pos)

Arguments

c: [UNKN ] Undocumented argument [codon]
pos: [UNKN ] Undocumented argument [int]
returns: [UNKN ] Undocumented return value [base]

Probably not the best function to use for this, but useful. Takes a codon and with pos being 1,2,3 gives you the firt,second of third base

28.1.7 codon_from_seq

External C: Wise2_codon_from_seq (seq)
Perl: &Wise2::codon_from_seq (seq)

Arguments

seq: [UNKN ] pointer to sequence of at least 3 chrs long. [char *]
returns: [UNKN ] Undocumented return value [codon]

takes an ASCII coded pointer to a 3 base pair sequence (it could be the part of a sequence: it only assummes that the seq points with 3 chars at pos 0,1,2 in C coordinates from seq. No NULL is required). It ives back the codon as made from standard mapping, ie, 25*base_1+5*base_2 + base3 being a number from 0-124 inc.

28.1.8 base4_codon_from_seq

External C: Wise2_base4_codon_from_seq (seq)
Perl: &Wise2::base4_codon_from_seq (seq)

Arguments

seq: [UNKN ] pointer to sequence of at least 3 chrs long [char *]
returns: [UNKN ] Undocumented return value [int]

Sometimes it is more useful to work in base64, ie, non N. this functions does the same thing as /codon_from_seq but produces a seq being 16*base1 + 4 *base2 + base3

28.1.9 char_from_base

External C: Wise2_char_from_base (b)
Perl: &Wise2::char_from_base (b)

Arguments

b: [UNKN ] Undocumented argument [base]
returns: [UNKN ] Undocumented return value [char]

maps a base number (-04 inc) to A,T,G,C,N

28.1.10 base_from_char

External C: Wise2_base_from_char (c)
Perl: &Wise2::base_from_char (c)

Arguments

c: [UNKN ] Undocumented argument [char]
returns: [UNKN ] Undocumented return value [base]

maps a char (atcgn) to number, case insensitive, returns BASE_N if not atcgn

28.1.11 char_complement_base

External C: Wise2_char_complement_base (c)
Perl: &Wise2::char_complement_base (c)

Arguments

c: [UNKN ] Undocumented argument [char]
returns: [UNKN ] Undocumented return value [char]

the char equivalent of /complement_base. this gives the complement in char of a base in char. Does not check for non ATGCN

28.1.12 complement_base

External C: Wise2_complement_base (b)
Perl: &Wise2::complement_base (b)

Arguments

b: [UNKN ] Undocumented argument [base]
returns: [UNKN ] Undocumented return value [base]

gives back the complement as a number ofthe base (given as a number)

28.2 Object CodonTable

The CodonTable object has the following fields. To see how to access them refer to 7.1

: codon_str[125] Type [aa : Scalar] No documentation
: name Type [char * : Scalar] No documentation

The codon table provides a mapping from the 64 codons to the 20 amino acids. The rest of the modules provides assorted codon<->base<->amino acid mappings.

Probably the trickiest thing is that there are two different types of representations of codons. One in base 5 (N being the 5th base), providing 0-124 inclusive codon numbers. These numbers are the ones going to be principly used in most calculations.

However, it is often very useful to use 0-63 numbers, for example in the precise definition of the codon table.

Member functions of CodonTable

28.2.1 read_CodonTable_file

External C: Wise2_read_CodonTable_file (file)
Perl: &Wise2::CodonTable::read_CodonTable_file (file)
Perl-OOP call: $obj->read_CodonTable_file()

Arguments

file: [READ ] filename to open [char *]
returns: [OWNER] A codon-table, NULL if error [CodonTable *]

Opens filename, reads it as if a Ewan style codon table and closes.

28.2.2 read_CodonTable

External C: Wise2_read_CodonTable (ifp)
Perl: &Wise2::CodonTable::read_CodonTable (ifp)
Perl-OOP call: $obj->read_CodonTable()

Arguments

ifp: [READ ] file input [FILE *]
returns: [UNKN ] Undocumented return value [CodonTable *]

reads a codon table from a filestream in Ewan format.

As Ewan format is really bad and has no start/stop this will effectively read to the end of the file. Ooops.

28.2.3 aminoacid_from_seq

External C: Wise2_aminoacid_from_seq (ct,seq)
Perl: &Wise2::CodonTable::aminoacid_from_seq (ct,seq)
Perl-OOP call: $obj->aminoacid_from_seq(seq)

Arguments

ct: [READ ] codon table [CodonTable *]
seq: [READ ] pointer to DNA chars [char *]
returns: [UNKN ] an amino acid char (A-Z) [aa]

Returns the amino acid for this position in the DNA sequence Takes the pointer +1 and +2 points.

No error checks implemented. Probably a mistake ;)

28.2.4 aminoacid_from_codon

External C: Wise2_aminoacid_from_codon (ct,c)
Perl: &Wise2::CodonTable::aminoacid_from_codon (ct,c)
Perl-OOP call: $obj->aminoacid_from_codon(c)

Arguments

ct: [READ ] codon table [CodonTable *]
c: [READ ] codon number [codon]
returns: [READ ] aminoacid that is this codon (X for ambiguous, * for stop) [aa]

returns amino acid for this codon number (NB codon numbers 0-125)

28.2.5 is_stop_codon

External C: Wise2_is_stop_codon (c,ct)
Perl: &Wise2::CodonTable::is_stop_codon (c,ct)
Perl-OOP call: $obj->is_stop_codon(ct)

Arguments

c: [READ ] codon number [codon]
ct: [READ ] codon table [CodonTable *]
returns: [UNKN ] TRUE if is stop, FALSE otherwise [boolean]

tells you whether this codon number is really a stop in this translation table

28.2.6 is_valid_aminoacid

External C: Wise2_is_valid_aminoacid (ct,c)
Perl: &Wise2::CodonTable::is_valid_aminoacid (ct,c)
Perl-OOP call: $obj->is_valid_aminoacid(c)

Arguments

ct: [READ ] Codon Table [CodonTable *]
c: [UNKN ] aminoacid [char]
returns: [UNKN ] TRUE if valid, FALSE if not. [boolean]

Tells you if this letter (c) is recognised as a valid amino acid in this codon table

29 codonmapper

This module contains the following objects

29.2 CodonMapper
This module also contains some factory methods

29.1 codonmapper factory methods

29.1.1 flat_CodonMapper

External C: Wise2_flat_CodonMapper (ct)
Perl: &Wise2::flat_CodonMapper (ct)

Arguments

ct: [UNKN ] Codon Table giving codon->aa info [CodonTable *]
returns: [UNKN ] Undocumented return value [CodonMapper *]

Makes a CodonMapper with no codon bias or error possiblities from codon table

29.2 Object CodonMapper

The CodonMapper object has the following fields. To see how to access them refer to 7.1

: ct Type [CodonTable * : Scalar] hard-linked!
: codon_map[125][26] Type [double : Scalar] No documentation

CodonMapper holds a matrix of 125 by 26 to provide a mapping between a probabilities calculated on amino acids to triplet codon probabilities. This mapping takes into account 3 things

  1) The CodonTable
  2) The distribution of synonmous codons (codon bias)
  3) substitution errors

Member functions of CodonMapper

29.2.1 sprinkle_errors_over_CodonMapper

External C: Wise2_sprinkle_errors_over_CodonMapper (cm,error)
Perl: &Wise2::CodonMapper::sprinkle_errors_over_CodonMapper (cm,error)
Perl-OOP call: $obj->sprinkle_errors_over_CodonMapper(error)

Arguments

cm: [READ ] CodonMapper to be sprinkled [CodonMapper *]
error: [UNKN ] substitution error rate [double]
returns: Nothing - no return value

Takes a codon mapper and assummes that the majority of errors are due to a single base change in the codon at probability error. Therefore, for each codon it adds error * prob(codon) * 0.25 to each other codon one base away, taking away therefore the result.

30 complexsequence

This module contains the following objects

30.1 ComplexSequence
30.2 ComplexSequenceEvalSet

30.1 Object ComplexSequence

The ComplexSequence object has the following fields. To see how to access them refer to 7.1

: type Type [int : Scalar] No documentation
: seq Type [Sequence * : Scalar] No documentation
: data Type [int * : Scalar] No documentation
: datastore Type [int * : Scalar] No documentation
: depth Type [int : Scalar] No documentation
: length Type [int : Scalar] No documentation
: creator Type [ComplexSequenceEvalSet * : Scalar] what made it

A ComplexSequence is an abstraction of a Sequence which can be handily made using ComplexSequenceEval functions and is efficiently laid out in memory.

Member functions of ComplexSequence

30.2 Object ComplexSequenceEvalSet

The ComplexSequenceEvalSet object has the following fields. To see how to access them refer to 7.1

: type Type [int : Scalar] No documentation
: has_been_prepared Type [boolean : Scalar] No documentation
: left_window Type [int : Scalar] overall sequence eval
: right_window Type [int : Scalar] overall sequence eval
: left_lookback Type [int : Scalar] overall sequence eval
: cse Type [ComplexSequenceEval ** : List] No documentation

This object holds a collection of ComplexSequenceEvals. Its role is to define the sequence specific parts of a dynamic programming algorithm as computable functions.

Ideally you should use pre-made ComplexSequenceEvalSets as it will save you alot of grief

Member functions of ComplexSequenceEvalSet

31 compmat

This module contains the following objects

31.1 CompMat

31.1 Object CompMat

The CompMat object has the following fields. To see how to access them refer to 7.1

: comp[26][26] Type [Score : Scalar] No documentation
: name Type [char * : Scalar] if any, could be NULL

This object stores BLOSUM and PAM comparison matrices. It stores them as scores: NB - this means probabilistically we are talking about some arbitary base of log which is really annoying.

Member functions of CompMat

31.1.1 fail_safe_CompMat_access

External C: Wise2_fail_safe_CompMat_access (cm,aa1,aa2)
Perl: &Wise2::CompMat::fail_safe_CompMat_access (cm,aa1,aa2)
Perl-OOP call: $obj->fail_safe_CompMat_access(aa1,aa2)

Arguments

cm: [UNKN ] compmat object [CompMat *]
aa1: [UNKN ] first amino acid [int]
aa2: [UNKN ] second amino acid [int]
returns: [UNKN ] Undocumented return value [Score]

gives the fail form of the macro CompMat_AAMATCH which checks that aa1 and a2 are sensible and that cm is not NULL.

31.1.2 write_Blast_CompMat

External C: Wise2_write_Blast_CompMat (cm,ofp)
Perl: &Wise2::CompMat::write_Blast_CompMat (cm,ofp)
Perl-OOP call: $obj->write_Blast_CompMat(ofp)

Arguments

cm: [UNKN ] CompMat object [CompMat *]
ofp: [UNKN ] file to output [FILE *]
returns: [UNKN ] Undocumented return value [boolean]

writes a protien CompMat with a standard alphabet.

31.1.3 read_Blast_file_CompMat

External C: Wise2_read_Blast_file_CompMat (filename)
Perl: &Wise2::CompMat::read_Blast_file_CompMat (filename)
Perl-OOP call: $obj->read_Blast_file_CompMat()

Arguments

filename: [UNKN ] Undocumented argument [char *]
returns: [UNKN ] Undocumented return value [CompMat *]

Opens file, reads matrix, closes file. calls /read_Blast_CompMat for the actual format reading. Uses /openfile to open the file, so will open from config files.

31.1.4 read_Blast_CompMat

External C: Wise2_read_Blast_CompMat (ifp)
Perl: &Wise2::CompMat::read_Blast_CompMat (ifp)
Perl-OOP call: $obj->read_Blast_CompMat()

Arguments

ifp: [UNKN ] Undocumented argument [FILE *]
returns: [UNKN ] Undocumented return value [CompMat *]

reads a BLAST format matrix and allocates a new ComMat structure.

32 dbsearchimpl

This module contains the following objects

32.2 DBSearchImpl
This module also contains some factory methods

32.1 dbsearchimpl factory methods

32.1.1 new_pthread_DBSearchImpl

External C: Wise2_new_pthread_DBSearchImpl (void)
Perl: &Wise2::new_pthread_DBSearchImpl ()

Arguments

returns: [UNKN ] Undocumented return value [DBSearchImpl *]

Makes a new pthreaded DBSearchImpl

For use mainly for apiś who dont want to initalize the object from the command line

32.1.2 new_serial_DBSearchImpl

External C: Wise2_new_serial_DBSearchImpl (void)
Perl: &Wise2::new_serial_DBSearchImpl ()

Arguments

returns: [UNKN ] Undocumented return value [DBSearchImpl *]

Makes a new serial DBSearchImpl

For use mainly for apiś who dont want to initalize the object from the command line

32.2 Object DBSearchImpl

The DBSearchImpl object has the following fields. To see how to access them refer to 7.1

: type Type [int : Scalar] No documentation
: trace_level Type [int : Scalar] how much debugging information to print
: trace_file Type [FILE * : Scalar] for writing out trace of db stuff
: suggest_thread_no Type [int : Scalar] default, -1, means the use a call to _SC_NPROC
: search_routine Type [int : Scalar] routine used for the calculation, exact/kbest

DBSearchImpl contains the information about the database search implementation used in a dynamite call. This is the only object which needs to specify say threads vs serial code

The construction of this object is from its own stripping of the command line. This way programs which do database searching need not know anything about the implementation that is being used by the dynamite compiler

From the API you should be looking to make it from one of the handy constructors.

For the threads code, if you leave the suggest_thread_no at (-1) - what it usually comes as for the constructors. The system then figures out the number of processors available on the machine and sets it to that.

Member functions of DBSearchImpl

32.2.1 impl_string_DBSearchImpl

External C: Wise2_impl_string_DBSearchImpl (dbsi)
Perl: &Wise2::DBSearchImpl::string (dbsi)
Perl-OOP call: $obj->string()

Arguments

dbsi: [UNKN ] Undocumented argument [DBSearchImpl *]
returns: [SOFT ] string of the search implementation [char *]

Gets a static text string out of the search implementation

33 dnamatrix

This module contains the following objects

33.2 DnaMatrix
33.3 DnaProbMatrix
This module also contains some factory methods

33.1 dnamatrix factory methods

33.1.1 identity_DnaMatrix

External C: Wise2_identity_DnaMatrix (id_score,mismatch)
Perl: &Wise2::identity_DnaMatrix (id_score,mismatch)

Arguments

id_score: [UNKN ] score of idenities [Score]
mismatch: [UNKN ] score of mistmatches [Score]
returns: [UNKN ] Undocumented return value [DnaMatrix *]

makes an idenity matrix wth id_score on the leading diagonal and mismatch elsewhere.

33.1.2 DnaProbMatrix_from_match

External C: Wise2_DnaProbMatrix_from_match (match,nmask_type)
Perl: &Wise2::DnaProbMatrix_from_match (match,nmask_type)

Arguments

match: [UNKN ] Undocumented argument [Probability]
nmask_type: [UNKN ] Undocumented argument [int]
returns: [UNKN ] Undocumented return value [DnaProbMatrix *]

Makes a probability matrix from simple match/mismatch probabilities.

33.1.3 DnaMatrix_from_DnaProbMatrix

External C: Wise2_DnaMatrix_from_DnaProbMatrix (dpm)
Perl: &Wise2::DnaMatrix_from_DnaProbMatrix (dpm)

Arguments

dpm: [UNKN ] Undocumented argument [DnaProbMatrix *]
returns: [UNKN ] Undocumented return value [DnaMatrix *]

Maps probabilities to scores

33.2 Object DnaMatrix

The DnaMatrix object has the following fields. To see how to access them refer to 7.1

: score[5][5] Type [Score : Scalar] No documentation

No documentation for DnaMatrix

Member functions of DnaMatrix

33.3 Object DnaProbMatrix

The DnaProbMatrix object has the following fields. To see how to access them refer to 7.1

: prob[5][5] Type [Probability : Scalar] No documentation

No documentation for DnaProbMatrix

Member functions of DnaProbMatrix

33.3.1 flat_null_DnaProbMatrix

External C: Wise2_flat_null_DnaProbMatrix (dpm)
Perl: &Wise2::DnaProbMatrix::flat_null_DnaProbMatrix (dpm)
Perl-OOP call: $obj->flat_null_DnaProbMatrix()

Arguments

dpm: [UNKN ] Undocumented argument [DnaProbMatrix *]
returns: Nothing - no return value

makes a odds of dpm via a 0.25 factor into each base.

34 gene

This module contains the following objects

34.1 Gene

34.1 Object Gene

The Gene object has the following fields. To see how to access them refer to 7.1

: start Type [int : Scalar] No documentation
: end Type [int : Scalar] No documentation
: parent Type [GenomicRegion * : Scalar] may not be here
: genomic Type [Genomic * : Scalar] may not be here!
: transcript Type [Transcript ** : List] No documentation
: name Type [char * : Scalar] ugly . Need a better system
: bits Type [double : Scalar] ditto...
: seqname Type [char * : Scalar] very bad! this is for keeping track of what sequence was used to make the gene
: ispseudo Type [boolean : Scalar] is a pseudogene or not

Gene is the datastructure which represents a single gene. At the moment this is considered to be a series of transcripts (the first transcript being "canonical") which are made from a certain start/stop region in genomic DNA.

The gene datastructure does not necessarily contain any DNA sequence. When someone askes for the gene sequence, via get_Genomic_from_Gene(), it first sees if there is anything in the Genomic * ćache. If this is null, it then looks at parent (the genomic region), and if that is null, complains and returns null. Otherwise it truncates its parentś dna in the correct way, places the data structure into the genomic * cache, and returns it.

The name, bits and seqname have put into this datastructure for convience of carrying around this information into some of the gene prediction output formats. Probabaly

  o they should be in transcript, not gene
  o they shouldn\'t be here at all.

<sigh>

Member functions of Gene

34.1.1 get_Genomic_from_Gene

External C: Wise2_get_Genomic_from_Gene (gene)
Perl: &Wise2::Gene::get_Genomic_from_Gene (gene)
Perl-OOP call: $obj->get_Genomic_from_Gene()

Arguments

gene: [READ ] gene to get Genomic from [Gene *]
returns: [SOFT ] Genomic DNA data structure [Genomic *]

Gives back a Genomic sequence type from a gene.

34.1.2 show_pretty_Gene

External C: Wise2_show_pretty_Gene (ge,show_supporting,ofp)
Perl: &Wise2::Gene::show_pretty_Gene (ge,show_supporting,ofp)
Perl-OOP call: $obj->show_pretty_Gene(show_supporting,ofp)

Arguments

ge: [UNKN ] Undocumented argument [Gene *]
show_supporting: [UNKN ] Undocumented argument [boolean]
ofp: [UNKN ] Undocumented argument [FILE *]
returns: Nothing - no return value

Shows a gene in the biologically accepted form

35 genomic

This module contains the following objects

35.2 Genomic
35.3 GenomicRepeat
This module also contains some factory methods

35.1 genomic factory methods

35.1.1 read_fasta_file_Genomic

External C: Wise2_read_fasta_file_Genomic (filename,length_of_N)
Perl: &Wise2::read_fasta_file_Genomic (filename,length_of_N)

Arguments

filename: [UNKN ] filename to be opened and read [char *]
length_of_N: [UNKN ] length of N to be considered repeat. -1 means none [int]
returns: [UNKN ] Undocumented return value [Genomic *]

Reads a fasta file assumming that it is Genomic. Will complain if it is not, and return NULL.

35.1.2 read_fasta_Genomic

External C: Wise2_read_fasta_Genomic (ifp,length_of_N)
Perl: &Wise2::read_fasta_Genomic (ifp,length_of_N)

Arguments

ifp: [UNKN ] file point to be read from [FILE *]
length_of_N: [UNKN ] length of N to be considered repeat. -1 means none [int]
returns: [UNKN ] Undocumented return value [Genomic *]

Reads a fasta file assumming that it is Genomic. Will complain if it is not, and return NULL.

35.1.3 Genomic_from_Sequence_Nheuristic

External C: Wise2_Genomic_from_Sequence_Nheuristic (seq,length_of_N)
Perl: &Wise2::Genomic_from_Sequence_Nheuristic (seq,length_of_N)

Arguments

seq: [UNKN ] Undocumented argument [Sequence *]
length_of_N: [UNKN ] Undocumented argument [int]
returns: [UNKN ] Undocumented return value [Genomic *]

makes a new genomic from a Sequence, but assummes that all the N runs greater than a certain level are actually repeats.

35.1.4 Genomic_from_Sequence

External C: Wise2_Genomic_from_Sequence (seq)
Perl: &Wise2::Genomic_from_Sequence (seq)

Arguments

seq: [OWNER] Sequence to make genomic from [Sequence *]
returns: [UNKN ] Undocumented return value [Genomic *]

makes a new genomic from a Sequence. It owns the Sequence memory, ie will attempt a /free_Sequence on the structure when /free_Genomic is called

If you want to give this genomic this Sequence and forget about it, then just hand it this sequence and set seq to NULL (no need to free it). If you intend to use the sequence object elsewhere outside of the Genomic datastructure then use Genomic_from_Sequence(/hard_link_Sequence(seq))

This is part of a strict typing system, and therefore is going to convert all non ATGCNs to Ns. You will lose information here.

35.2 Object Genomic

The Genomic object has the following fields. To see how to access them refer to 7.1

: baseseq Type [Sequence * : Scalar] No documentation
: repeat Type [GenomicRepeat ** : List] No documentation

No documentation for Genomic

Member functions of Genomic

35.2.1 truncate_Genomic

External C: Wise2_truncate_Genomic (gen,start,stop)
Perl: &Wise2::Genomic::truncate_Genomic (gen,start,stop)
Perl-OOP call: $obj->truncate_Genomic(start,stop)

Arguments

gen: [READ ] Genomic that is truncated [Genomic *]
start: [UNKN ] Undocumented argument [int]
stop: [UNKN ] Undocumented argument [int]
returns: [UNKN ] Undocumented return value [Genomic *]

Truncates a Genomic sequence. Basically uses the /magic_trunc_Sequence function (of course!)

It does not alter gen, rather it returns a new sequence with that truncation

Handles repeat information correctly.

35.2.2 Genomic_name

External C: Wise2_Genomic_name (gen)
Perl: &Wise2::Genomic::Genomic_name (gen)
Perl-OOP call: $obj->Genomic_name()

Arguments

gen: [UNKN ] Undocumented argument [Genomic *]
returns: [UNKN ] Undocumented return value [char *]

Returns the name of the Genomic

35.2.3 Genomic_length

External C: Wise2_Genomic_length (gen)
Perl: &Wise2::Genomic::Genomic_length (gen)
Perl-OOP call: $obj->Genomic_length()

Arguments

gen: [UNKN ] Undocumented argument [Genomic *]
returns: [UNKN ] Undocumented return value [int]

Returns the length of the Genomic

35.2.4 Genomic_seqchar

External C: Wise2_Genomic_seqchar (gen,pos)
Perl: &Wise2::Genomic::Genomic_seqchar (gen,pos)
Perl-OOP call: $obj->Genomic_seqchar(pos)

Arguments

gen: [UNKN ] Genomic [Genomic *]
pos: [UNKN ] position in Genomic to get char [int]
returns: [UNKN ] Undocumented return value [char]

Returns sequence character at this position.

35.2.5 show_Genomic

External C: Wise2_show_Genomic (gen,ofp)
Perl: &Wise2::Genomic::show (gen,ofp)
Perl-OOP call: $obj->show(ofp)

Arguments

gen: [UNKN ] Undocumented argument [Genomic *]
ofp: [UNKN ] Undocumented argument [FILE *]
returns: Nothing - no return value

For debugging

35.3 Object GenomicRepeat

The GenomicRepeat object has the following fields. To see how to access them refer to 7.1

: start Type [int : Scalar] No documentation
: end Type [int : Scalar] No documentation
: type Type [char * : Scalar] No documentation

No documentation for GenomicRepeat

Member functions of GenomicRepeat

36 genomicdb

This module contains the following objects

36.2 GenomicDB
This module also contains some factory methods

36.1 genomicdb factory methods

36.1.1 new_GenomicDB_from_single_seq

External C: Wise2_new_GenomicDB_from_single_seq (gen,cses,score_in_repeat_coding)
Perl: &Wise2::new_GenomicDB_from_single_seq (gen,cses,score_in_repeat_coding)

Arguments

gen: [UNKN ] sequence which as placed into GenomicDB structure. [Genomic *]
cses: [UNKN ] Undocumented argument [ComplexSequenceEvalSet *]
score_in_repeat_coding: [UNKN ] Undocumented argument [int]
returns: [UNKN ] Undocumented return value [GenomicDB *]

To make a new genomic database from a single Genomic Sequence with a eval system

36.1.2 new_GenomicDB

External C: Wise2_new_GenomicDB (seqdb,cses,length_of_N,repeat_in_cds_score)
Perl: &Wise2::new_GenomicDB (seqdb,cses,length_of_N,repeat_in_cds_score)

Arguments

seqdb: [UNKN ] sequence database [SequenceDB *]
cses: [UNKN ] protein evaluation set [ComplexSequenceEvalSet *]
length_of_N: [UNKN ] Undocumented argument [int]
repeat_in_cds_score: [UNKN ] Undocumented argument [int]
returns: [UNKN ] Undocumented return value [GenomicDB *]

To make a new genomic database

36.2 Object GenomicDB

The GenomicDB object has the following fields. To see how to access them refer to 7.1

: is_single_seq Type [boolean : Scalar] No documentation
: done_forward Type [boolean : Scalar] No documentation
: forw Type [ComplexSequence * : Scalar] No documentation
: rev Type [ComplexSequence * : Scalar] No documentation
: sdb Type [SequenceDB * : Scalar] No documentation
: current Type [Genomic * : Scalar] No documentation
: cses Type [ComplexSequenceEvalSet * : Scalar] No documentation
: error_handling Type [GenDBErrorType : Scalar] No documentation
: single Type [Genomic * : Scalar] for single sequence cases, so we can índexón it
: revsingle Type [Genomic * : Scalar] No documentation
: length_of_N Type [int : Scalar] No documentation
: repeat_in_cds_score Type [int : Scalar] No documentation

This object hold a database of genomic sequences.

You will probably use it in one of two ways

1 A sequence formatted database, which is provided by a /SequenceDB object is used to provide the raw sequences

2 A single Genomic sequence is used.

In each case this database provides both the forward and reverse strands into the system.

Notice that what is exported are /ComplexSequence objects, not genomic dna, as this is what is generally needed. These are things with splice sites calculated etc. This is why for initialisation this needs a /ComplexSequenceEvalSet of the correct type.

Member functions of GenomicDB

36.2.1 get_Genomic_from_GenomicDB

External C: Wise2_get_Genomic_from_GenomicDB (gendb,de)
Perl: &Wise2::GenomicDB::get_Genomic_from_GenomicDB (gendb,de)
Perl-OOP call: $obj->get_Genomic_from_GenomicDB(de)

Arguments

gendb: [UNKN ] Undocumented argument [GenomicDB *]
de: [UNKN ] Undocumented argument [DataEntry *]
returns: [UNKN ] Undocumented return value [Genomic *]

Gets Genomic sequence out from the GenomicDB using the information stored in dataentry

37 genomicregion

This module contains the following objects

37.1 GenomicRegion

37.1 Object GenomicRegion

The GenomicRegion object has the following fields. To see how to access them refer to 7.1

: gene Type [Gene ** : List] No documentation
: genomic Type [Genomic * : Scalar] NB, may not be here. Be careful!

GenomicRegion is structure which represents information on a region of genomic DNA. It *may not* have the actual DNA sequence in there, and it is important to realise that.

The numbering scheme of many other things (eg, genes) are going to be represented in this guys coordinates

Member functions of GenomicRegion

37.1.1 new_GenomicRegion

External C: Wise2_new_GenomicRegion (gen)
Perl: &Wise2::GenomicRegion::new_GenomicRegion (gen)
Perl-OOP call: $obj->new_GenomicRegion()

Arguments

gen: [UNKN ] Undocumented argument [Genomic *]
returns: [UNKN ] Undocumented return value [GenomicRegion *]

makes a genomicregion from a genomic sequence

37.1.2 read_EMBL_GenomicRegion_file

External C: Wise2_read_EMBL_GenomicRegion_file (filename)
Perl: &Wise2::GenomicRegion::read_EMBL_GenomicRegion_file (filename)
Perl-OOP call: $obj->read_EMBL_GenomicRegion_file()

Arguments

filename: [UNKN ] Undocumented argument [char *]
returns: [UNKN ] Undocumented return value [GenomicRegion *]

Reads in both EMBL sequence and features

37.1.3 read_EMBL_GenomicRegion

External C: Wise2_read_EMBL_GenomicRegion (ifp)
Perl: &Wise2::GenomicRegion::read_EMBL_GenomicRegion (ifp)
Perl-OOP call: $obj->read_EMBL_GenomicRegion()

Arguments

ifp: [UNKN ] Undocumented argument [FILE *]
returns: [UNKN ] Undocumented return value [GenomicRegion *]

Reads in both EMBL sequence and features

37.1.4 add_Gene_to_GenomicRegion

External C: Wise2_add_Gene_to_GenomicRegion (gr,gene)
Perl: &Wise2::GenomicRegion::add_Gene_to_GenomicRegion (gr,gene)
Perl-OOP call: $obj->add_Gene_to_GenomicRegion(gene)

Arguments

gr: [UNKN ] GenomicRegion to be added to [GenomicRegion *]
gene: [UNKN ] Gene to be added [Gene *]
returns: [UNKN ] Undocumented return value [boolean]

adds a Gene to this GenomicRegion, making sure that it parent/son relationship is ok

37.1.5 show_ace_GenomicRegion

External C: Wise2_show_ace_GenomicRegion (gr,seq_name,ofp)
Perl: &Wise2::GenomicRegion::show_ace_GenomicRegion (gr,seq_name,ofp)
Perl-OOP call: $obj->show_ace_GenomicRegion(seq_name,ofp)

Arguments

gr: [UNKN ] Undocumented argument [GenomicRegion *]
seq_name: [UNKN ] Undocumented argument [char *]
ofp: [UNKN ] Undocumented argument [FILE *]
returns: Nothing - no return value

shows ACeDB subsequence source.

Assummes

  a only one transcript per gene
  b only cds exons are used

37.1.6 show_pretty_GenomicRegion

External C: Wise2_show_pretty_GenomicRegion (gr,show_supporting,ofp)
Perl: &Wise2::GenomicRegion::show_pretty_GenomicRegion (gr,show_supporting,ofp)
Perl-OOP call: $obj->show_pretty_GenomicRegion(show_supporting,ofp)

Arguments

gr: [UNKN ] Undocumented argument [GenomicRegion *]
show_supporting: [UNKN ] Undocumented argument [boolean]
ofp: [UNKN ] Undocumented argument [FILE *]
returns: Nothing - no return value

shows ṕrettybio type gene

37.1.7 write_Diana_FT_GenomicRegion

External C: Wise2_write_Diana_FT_GenomicRegion (gr,ofp)
Perl: &Wise2::GenomicRegion::write_Diana_FT_GenomicRegion (gr,ofp)
Perl-OOP call: $obj->write_Diana_FT_GenomicRegion(ofp)

Arguments

gr: [UNKN ] Undocumented argument [GenomicRegion *]
ofp: [UNKN ] Undocumented argument [FILE *]
returns: Nothing - no return value

Writes Embl feature table for diana use. Does assumme that there is only one transcript per gene and only cds exons are used

Output like

   FT   misc\_feature       join(100..200)

37.1.8 write_Embl_FT_GenomicRegion

External C: Wise2_write_Embl_FT_GenomicRegion (gr,ofp)
Perl: &Wise2::GenomicRegion::write_Embl_FT_GenomicRegion (gr,ofp)
Perl-OOP call: $obj->write_Embl_FT_GenomicRegion(ofp)

Arguments

gr: [UNKN ] Undocumented argument [GenomicRegion *]
ofp: [UNKN ] Undocumented argument [FILE *]
returns: Nothing - no return value

Writes Embl feature table. Does assumme that there is only one transcript per gene and only cds exons are used

Output like

   FT   CDS          join(100..200)

38 histogram

This module contains the following objects

38.2 Histogram
This module also contains some factory methods

38.1 histogram factory methods

38.1.1 new_Histogram

External C: Wise2_new_Histogram (min,max,lumpsize)
Perl: &Wise2::new_Histogram (min,max,lumpsize)

Arguments

min: [UNKN ] minimum score (integer) [int]
max: [UNKN ] maximum score (integer) [int]
lumpsize: [UNKN ] when reallocating histogram, the reallocation amount [int]
returns: [UNKN ] Undocumented return value [Histogram *]

This function was written by Sean Eddy as part of his HMMer2 histogram.c module

Converted by Ewan Birney to Dynamite source June 98. Copyright is LGPL. For more info read READMEs

Documentation:

Allocate and return a histogram structure. min and max are your best guess. They need not be absolutely correct; the histogram will expand dynamically to accomodate scores that exceed these suggested bounds. The amount that the histogram grows by is set by "lumpsize"

Was called AllocHistorgram new_Historgram is more wise2-ish

38.2 Object Histogram

The Histogram object has the following fields. To see how to access them refer to 7.1

: histogram Type [int* : Scalar] counts of hits
: min Type [int : Scalar] elem 0 of histogram == min
: max Type [int : Scalar] last elem of histogram == max
: highscore Type [int : Scalar] highest active elem has this score
: lowscore Type [int : Scalar] lowest active elem has this score
: lumpsize Type [int : Scalar] when resizing, overalloc by this
: total Type [int : Scalar] total # of hits counted
: expect Type [float* : Scalar] expected counts of hits
: fit_type Type [int : Scalar] flag indicating distribution type
: param[3] Type [float : Scalar] parameters used for fits
: chisq Type [float : Scalar] chi-squared val for goodness of fit
: chip Type [float : Scalar] P value for chisquared

This Object came from Sean Eddy excellent histogram package. He donated it free of all restrictions to allow it to be used in the Wise2 package without complicated licensing terms. He is a *very* nice man.

It was made into a dynamite object so that

   a) External ports to scripting languages would be trivial
   b) cooperation with future versions of histogram.c would be possible.

Here is the rest of the documentation from sean.

Keep a score histogram.

The main implementation issue here is that the range of scores is unknown, and will go negative. histogram is a 0..max-min array that represents the range min..max. A given score is indexed in histogram array as score-min. The AddToHistogram function deals with dynamically resizing the histogram array when necessary.

Member functions of Histogram

38.2.1 UnfitHistogram

External C: Wise2_UnfitHistogram (h)
Perl: &Wise2::Histogram::UnfitHistogram (h)
Perl-OOP call: $obj->UnfitHistogram()

Arguments

h: [UNKN ] Undocumented argument [Histogram *]
returns: Nothing - no return value

This function was written by Sean Eddy as part of his HMMer2 histogram.c module

Converted by Ewan Birney to Dynamite source June 98. Copyright is LGPL. For more info read READMEs

Documentation:

Free only the theoretical fit part of a histogram.

38.2.2 AddToHistogram

External C: Wise2_AddToHistogram (h,sc)
Perl: &Wise2::Histogram::add (h,sc)
Perl-OOP call: $obj->add(sc)

Arguments

h: [UNKN ] Undocumented argument [Histogram *]
sc: [UNKN ] Undocumented argument [float]
returns: Nothing - no return value

This function was written by Sean Eddy as part of his HMMer2 histogram.c module

Converted by Ewan Birney to Dynamite source June 98. Copyright is LGPL. For more info read READMEs

Documentation:

Bump the appropriate counter in a histogram structure, given a score. The score is rounded off from float precision to the next lower integer.

38.2.3 PrintASCIIHistogram

External C: Wise2_PrintASCIIHistogram (h,fp)
Perl: &Wise2::Histogram::show (h,fp)
Perl-OOP call: $obj->show(fp)

Arguments

h: [UNKN ] histogram to print [Histogram *]
fp: [UNKN ] open file to print to (stdout works) [FILE *]
returns: Nothing - no return value

This function was written by Sean Eddy as part of his HMMer2 histogram.c module

Converted by Ewan Birney to Dynamite source June 98. Copyright is LGPL. For more info read READMEs

Documentation:

Print a "prettified" histogram to a file pointer. Deliberately a look-and-feel clone of Bill Pearsonś excellent FASTA output.

38.2.4 EVDBasicFit

External C: Wise2_EVDBasicFit (h)
Perl: &Wise2::Histogram::EVDBasicFit (h)
Perl-OOP call: $obj->EVDBasicFit()

Arguments

h: [UNKN ] histogram to fit [Histogram *]
returns: Nothing - no return value

This function was written by Sean Eddy as part of his HMMer2 histogram.c module

Converted by Ewan Birney to Dynamite source June 98. Copyright is LGPL. For more info read READMEs

Documentation:

Fit a score histogram to the extreme value distribution. Set the parameters lambda and mu in the histogram structure. Fill in the expected values in the histogram. Calculate a chi-square test as a measure of goodness of fit.

This is the basic version of ExtremeValueFitHistogram(), in a nonrobust form: simple linear regression with no outlier pruning.

Methods: Uses a linear regression fitting method [Collins88,Lawless82]

38.2.5 ExtremeValueFitHistogram

External C: Wise2_ExtremeValueFitHistogram (h,censor,high_hint)
Perl: &Wise2::Histogram::fit_EVD (h,censor,high_hint)
Perl-OOP call: $obj->fit_EVD(censor,high_hint)

Arguments

h: [UNKN ] histogram to fit [Histogram *]
censor: [UNKN ] TRUE to censor data left of the peak [int]
high_hint: [UNKN ] score cutoff; above this are real hits that arent fit [float]
returns: [UNKN ] if fit is judged to be valid else 0 if fit is invalid (too few seqs.) [int]

This function was written by Sean Eddy as part of his HMMer2 histogram.c module

Converted by Ewan Birney to Dynamite source June 98. Copyright is LGPL. For more info read READMEs

Documentation:

Purpose: Fit a score histogram to the extreme value distribution. Set the parameters lambda and mu in the histogram structure. Calculate a chi-square test as a measure of goodness of fit.

Methods: Uses a maximum likelihood method [Lawless82]. Lower outliers are removed by censoring the data below the peak. Upper outliers are removed iteratively using method described by [Mott92].

38.2.6 ExtremeValueSetHistogram

External C: Wise2_ExtremeValueSetHistogram (h,mu,lambda,lowbound,highbound,wonka,ndegrees)
Perl: &Wise2::Histogram::set_EVD (h,mu,lambda,lowbound,highbound,wonka,ndegrees)
Perl-OOP call: $obj->set_EVD(mu,lambda,lowbound,highbound,wonka,ndegrees)

Arguments

h: [UNKN ] the histogram to set [Histogram *]
mu: [UNKN ] mu location parameter [float]
lambda: [UNKN ] lambda scale parameter [float]
lowbound: [UNKN ] low bound of the histogram that was fit [float]
highbound: [UNKN ] high bound of histogram that was fit [float]
wonka: [UNKN ] fudge factor; fraction of hits estimated to be "EVD-like" [float]
ndegrees: [UNKN ] extra degrees of freedom to subtract in chi2 test: [int]
returns: Nothing - no return value

This function was written by Sean Eddy as part of his HMMer2 histogram.c module

Converted by Ewan Birney to Dynamite source June 98. Copyright is LGPL. For more info read READMEs

Documentation:

Instead of fitting the histogram to an EVD, simply set the EVD parameters from an external source.

Note that the fudge factor "wonka" is used /only/ for prettification of expected/theoretical curves in PrintASCIIHistogram displays.

38.2.7 GaussianFitHistogram

External C: Wise2_GaussianFitHistogram (h,high_hint)
Perl: &Wise2::Histogram::fit_Gaussian (h,high_hint)
Perl-OOP call: $obj->fit_Gaussian(high_hint)

Arguments

h: [UNKN ] histogram to fit [Histogram *]
high_hint: [UNKN ] score cutoff; above this are ‘realhits that arent fit [float]
returns: [UNKN ] if fit is judged to be valid else 0 if fit is invalid (too few seqs.) [int]

This function was written by Sean Eddy as part of his HMMer2 histogram.c module

Converted by Ewan Birney to Dynamite source June 98. Copyright is LGPL. For more info read READMEs

Documentation:

Fit a score histogram to a Gaussian distribution. Set the parameters mean and sd in the histogram structure, as well as a chi-squared test for goodness of fit.

38.2.8 GaussianSetHistogram

External C: Wise2_GaussianSetHistogram (h,mean,sd)
Perl: &Wise2::Histogram::set_Gaussian (h,mean,sd)
Perl-OOP call: $obj->set_Gaussian(mean,sd)

Arguments

h: [UNKN ] Undocumented argument [Histogram *]
mean: [UNKN ] Undocumented argument [float]
sd: [UNKN ] Undocumented argument [float]
returns: Nothing - no return value

This function was written by Sean Eddy as part of his HMMer2 histogram.c module

Converted by Ewan Birney to Dynamite source June 98. Copyright is LGPL. For more info read READMEs

Documentation:

Instead of fitting the histogram to a Gaussian, simply set the Gaussian parameters from an external source.

38.2.9 Evalue_from_Histogram

External C: Wise2_Evalue_from_Histogram (his,score)
Perl: &Wise2::Histogram::evalue (his,score)
Perl-OOP call: $obj->evalue(score)

Arguments

his: [UNKN ] Histogram object [Histogram *]
score: [UNKN ] score you want the evalue for [double]
returns: [UNKN ] Undocumented return value [double]

This is a convient short cut for calculating expected values from the histogram of results

39 hscore

This module contains the following objects

39.2 Hscore
39.3 DataScore
39.4 DataEntry
This module also contains some factory methods

39.1 hscore factory methods

39.1.1 std_score_Hscore

External C: Wise2_std_score_Hscore (cut_off,report_stagger)
Perl: &Wise2::std_score_Hscore (cut_off,report_stagger)

Arguments

cut_off: [UNKN ] Undocumented argument [int]
report_stagger: [UNKN ] Undocumented argument [int]
returns: [UNKN ] Undocumented return value [Hscore *]

This gives you a standard Hscore module with a cutoff in score

39.2 Object Hscore

The Hscore object has the following fields. To see how to access them refer to 7.1

: ds Type [DataScore ** : List] No documentation
: store Type [DataScoreStorage ** : List] No documentation
: his Type [Histogram * : Scalar] No documentation
: score_level Type [double : Scalar] passed into should_store function
: should_store Type [boolean (*should_store)(int given_score,double internal_score_level) : Scalar] No documentation
: score_to_his Type [float (*score_to_his)(int given_score) : Scalar] No documentation
: report_level Type [int : Scalar] number of sequences to report on
: total Type [long : Scalar] total number of scores (duplicated info in histogram)

Holds the information about a db search. Meant to be very lightweight

The histogram is carried for on-the-fly histogram storage outside of the database. The idea is that the function should_store will tell whether the datascore structure should be stored (if it is NULL, it is always stored). The score_to_his function maps the score in datascore to the float in Histogram, allowing the scoring system of the search method to be on a different basis to the scoring system of the histogram. For most times, this is going to be Score2Bits

To prevent too much dependency, the śtandardẃay of making a histogram that has a bits cut off level is done using functions in the dynlibcross module (cross code), as it needs both Hscore and Probability. You should read dynlibcross module for the constructors for Hscore

Member functions of Hscore

39.2.1 minimum_score_Hscore

External C: Wise2_minimum_score_Hscore (hs)
Perl: &Wise2::Hscore::minimum_score_Hscore (hs)
Perl-OOP call: $obj->minimum_score_Hscore()

Arguments

hs: [UNKN ] Undocumented argument [Hscore *]
returns: [UNKN ] Undocumented return value [int]

gets the minimum score from Hscore

39.2.2 maximum_score_Hscore

External C: Wise2_maximum_score_Hscore (hs)
Perl: &Wise2::Hscore::maximum_score_Hscore (hs)
Perl-OOP call: $obj->maximum_score_Hscore()

Arguments

hs: [UNKN ] Undocumented argument [Hscore *]
returns: [UNKN ] Undocumented return value [int]

gets the maximum score from Hscore

39.2.3 sort_Hscore_by_score

External C: Wise2_sort_Hscore_by_score (hs)
Perl: &Wise2::Hscore::sort_Hscore_by_score (hs)
Perl-OOP call: $obj->sort_Hscore_by_score()

Arguments

hs: [UNKN ] Hscore to be sorted [Hscore *]
returns: Nothing - no return value

As it says, sorts the high score by its score

39.2.4 length_datascore_Hscore

External C: Wise2_length_datascore_Hscore (obj)
Perl: &Wise2::Hscore::length (obj)
Perl-OOP call: $obj->length()

Arguments

obj: [READ ] Hscore object [Hscore *]
returns: [UNKN ] Undocumented return value [int]

Returns the number of datascores in the hscore structure

39.2.5 get_datascore_Hscore

External C: Wise2_get_datascore_Hscore (hs,i)
Perl: &Wise2::Hscore::datascore (hs,i)
Perl-OOP call: $obj->datascore(i)

Arguments

hs: [READ ] Hscore object [Hscore *]
i: [UNKN ] position to be read [int]
returns: [OWNER] New datascore object [DataScore *]

Returns the specific datascore held at this position.

This requires a considerable amount of memory duplication, so please dont process all your results by looping through this.

39.2.6 get_score_Hscore

External C: Wise2_get_score_Hscore (hs,i)
Perl: &Wise2::Hscore::score (hs,i)
Perl-OOP call: $obj->score(i)

Arguments

hs: [READ ] Hscore object [Hscore *]
i: [UNKN ] position to be read [int]
returns: [UNKN ] score [int]

Returns the score of the specific datascore held at this position.

39.2.7 get_evalue_Hscore

External C: Wise2_get_evalue_Hscore (hs,i)
Perl: &Wise2::Hscore::evalue (hs,i)
Perl-OOP call: $obj->evalue(i)

Arguments

hs: [READ ] Hscore object [Hscore *]
i: [UNKN ] position to be read [int]
returns: [UNKN ] evalue [double]

Returns the evalue of the specific datascore held at this position.

39.2.8 basic_show_Hscore

External C: Wise2_basic_show_Hscore (hs,ofp)
Perl: &Wise2::Hscore::show (hs,ofp)
Perl-OOP call: $obj->show(ofp)

Arguments

hs: [UNKN ] Undocumented argument [Hscore *]
ofp: [UNKN ] Undocumented argument [FILE *]
returns: Nothing - no return value

The most baby-talk showing of Hscore

39.3 Object DataScore

The DataScore object has the following fields. To see how to access them refer to 7.1

: query Type [DataEntry * : Scalar] No documentation
: target Type [DataEntry * : Scalar] No documentation
: score Type [int : Scalar] No documentation
: evalue Type [double : Scalar] No documentation
: is_stored Type [int : Scalar] No documentation

The basic one entry vs one entry structure. Each of the DataEntry datastructures actually store the information about the indexing etc.

Member functions of DataScore

39.4 Object DataEntry

The DataEntry object has the following fields. To see how to access them refer to 7.1

: name Type [char * : Scalar] name of the entry
: data[DATAENTRYSTDPOINTS] Type [int : Scalar] space for algorithms to use
: is_reversed Type [boolean : Scalar] for sequences. handy
: byte_position Type [BytePosition : Scalar] useful for indexers - hopefully long enough!
: filename Type [char * : Scalar] useful for indexers etc.

A lightweight structure to represent the information a db search algorithm will want to store and *nothing* more about a single entry

This object will be stored twice (once for the target and once for the query) for each comparison: they probably will be on different databases and different objects.

The data field is just a number (8 at the moment) of intś available for databases to store useful information (eg, length) of the object.

A number of extra fields are provided for convience sake for indexers, including byte_position and filename.

Member functions of DataEntry

40 packaln

This module contains the following objects

40.1 PackAln
40.2 PackAlnUnit

40.1 Object PackAln

The PackAln object has the following fields. To see how to access them refer to 7.1

: pau Type [PackAlnUnit ** : List] list of PackAlnUnits from start to end
: score Type [int : Scalar] score over the entire alignment

This is the lowest-level of representation of a DP alignment, being the list of (i,j,state) triples taken through the DP matrix. The score for the transition to this point is held as well.

This object is very low level and often a much better choice for representation is in /AlnBlock objects

Member functions of PackAln

40.1.1 show_simple_PackAln

External C: Wise2_show_simple_PackAln (pal,ofp)
Perl: &Wise2::PackAln::show_simple_PackAln (pal,ofp)
Perl-OOP call: $obj->show_simple_PackAln(ofp)

Arguments

pal: [UNKN ] Undocumented argument [PackAln *]
ofp: [UNKN ] Undocumented argument [FILE *]
returns: Nothing - no return value

shows packaln with a pretty verbose debugging format

40.1.2 show_bits_and_cumlative_PackAln

External C: Wise2_show_bits_and_cumlative_PackAln (pal,ofp)
Perl: &Wise2::PackAln::show_bits_and_cumlative_PackAln (pal,ofp)
Perl-OOP call: $obj->show_bits_and_cumlative_PackAln(ofp)

Arguments

pal: [UNKN ] Undocumented argument [PackAln *]
ofp: [UNKN ] Undocumented argument [FILE *]
returns: Nothing - no return value

Shows packaln as:

i,j,state,score,bits,cumlative-score,cumlative-bits

cumlative score and cumlative bits are useful sometimes

40.2 Object PackAlnUnit

The PackAlnUnit object has the following fields. To see how to access them refer to 7.1

: i Type [int : Scalar] position in query
: j Type [int : Scalar] position in target
: state Type [int : Scalar] state in FSM
: score Type [int : Scalar] score of the transition that reached this state

Internal object for /PackAln: A single position of an alignment as the (i,j,state) triple

Member functions of PackAlnUnit

41 probability

This module only contains factory methods

41.1 probability factory methods

41.1.1 Probability_from_average_state_occupancy

External C: Wise2_Probability_from_average_state_occupancy (length)
Perl: &Wise2::Probability_from_average_state_occupancy (length)

Arguments

length: [UNKN ] average length of state [double]
returns: [UNKN ] Undocumented return value [Probability]

for single state (exponetial decays) takes an average length and converts that to a probability that will produce that length (on average) for the state. NB... this *assumes* that you want a single state exp decay.

41.1.2 state_occupancy_from_Probability

External C: Wise2_state_occupancy_from_Probability (p)
Perl: &Wise2::state_occupancy_from_Probability (p)

Arguments

p: [UNKN ] probability of staying in the state [double]
returns: [UNKN ] Undocumented return value [double]

If you have a single state then this will tak the probability for the state->state transition and give you back the average length in the state

41.1.3 Probability2Score

External C: Wise2_Probability2Score (p)
Perl: &Wise2::Probability2Score (p)

Arguments

p: [UNKN ] Undocumented argument [Probability]
returns: [UNKN ] Undocumented return value [Score]

maps probabilities to scores. Deals sensibly with 0ś.

41.1.4 Score2Probability

External C: Wise2_Score2Probability (s)
Perl: &Wise2::Score2Probability (s)

Arguments

s: [UNKN ] Undocumented argument [Score]
returns: [UNKN ] Undocumented return value [Probability]

maps scores to probabilities

41.1.5 Score2Bits

External C: Wise2_Score2Bits (s)
Perl: &Wise2::Score2Bits (s)

Arguments

s: [UNKN ] Undocumented argument [Score]
returns: [UNKN ] Undocumented return value [Bits]

maps scores to bits

41.1.6 halfbit2Probability

External C: Wise2_halfbit2Probability (half_bit)
Perl: &Wise2::halfbit2Probability (half_bit)

Arguments

half_bit: [UNKN ] Undocumented argument [double]
returns: [UNKN ] Undocumented return value [Probability]

maps halfbits (log2(prob*2) to probabilities

42 protein

This module contains the following objects

42.1 Protein

42.1 Object Protein

The Protein object has the following fields. To see how to access them refer to 7.1

: baseseq Type [Sequence * : Scalar] No documentation

The protein object is a typed example of a sequence object.

It does nothing more than a sequence object but is typed

Member functions of Protein

42.1.1 Protein_from_Sequence

External C: Wise2_Protein_from_Sequence (seq)
Perl: &Wise2::Protein::Protein_from_Sequence (seq)
Perl-OOP call: $obj->Protein_from_Sequence()

Arguments

seq: [OWNER] Sequence to make protein from [Sequence *]
returns: [UNKN ] Undocumented return value [Protein *]

makes a new protein from a Sequence. It owns the Sequence memory, ie will attempt a /free_Sequence on the structure when /free_Protein is called

If you want to give this protein this Sequence and forget about it, then just hand it this sequence and set seq to NULL (no need to free it). If you intend to use the sequecne object elsewhere outside of the Protein datastructure then use Protein_from_Sequence(/hard_link_Sequence(seq))

43 proteindb

This module contains the following objects

43.2 ProteinDB
This module also contains some factory methods

43.1 proteindb factory methods

43.1.1 new_ProteinDB_from_single_seq

External C: Wise2_new_ProteinDB_from_single_seq (seq)
Perl: &Wise2::new_ProteinDB_from_single_seq (seq)

Arguments

seq: [UNKN ] sequence which as placed into ProteinDB structure. [Sequence *]
returns: [UNKN ] Undocumented return value [ProteinDB *]

To make a new protein database from a single Sequence with default amino acid mapping

43.1.2 single_fasta_ProteinDB

External C: Wise2_single_fasta_ProteinDB (filename)
Perl: &Wise2::single_fasta_ProteinDB (filename)

Arguments

filename: [UNKN ] name of fasta file [char *]
returns: [UNKN ] Undocumented return value [ProteinDB *]

pre-packed single fasta protein database

43.1.3 new_ProteinDB

External C: Wise2_new_ProteinDB (seqdb,cses)
Perl: &Wise2::new_ProteinDB (seqdb,cses)

Arguments

seqdb: [UNKN ] sequence database [SequenceDB *]
cses: [UNKN ] protein evaluation set [ComplexSequenceEvalSet *]
returns: [UNKN ] Undocumented return value [ProteinDB *]

To make a new protein database

43.2 Object ProteinDB

The ProteinDB object has the following fields. To see how to access them refer to 7.1

: is_single_seq Type [boolean : Scalar] No documentation
: is_random_db Type [boolean : Scalar] No documentation
: single Type [ComplexSequence * : Scalar] No documentation
: sdb Type [SequenceDB * : Scalar] No documentation
: cses Type [ComplexSequenceEvalSet * : Scalar] No documentation
: rnd Type [RandomProteinDB * : Scalar] No documentation
: test_dna Type [boolean : Scalar] No documentation

A database of proteins. This takes either a single sequence or a sequence database and allows a method expecting protein complexsequences to loop over it.

It also provides generic database indexing for this

Horrible physical dependency in this generated by the random protein db requiring histogram and randommodel stuff. Yuk!

Member functions of ProteinDB

44 randomdb

This module contains the following objects

44.1 RandomProteinDB
44.2 RandomDNADB

44.1 Object RandomProteinDB

The RandomProteinDB object has the following fields. To see how to access them refer to 7.1

: use_flat_length Type [boolean : Scalar] No documentation
: length Type [int : Scalar] No documentation
: length_dist Type [Histogram * : Scalar] No documentation
: emission Type [RandomModel * : Scalar] No documentation
: num Type [int : Scalar] No documentation

No documentation for RandomProteinDB

Member functions of RandomProteinDB

44.2 Object RandomDNADB

The RandomDNADB object has the following fields. To see how to access them refer to 7.1

: use_flat_length Type [boolean : Scalar] No documentation
: length Type [int : Scalar] No documentation
: length_dist Type [Histogram * : Scalar] No documentation
: emission Type [RandomModelDNA * : Scalar] No documentation
: num Type [int : Scalar] No documentation

No documentation for RandomDNADB

Member functions of RandomDNADB

45 randommodel

This module contains the following objects

45.2 RandomModelDNA
45.3 RandomModel
This module also contains some factory methods

45.1 randommodel factory methods

45.1.1 RandomModelDNA_std

External C: Wise2_RandomModelDNA_std (void)
Perl: &Wise2::RandomModelDNA_std ()

Arguments

returns: [UNKN ] Undocumented return value [RandomModelDNA *]

Returns a structure with 0.25 set in each place

45.1.2 default_RandomModel

External C: Wise2_default_RandomModel (void)
Perl: &Wise2::default_RandomModel ()

Arguments

returns: [UNKN ] Undocumented return value [RandomModel *]

Gives a default random model numbers from swissprot34- via the HMMEr1 package

45.2 Object RandomModelDNA

The RandomModelDNA object has the following fields. To see how to access them refer to 7.1

: base[5] Type [Probability : Scalar] No documentation
: name Type [char * : Scalar] No documentation

No documentation for RandomModelDNA

Member functions of RandomModelDNA

45.3 Object RandomModel

The RandomModel object has the following fields. To see how to access them refer to 7.1

: aminoacid[26] Type [Probability : Scalar] No documentation
: name Type [char * : Scalar] No documentation

No documentation for RandomModel

Member functions of RandomModel

46 sequence

This module contains the following objects

46.2 Sequence
46.3 SequenceSet
This module also contains some factory methods

46.1 sequence factory methods

46.1.1 Sequence_type_to_string

External C: Wise2_Sequence_type_to_string (type)
Perl: &Wise2::Sequence_type_to_string (type)

Arguments

type: [UNKN ] type eg SEQUENCE_PROTEIN [int]
returns: [UNKN ] Undocumented return value [char *]

Converts sequence type (SEQUENCE_*) to a string

46.1.2 new_Sequence_from_strings

External C: Wise2_new_Sequence_from_strings (name,seq)
Perl: &Wise2::new_Sequence_from_strings (name,seq)

Arguments

name: [READ ] name of sequence, memory is allocated for it. [char *]
seq: [READ ] char * of sequence, memory is allocated for it. [char *]
returns: [UNKN ] Undocumented return value [Sequence *]

Makes a new sequence from strings given. Separate memory will be allocated for them and them copied into it.

They can be NULL, in which case o a dummy name SequenceName will be assigned o No sequence placed and length of zero.

Though this is dangerous later on.

The sequence type is calculated automatically using /best_guess_type. If you want a DNA sequence but are unsure of the content of, for example, IUPAC codes, please use /force_to_dna_Sequence before using the sequence. Most of the rest of dynamite relies on a five letter A,T,G,C,N alphabet, but this function will allow any sequence type to be stored, so please check if you want to save yourself alot of grief.

In perl and other interfaces, this is a much safer constructor than the raw "new" type

46.2 Object Sequence

The Sequence object has the following fields. To see how to access them refer to 7.1

: name Type [char * : Scalar] name of the sequence
: seq Type [char * : Scalar] actual sequence
: len Type [int : Scalar] length of the sequence
: maxlen Type [int : Scalar] internal counter, indicating how much space in seq there is
: offset Type [int : Scalar] start (in bio-coords) of the sequence. Not called start due to weird legacy.
: end Type [int : Scalar] end (in bio-coords == C coords) of the sequence
: type Type [char : Scalar] guess of protein/dna type
: tax_id Type [int : Scalar] taxonimic id of this
: weight Type [double : Scalar] No documentation
: desc Type [char * : Scalar] description - often this will be NULL

This object is the basic sequence object, trying to hold little more than the name and sequence of the DNA/protein.

The len/maxlen is the actual length of the sequence (strlen(obj->seq)) and amount of memory allocated in obj->seq mainly for parsing purposes.

Member functions of Sequence

46.2.1 uppercase_Sequence

External C: Wise2_uppercase_Sequence (seq)
Perl: &Wise2::Sequence::uppercase (seq)
Perl-OOP call: $obj->uppercase()

Arguments

seq: [RW ] Sequence to be uppercased [Sequence *]
returns: Nothing - no return value

makes all the sequence uppercase

46.2.2 force_to_dna_Sequence

External C: Wise2_force_to_dna_Sequence (seq,fraction,number_of_conver)
Perl: &Wise2::Sequence::force_to_dna (seq,fraction)
Perl-OOP call: $obj->force_to_dna(fraction)

Arguments

seq: [RW ] sequence object read and converted [Sequence *]
fraction: [READ ] number 0..1 for tolerance of conversion [double]
number_of_conver: only for C api [WRITE] number of conversions actually made [int *]
returns: [READ ] TRUE for conversion to DNA, FALSE if not [boolean]

This

 a) sees how many non ATGCN characters there are in Seq
 b) If the level is below fraction
    a) flips non ATGC chars to N
    b) writes number of conversions to number\_of\_conver
    c) returns TRUE
 c) else returns FALSE

fraction of 0.0 means completely intolerant of errors fraction of 1.0 means completely tolerant of errors

46.2.3 is_reversed_Sequence

External C: Wise2_is_reversed_Sequence (seq)
Perl: &Wise2::Sequence::is_reversed (seq)
Perl-OOP call: $obj->is_reversed()

Arguments

seq: [READ ] sequence to test [Sequence *]
returns: [UNKN ] Undocumented return value [boolean]

Currently the sequence object stores reversed sequences as start > end.

This tests that and returns true if it is

46.2.4 translate_Sequence

External C: Wise2_translate_Sequence (dna,ct)
Perl: &Wise2::Sequence::translate (dna,ct)
Perl-OOP call: $obj->translate(ct)

Arguments

dna: [READ ] DNA sequence to be translated [Sequence *]
ct: [READ ] Codon table to do codon->aa mapping [CodonTable *]
returns: [OWNER] new protein sequence [Sequence *]

This translates a DNA sequence to a protein. It assummes that it starts at first residue (use trunc_Sequence to chop a sequence up).

46.2.5 reverse_complement_Sequence

External C: Wise2_reverse_complement_Sequence (seq)
Perl: &Wise2::Sequence::revcomp (seq)
Perl-OOP call: $obj->revcomp()

Arguments

seq: [READ ] Sequence to that is used to reverse (makes a new Sequence) [Sequence *]
returns: [OWNER] new Sequence which is reversed [Sequence *]

This both complements and reverses a sequence, - a common wish!

The start/end are correct with respect to the start/end of the sequence (ie start = end, end = start).

46.2.6 magic_trunc_Sequence

External C: Wise2_magic_trunc_Sequence (seq,start,end)
Perl: &Wise2::Sequence::magic_trunc (seq,start,end)
Perl-OOP call: $obj->magic_trunc(start,end)

Arguments

seq: [READ ] sequence that is the source to be truncated [Sequence *]
start: [READ ] start point [int]
end: [READ ] end point [int]
returns: [OWNER] new Sequence which is truncated/reversed [Sequence *]

Clever function for dna sequences.

When start < end, truncates normally

when start > end, truncates end,start and then reverse complements.

ie. If you have a coordinate system where reverse sequences are labelled in reverse start/end way, then this routine produces the correct sequence.

46.2.7 trunc_Sequence

External C: Wise2_trunc_Sequence (seq,start,end)
Perl: &Wise2::Sequence::trunc (seq,start,end)
Perl-OOP call: $obj->trunc(start,end)

Arguments

seq: [READ ] object holding the sequence to be truncated [Sequence *]
start: [READ ] start point of truncation [int]
end: [READ ] end point of truncation [int]
returns: [OWNER] newly allocated sequence structure [Sequence *]

truncates a sequence. It produces a new memory structure which is filled from sequence start to end.

Please notice

  
  Truncation is in C coordinates. That is

the first residue is 0 and end is the number of the residue after the cut-point. In otherwords to 2 - 3 would be a single residue truncation. So - if you want to work in more usual, ínclusiveḿolecular biology numbers, which start at 1, then you need to say

  trunc\_Sequence(seq,start-1,end);

(NB, should be (end - 1 + 1) = end for the last coordinate).

  Truncation occurs against the *absolute* coordinate

system of the Sequence, not the offset/end pair inside. So, this is a very bad error

 
  ** wrong code, and also leaks memory **


  tru = trunc\_Sequence(trunc\_Sequence(seq,50,80),55,75);

This the most portable way of doing this

  temp = trunc\_Sequence(seq,50,80);


  tru  = trunc\_Sequence(temp,55-temp->offset,75-temp->offset);


  free\_Sequence(temp);

46.2.8 read_fasta_file_Sequence

External C: Wise2_read_fasta_file_Sequence (filename)
Perl: &Wise2::Sequence::read_fasta_file_Sequence (filename)
Perl-OOP call: $obj->read_fasta_file_Sequence()

Arguments

filename: [READ ] filename to open [char *]
returns: [UNKN ] Undocumented return value [Sequence *]

Just a call

  a) open filename
  b) read sequence with /read\_fasta\_Sequence
  c) close file.

46.2.9 read_Sequence_EMBL_seq

External C: Wise2_read_Sequence_EMBL_seq (buffer,maxlen,ifp)
Perl: &Wise2::Sequence::read_Sequence_EMBL_seq (buffer,maxlen,ifp)
Perl-OOP call: $obj->read_Sequence_EMBL_seq(maxlen,ifp)

Arguments

buffer: [RW ] buffer containing the first line. [char *]
maxlen: [READ ] length of buffer [int]
ifp: [READ ] input file to read from [FILE *]
returns: [UNKN ] Undocumented return value [Sequence *]

reads the sequence part of an EMBL file.

This function can either take a file which starts

46.2.10 read_fasta_Sequence

External C: Wise2_read_fasta_Sequence (ifp)
Perl: &Wise2::Sequence::read_fasta_Sequence (ifp)
Perl-OOP call: $obj->read_fasta_Sequence()

Arguments

ifp: [UNKN ] Undocumented argument [FILE *]
returns: [UNKN ] Undocumented return value [Sequence *]

reads a fasta file assumming pretty standard sanity for the file layout

46.2.11 show_Sequence_residue_list

External C: Wise2_show_Sequence_residue_list (seq,start,end,ofp)
Perl: &Wise2::Sequence::show_debug (seq,start,end,ofp)
Perl-OOP call: $obj->show_debug(start,end,ofp)

Arguments

seq: [READ ] Sequence to show [Sequence *]
start: [READ ] start of list [int]
end: [READ ] end of list [int]
ofp: [UNKN ] Undocumented argument [FILE *]
returns: Nothing - no return value

shows a region of a sequence as

   124  A
   125  T

etc from start to end. The numbers are in C coordinates (ie, 0 is the first letter).

useful for debugging

46.2.12 write_fasta_Sequence

External C: Wise2_write_fasta_Sequence (seq,ofp)
Perl: &Wise2::Sequence::write_fasta (seq,ofp)
Perl-OOP call: $obj->write_fasta(ofp)

Arguments

seq: [READ ] sequence to be written [Sequence *]
ofp: [UNKN ] file to write to [FILE *]
returns: Nothing - no return value

writes a fasta file of the form >name Sequence

46.2.13 make_len_type_Sequence

External C: Wise2_make_len_type_Sequence (seq)
Perl: &Wise2::Sequence::validate (seq)
Perl-OOP call: $obj->validate()

Arguments

seq: [RW ] Sequence object [Sequence *]
returns: Nothing - no return value

makes seq->len and seq->end match the seq->seq length number.

It also checks the type of the sequence with /best_guess_type

46.3 Object SequenceSet

The SequenceSet object has the following fields. To see how to access them refer to 7.1

: set Type [Sequence ** : List] No documentation

A list of sequences. Not a database (you should use the db stuff for that!). But useful anyway

Member functions of SequenceSet

47 sequence_codon

This module only contains factory methods

47.1 sequence_codon factory methods

47.1.1 reverse_complement_Sequence

External C: Wise2_reverse_complement_Sequence (seq)
Perl: &Wise2::reverse_complement_Sequence (seq)

Arguments

seq: [READ ] Sequence to that is used to reverse (makes a new Sequence) [Sequence *]
returns: [OWNER] new Sequence which is reversed [Sequence *]

This both complements and reverses a sequence, - a common wish!

The start/end are correct with respect to the start/end of the sequence (ie start = end, end = start).

47.1.2 magic_trunc_Sequence

External C: Wise2_magic_trunc_Sequence (seq,start,end)
Perl: &Wise2::magic_trunc_Sequence (seq,start,end)

Arguments

seq: [READ ] sequence that is the source to be truncated [Sequence *]
start: [READ ] start point [int]
end: [READ ] end point [int]
returns: [OWNER] new Sequence which is truncated/reversed [Sequence *]

Clever function for dna sequences.

When start < end, truncates normally

when start > end, truncates end,start and then reverse complements.

ie. If you have a coordinate system where reverse sequences are labelled in reverse start/end way, then this routine produces the correct sequence.

47.1.3 translate_Sequence

External C: Wise2_translate_Sequence (dna,ct)
Perl: &Wise2::translate (dna,ct)

Arguments

dna: [READ ] DNA sequence to be translated [Sequence *]
ct: [READ ] Codon table to do codon->aa mapping [CodonTable *]
returns: [OWNER] new protein sequence [Sequence *]

This translates a DNA sequence to a protein. It assummes that it starts at first residue (use trunc_Sequence to chop a sequence up).

48 sequencedb

This module contains the following objects

48.2 SequenceDB
48.3 FileSource
This module also contains some factory methods

48.1 sequencedb factory methods

48.1.1 single_fasta_SequenceDB

External C: Wise2_single_fasta_SequenceDB (filename)
Perl: &Wise2::single_fasta_SequenceDB (filename)

Arguments

filename: [UNKN ] name of fastadb [char *]
returns: [UNKN ] Undocumented return value [SequenceDB *]

pre-packed single fasta file db

48.2 Object SequenceDB

The SequenceDB object has the following fields. To see how to access them refer to 7.1

: name Type [char * : Scalar] No documentation
: fs Type [FileSource ** : List] No documentation
: current_source Type [int : Scalar] No documentation
: current_file Type [FILE * : Scalar] No documentation
: sequence_no Type [int : Scalar] No documentation
: byte_position Type [int : Scalar] No documentation
: has_warned_single Type [int : Scalar] No documentation
: seq_start Type [int : Scalar] No documentation
: seq_end Type [int : Scalar] No documentation

This is the basic Sequence database wrapper - it handles all the formats and the on-the-fly indexing.

Generally it wont be directly used by an algorithm, which will be using something specific to the sequence type produce complex sequence type objects: ie you will be using proteindb or cdnadb which internally will be using this object

Member functions of SequenceDB

48.2.1 close_SequenceDB

External C: Wise2_close_SequenceDB (last,sdb)
Perl: &Wise2::SequenceDB::close_SequenceDB (last,sdb)
Perl-OOP call: $obj->close_SequenceDB(sdb)

Arguments

last: [WRITE] Sequence object to be freed [Sequence *]
sdb: [READ ] database to be closed [SequenceDB *]
returns: [UNKN ] Undocumented return value [boolean]

top level function that closes the SequenceDB after the last sequence is read.

48.3 Object FileSource

The FileSource object has the following fields. To see how to access them refer to 7.1

: filename Type [char * : Scalar] No documentation
: input Type [FILE * : Scalar] could be stdin!
: format Type [int : Scalar] No documentation
: type Type [int : Scalar] No documentation

This object represents a single file source for a database. At the moment only multiple fasta files are catered for

Member functions of FileSource

49 transcript

This module contains the following objects

49.1 Exon
49.2 Transcript

49.1 Object Exon

The Exon object has the following fields. To see how to access them refer to 7.1

: start Type [int : Scalar] No documentation
: end Type [int : Scalar] No documentation
: used Type [boolean : Scalar] used by some prediction programs etc
: score Type [double : Scalar] No documentation
: sf Type [SupportingFeature ** : List] No documentation
: phase Type [int : Scalar] No documentation

No documentation for Exon

Member functions of Exon

49.2 Object Transcript

The Transcript object has the following fields. To see how to access them refer to 7.1

: exon Type [Exon ** : List] No documentation
: parent Type [Gene * : Scalar] No documentation
: translation Type [Translation ** : List] No documentation
: cDNA Type [cDNA * : Scalar] may not be here!

Transcript represents a single spliced product from a gene. The transcript is considered to be a series of exons and it contains, in addition a series of translations. Most genes will only have one translation.

Like gene before it, transcript does not necessarily contain DNA. When some DNA is asked from it, via get_cDNA_from_Transcript (notice the change from Genomic typed sequence in Gene to cDNA typed sequence in Transcript) it first checkes the ćache. If it is not there, it asks for its parents genomic DNA, and then assemblies the cDNA using the exon coordinates. The exon coordinates are such that 0 means the start of the gene, not the start of the genomic region. (makes some outputs a pain).

Supporting Features are added to exons, and, keeping in the spirit of this module, are relative to the exon start. The strand is inherieted from the exon

Member functions of Transcript

49.2.1 get_cDNA_from_Transcript

External C: Wise2_get_cDNA_from_Transcript (trs)
Perl: &Wise2::Transcript::get_cDNA_from_Transcript (trs)
Perl-OOP call: $obj->get_cDNA_from_Transcript()

Arguments

trs: [READ ] transcript to get cDNA from [Transcript *]
returns: [SOFT ] cDNA of the transcript [cDNA *]

gets the cDNA associated with this transcript, if necessary, building it from the exon information provided.

returns a soft-linked object. If you want to ensure that this cDNA object remains in memory use /hard_link_cDNA on the object.

50 translation

This module contains the following objects

50.1 Translation

50.1 Object Translation

The Translation object has the following fields. To see how to access them refer to 7.1

: start Type [int : Scalar] No documentation
: end Type [int : Scalar] No documentation
: parent Type [Transcript * : Scalar] No documentation
: protein Type [Protein * : Scalar] No documentation

Translation represents a single translation from a cDNA. Although most cDNAs will have one translation, this does allow us to deal with alternative translation points etc.

As with Transcript and Gene before it, the translation does not necessarily have any sequence in it. When sequence is asked for by get_Protein_from_Translation() the cache is checked, and if it is empty, then the transcriptś DNA is called for, and the converted into the translation with appropiate start and stops.

Of course, get_Protein_from_Translation can potentially trigger the construction of an entire gene upstairs, but that need not worry you here

Member functions of Translation

50.1.1 get_Protein_from_Translation

External C: Wise2_get_Protein_from_Translation (ts,ct)
Perl: &Wise2::Translation::get_Protein_from_Translation (ts,ct)
Perl-OOP call: $obj->get_Protein_from_Translation(ct)

Arguments

ts: [UNKN ] translation [Translation *]
ct: [UNKN ] codon table to use [CodonTable *]
returns: [SOFT ] Protein sequence [Protein *]

Gets the protein

This document was translated from L^AT_EX by H^EV^EA.

Wise2 API (version 2.1.19b)

Ewan Birney Sanger Centre Wellcome Trust Genome Campus Hinxton, Cambridge CB10 1SA, England. Email: birney@sanger.ac.uk

Contents

1 Overview

2 WARNING - still in alpha

3 API generation

4 Getting Started for the impatient

4.1 Perl reverse complement

4.2 Wise2 external API calls

4.3 Wise2 internal API calls

5 Navigating the source code

5.1 Making a translation of a DNA sequence

5.2 Comparing two sequences using smith waterman

5.3 Running a smith waterman search of single protein sequence vs a db

5.4 Running a genewise on a single protein vs a single DNA sequence

6 Concepts and overview of the API

7 The reference section

7.1 Accessing fields in objects

7.1.1 Perl scalar accessors

7.1.2 External C scalar accessors

7.1.3 Perl List accessors

7.1.4 External API List accessors

7.2 Object Construction and handling

7.2.1 Low level Object Constructing

7.2.2 Object deconstructors

8 Wise2 Specific Modules

9 Dynamite library modules

9.1 Sequence modules

9.2 Generic probabilistic modelling support

9.3 Generic Database Searching

9.4 Generic Dynamite algorithm support

10 cdparser

10.1 cdparser factory methods

10.1.1 flat_cDNAParser

10.2 Object cDNAParser

11 dnaalign

11.1 dnaalign factory methods

11.1.1 make_align_dnaalign

11.1.2 DnaStartEnd_from_policy

11.2 Object DnaStartEnd

12 estwrap

12.1 estwrap factory methods

12.1.1 Hscore_from_TSM_estwise

12.1.2 AlnBlock_from_Protein_estwise_wrap

12.1.3 AlnBlock_from_TSM_estwise_wrap

12.1.4 alg_estwrap_from_string

13 genedisplay

13.1 genedisplay factory methods

13.1.1 protein2genomic_ascii_display

13.1.2 protcdna_ascii_display

14 genefrequency

14.1 genefrequency factory methods

14.1.1 read_GeneFrequency21_file

14.1.2 read_GeneFrequency21

14.2 Object GeneFrequency21

14.3 Object GeneConsensus

14.4 Object GeneSingleCons

15 geneparameter

15.1 Object GeneParameter21

16 matchsum

16.1 Object MatchSummarySet

16.1.1 MatchSummarySet_from_AlnBlock_estwise

16.1.2 MatchSummarySet_from_AlnBlock_genewise

16.2 Object MatchSummary

17 pfamhmmer1db

17.1 Object PfamHmmer1DB

17.2 Object PfamHmmer1Entry

18 seqaligndisplay

18.1 seqaligndisplay factory methods

18.1.1 write_pretty_str_align

18.1.2 write_pretty_seq_align

18.1.3 write_pretty_Protein_align

19 seqhit

19.1 Object DnaSequenceHitList

19.1.1 show_DnaSequenceHitList

19.1.2 read_MSPcrunch_DnaSequenceHitList

19.2 Object SegmentHitList

19.3 Object SegmentHit

20 sw_wrap

20.1 sw_wrap factory methods

Ewan Birney
Sanger Centre
Wellcome Trust Genome Campus
Hinxton, Cambridge CB10 1SA,
England.
Email: birney@sanger.ac.uk