Lale Class Reference

#include <kzblast.h>

List of all members.

Public Member Functions

 Lale ()
 Lale (int ws, const BlastParameter &bp)
void readQuery (const string &file)
void buildfsa ()
void showfsa (ostream &ous) const
list< Alignseg > & findMatch (const bioseq &s)
int extendhit (Alignseg &as, const int *tc, const int tclen)
void showseg_debug (const Alignseg &as, const bioseq &ss) const
void showmaxss_debug (const bioseq *tseq) const
bool isComplete (const Alignseg &as, const bioseq &ss)
void setMatrix (const Matrix &mt)
const MatrixgetMatrix () const
const bioseqgetQuery () const

Private Member Functions

void mergeNeighbor (list< Alignseg > &segs, const int *qsc, const int *tsc)

Private Attributes

bioseqqseq
Matrix mtx
BlastParameter param
list< Alignsegmaxss
ltriefsa


Detailed Description

Local Alignment by look up and extension LALE When this object is build it constructs a default BlastParameter object as member for useage. This parameter can be altered by setting to a new value.

Look up used a trie data structure.


Constructor & Destructor Documentation

Lale::Lale (  )  [inline]

Lale::Lale ( int  ws,
const BlastParameter bp 
) [inline]


Member Function Documentation

void Lale::readQuery ( const string &  file  ) 

References bioseq::length(), qseq, and bioseq::read().

Referenced by main().

void Lale::buildfsa (  )  [inline]

void Lale::showfsa ( ostream &  ous  )  const [inline]

debug function

References fsa, and ltrie::showLeafs().

list< Alignseg > & Lale::findMatch ( const bioseq s  ) 

find matches to the target sequence s

this method should return a data structure for further processing Lale is X, target is Y

References BlastParameter::extCutoff, extendhit(), ltrie::find(), fsa, bioseq::getcode(), isComplete(), bioseq::length(), maxss, mergeNeighbor(), mtx, intnode::next, param, lnode::pos, qseq, Matrix::score(), bioseq::toString(), intnode::value, BlastParameter::wordjoin, and BlastParameter::wordSize.

Referenced by main().

int Lale::extendhit ( Alignseg as,
const int *  tc,
const int  tclen 
)

extend the hit Alignseg as. check to see the Alignseg has score > extCutoff then extend in both directions to see that the score reaches the max value

return 0 for not extended 1 extend end 2 extend start 3 extend both start and end

References Alignseg::addLength(), BlastParameter::extCutoff, BlastParameter::extDropoff, bioseq::getcode(), Alignseg::getLength(), Alignseg::getScore(), Alignseg::getX(), Alignseg::getY(), length, Matrix::lookup(), mtx, param, qseq, Alignseg::setScore(), Alignseg::setX(), and Alignseg::setY().

Referenced by findMatch().

void Lale::showseg_debug ( const Alignseg as,
const bioseq ss 
) const

void Lale::showmaxss_debug ( const bioseq tseq  )  const

References maxss, qseq, and bioseq::substring().

bool Lale::isComplete ( const Alignseg as,
const bioseq ss 
)

if it is complete, it will change the maxss to the only element that is complete

References Alignseg::getLength(), bioseq::length(), maxss, and qseq.

Referenced by findMatch().

void Lale::setMatrix ( const Matrix mt  )  [inline]

References mtx.

const Matrix* Lale::getMatrix (  )  const [inline]

References mtx.

Referenced by main().

const bioseq* Lale::getQuery (  )  const [inline]

References qseq.

Referenced by main().

void Lale::mergeNeighbor ( list< Alignseg > &  segs,
const int *  qsc,
const int *  tsc 
) [private]

References mtx, and Matrix::score().

Referenced by findMatch().


Member Data Documentation

bioseq* Lale::qseq [private]

Matrix Lale::mtx [private]

Referenced by buildfsa(), extendhit(), and findMatch().

list<Alignseg> Lale::maxss [private]

ltrie* Lale::fsa [private]

Referenced by buildfsa(), findMatch(), and showfsa().


The documentation for this class was generated from the following files:

Generated on Wed Aug 10 11:57:12 2011 for Softwares from Orpara by  doxygen 1.5.6