embest.c


Function embEstSetDebug

Sets debugging calls on

Prototype

void embEstSetDebug (
      void
);

TypeNameRead/WriteDescription
void RETURN

From EMBOSS 1.13.0


Function embEstSetVerbose

Sets verbose debugging calls on

Prototype

void embEstSetVerbose (
      void
);

TypeNameRead/WriteDescription
void RETURN

From EMBOSS 1.13.0


Function embEstGetSeed

Returns a seed for the random number generator, using the system clock.

Prototype

ajint embEstGetSeed (
      void
);

TypeNameRead/WriteDescription
ajint RETURNseed.

From EMBOSS 1.13.0


Function embEstMatInit

Comparison matrix initialisation.

Prototype

void embEstMatInit (
      ajint match,
      ajint mismatch,
      ajint gap,
      ajint neutral,
      char pad_char
);

TypeNameRead/WriteDescription
ajintmatchInputMatch code
ajintmismatchInputMismatch penalty
ajintgapInputGap penalty
ajintneutralInputScore for ambiguous base positions.
charpad_charInputPad character for gaps in input sequences
void RETURN

From EMBOSS 1.13.0


Function embEstFindSpliceSites

Finds all putative DONOR and ACCEPTOR splice sites in the genomic sequence.

Prototype

AjPSeq embEstFindSpliceSites (
      const AjPSeq genome,
      ajint forward 
);

TypeNameRead/WriteDescription
const AjPSeqgenomeInputGenomic sequence
ajintforwardInputBoolean. 1 = forward direction
AjPSeq RETURNSequence of bitmask codes for splice sites.

From EMBOSS 1.13.0


Function embEstShuffleSeq

Shuffle the sequence.

Prototype

AjPSeq embEstShuffleSeq (
      AjPSeq seq,
      ajint in_place,
      ajint* seed 
);

TypeNameRead/WriteDescription
AjPSeqseqModifyOriginal sequence
ajintin_placeInputBoolean 1=shuffle in place
ajint*seedModifyRandom number seed.
AjPSeq RETURNshuffled sequence.

From EMBOSS 1.13.0


Function embEstFreeAlign

Free a genomic EST alignment structure

Prototype

void embEstFreeAlign (
      EmbPEstAlign* ge
);

TypeNameRead/WriteDescription
EmbPEstAlign*geDeleteGenomic EST alignment data structure
void RETURN

From EMBOSS 1.13.0


Function embEstPrintAlign

Print the alignment

Prototype

void embEstPrintAlign (
      AjPFile ofile,
      const AjPSeq genome,
      const AjPSeq est,
      const EmbPEstAlign ge,
      ajint width 
);

TypeNameRead/WriteDescription
AjPFileofileModifyOutput file
const AjPSeqgenomeInputGenomic sequence
const AjPSeqestInputEST sequence
const EmbPEstAligngeInputGenomic EST alignment
ajintwidthInputOutput width (in bases)
void RETURN

From EMBOSS 1.13.0


Function embEstAlignNonRecursive

Modified Smith-Waterman/Needleman to align an EST or mRNA to a Genomic sequence, allowing for introns.

Prototype

EmbPEstAlign embEstAlignNonRecursive (
      const AjPSeq est,
      const AjPSeq genome,
      ajint gap_penalty,
      ajint intron_penalty,
      ajint splice_penalty,
      const AjPSeq splice_sites,
      ajint backtrack,
      ajint needleman,
      ajint init_path
);

TypeNameRead/WriteDescription
const AjPSeqestInputSequence of EST
const AjPSeqgenomeInputSequence of genomic region
ajintgap_penaltyInputGap penalty
ajintintron_penaltyInputIntron penalty
ajintsplice_penaltyInputSplice site penalty
const AjPSeqsplice_sitesInputMarked splice sites. The intron_penalty may be modified to splice_penalty if splice_sites is non-null and there are DONOR and ACCEPTOR sites at the start and end of the intron.
ajintbacktrackInputBoolean. If backtrack is 0 then only the start and end points and the score are computed, and no path matrix is allocated.
ajintneedlemanInputBoolean 1 = global alignment 0 = local alignment
ajintinit_pathInputType of initialization for the path. If init_path is DIAGONAL then the boundary conditions are adjusted so that the optimal path enters the cell (0,0) diagonally. Otherwise it enters from the left (ie as a deletion in the EST)
EmbPEstAlign RETURNResulting genomic EST alignment

From EMBOSS 1.13.0


Function embEstAlignLinearSpace

Align EST sequence to genomic in linear space

Prototype

EmbPEstAlign embEstAlignLinearSpace (
      const AjPSeq est,
      const AjPSeq genome,
      ajint match,
      ajint mismatch,
      ajint gap_penalty,
      ajint intron_penalty,
      ajint splice_penalty,
      const AjPSeq splice_sites,
      float megabytes 
);

TypeNameRead/WriteDescription
const AjPSeqestInputSequence of EST
const AjPSeqgenomeInputSequence of genomic region
ajintmatchInputMatch score
ajintmismatchInputMismatch penalty (positive)
ajintgap_penaltyInputGap penalty
ajintintron_penaltyInputIntron penalty
ajintsplice_penaltyInputSplice site penalty
const AjPSeqsplice_sitesInputMarked splice sites. The intron_penalty may be modified to splice_penalty if splice_sites is non-null and there are DONOR and ACCEPTOR sites at the start and end of the intron.
floatmegabytesInputMaximum memory allowed in Mbytes for alignment by standard methods.
EmbPEstAlign RETURNGenomic EST alignment

From EMBOSS 1.13.0


Function embEstOutBlastStyle

output in blast style.

Prototype

void embEstOutBlastStyle (
      AjPFile blast,
      const AjPSeq genome,
      const AjPSeq est,
      const EmbPEstAlign ge,
      ajint gap_penalty,
      ajint intron_penalty,
      ajint splice_penalty,
      ajint gapped,
      ajint reverse
);

TypeNameRead/WriteDescription
AjPFileblastModifyOutput file
const AjPSeqgenomeInputGenomic sequence
const AjPSeqestInputEST sequence
const EmbPEstAligngeInputGenomic EST alignment
ajintgap_penaltyInputGap penalty
ajintintron_penaltyInputIntron penalty
ajintsplice_penaltyInputSplice site penalty
ajintgappedInputBoolean. 1 = write a gapped alignment
ajintreverseInputBoolean. 1 = reverse alignment
void RETURN

From EMBOSS 1.13.0