inc/pairHmm.h File Reference

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Data Structures
struct	phmmMommy
struct	phmmState
struct	phmmMatrix
Defines
#define	phmmPackMommy(stateIx, qOff, tOff)   ((UBYTE)((stateIx) + ((-(qOff))<<5) + ((-(tOff))<<6)))
#define	phmmMommyTraceBit   (1<<7)
Functions
phmmMatrix *	phmmMatrixNew (int stateCount, char *query, int querySize, char *target, int targetSize)
void	phmmMatrixFree (struct phmmMatrix **pAm)
phmmState *	phmmNameState (struct phmmMatrix *am, int stateIx, char *name, char emitLetter)
phmmAliPair *	phmmTraceBack (struct phmmMatrix *am, struct phmmMommy *end)
void	phmmPrintTrace (struct phmmMatrix *am, struct phmmAliPair *pairList, boolean showStates, FILE *f, boolean extraAtEnds)
axt *	phhmTraceToAxt (struct phmmMatrix *am, struct phmmAliPair *pairList, int score, char *qName, char *tName)
Variables
UBYTE	phmmNullMommy

---

Define Documentation

#define phmmMommyTraceBit   (1<<7)

Definition at line 29 of file pairHmm.h.

Referenced by phmmTraceBack().

#define phmmPackMommy	(	stateIx,
		qOff,
		tOff		)	((UBYTE)((stateIx) + ((-(qOff))<<5) + ((-(tOff))<<6)))

Definition at line 26 of file pairHmm.h.

---

Function Documentation

struct axt* phhmTraceToAxt	(	struct phmmMatrix *	am,
		struct phmmAliPair *	pairList,
		int	score,
		char *	qName,
		char *	tName
	)			[read]

Definition at line 304 of file pairHmm.c.

References AllocArray, AllocVar, cloneString(), phmmAliPair::next, axt::qEnd, axt::qName, axt::qStart, axt::qStrand, axt::qSym, phmmAliPair::queryIx, phmmAliPair::querySym, axt::score, slCount(), axt::symCount, phmmAliPair::targetIx, phmmAliPair::targetSym, axt::tEnd, axt::tName, axt::tStart, axt::tStrand, and axt::tSym.

Referenced by axtAffine().

{
struct axt *axt;
struct phmmAliPair *pair;
int qEnd, tEnd;
char *qSym, *tSym;
int i;

/* Make sure got something real. */
if ((pair = pairList) == NULL)
    return NULL;

/* Allocate memory for axt. */
AllocVar(axt);
axt->symCount = slCount(pairList);
axt->qSym = AllocArray(qSym, axt->symCount+1);
axt->tSym = AllocArray(tSym, axt->symCount+1);

/* Fill in basic fields. */
axt->qName = cloneString(qName);
axt->tName = cloneString(tName);
axt->qStrand = axt->tStrand = '+';
axt->qStart = qEnd = pair->queryIx;
axt->tStart = tEnd = pair->targetIx;
axt->score = score;

/* Store alignment symbols and keep track of last symbol used. */
for (i=0, pair = pairList; pair != NULL; pair = pair->next, ++i)
    {
    qSym[i] = pair->querySym;
    tSym[i] = pair->targetSym;
    qEnd = pair->queryIx;
    tEnd = pair->targetIx;
    }

/* Store end and return. */
axt->qEnd = qEnd + 1;
axt->tEnd = tEnd + 1;
return axt;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void phmmMatrixFree

(

struct phmmMatrix **

pAm

)

Definition at line 70 of file pairHmm.c.

References phmmMatrix::allCells, phmmMatrix::allScores, freeMem(), freez(), and phmmMatrix::states.

Referenced by axtAffine().

{
struct phmmMatrix *am = *pAm;
if (am != NULL)
    {
    freeMem(am->states);
    freeMem(am->allCells);
    freeMem(am->allScores);
    freez(pAm);
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

struct phmmMatrix* phmmMatrixNew	(	int	stateCount,
		char *	query,
		int	querySize,
		char *	target,
		int	targetSize
	)			[read]

Definition at line 24 of file pairHmm.c.

References phmmMatrix::allCells, AllocVar, phmmMatrix::allScores, phmmState::cells, phmmState::lastScores, needLargeMem(), needMem(), phmmMatrix::qDim, phmmMatrix::query, phmmMatrix::querySize, phmmState::scores, phmmMatrix::stateByteSize, phmmMatrix::stateCount, phmmState::stateIx, phmmMatrix::states, phmmMatrix::stateSize, phmmMatrix::target, phmmMatrix::targetSize, and phmmMatrix::tDim.

Referenced by affineAlign(), and xenAlignSmall().

{
int i;
struct phmmMommy *allCells;
int allCellSize;
int rowSize;
int *allScores;
struct phmmMatrix *am;

AllocVar(am);
am->query = query;
am->target = target;
am->querySize = querySize;
am->targetSize = targetSize;
am->qDim = rowSize = am->querySize + 1;
am->tDim = am->targetSize + 1;
am->stateCount = stateCount;
am->stateSize = rowSize * am->tDim;
am->stateByteSize = am->stateSize * sizeof(struct phmmMommy);
am->states = needMem(stateCount * sizeof(struct phmmState));

/* Initialize matrix of cells for each state. */
allCellSize = stateCount * am->stateByteSize;
am->allCells = allCells = needLargeMem(allCellSize); 
memset(allCells, 0, allCellSize);
for (i=0; i<stateCount; ++i)
    {
    am->states[i].cells = allCells;
    allCells += am->stateSize;
    am->states[i].stateIx = i;
    }

/* Initialize two rows of scores for each state. */
allScores = am->allScores = needMem(rowSize * 2 * stateCount * sizeof(int) );
for (i=0; i<stateCount; ++i)
    {
    am->states[i].scores = allScores;
    allScores += rowSize;
    am->states[i].lastScores = allScores;
    allScores += rowSize;
    }
return am;
}

Here is the call graph for this function:

Here is the caller graph for this function:

struct phmmState* phmmNameState	(	struct phmmMatrix *	am,
		int	stateIx,
		char *	name,
		char	emitLetter
	)			[read]

Definition at line 84 of file pairHmm.c.

References phmmState::emitLetter, phmmState::name, and phmmMatrix::states.

Referenced by affineAlign(), and xenAlignSmall().

{
struct phmmState *state;
assert(stateIx < am->stateCount);
state = &am->states[stateIx];
state->name = name;
state->emitLetter = emitLetter;
return state;
}

Here is the caller graph for this function:

void phmmPrintTrace	(	struct phmmMatrix *	am,
		struct phmmAliPair *	pairList,
		boolean	showStates,
		FILE *	f,
		boolean	extraAtEnds
	)

Definition at line 187 of file pairHmm.c.

References FALSE, phmmAliPair::hiddenSym, lineLen, phmmAliPair::next, phmmMatrix::query, phmmAliPair::queryIx, phmmMatrix::querySize, phmmAliPair::querySym, phmmMatrix::target, phmmAliPair::targetIx, phmmMatrix::targetSize, phmmAliPair::targetSym, and TRUE.

{
struct phmmAliPair *pair;
#define lineLen 50
char asyms[lineLen+1];
char qsyms[lineLen+1];
char tsyms[lineLen+1];
char hsyms[lineLen+1];
int lineIx = 0;
int qs, ts;
boolean gotQs = FALSE;

if ((pair = pairList) == NULL)
    {
    fprintf(f, "Empty pair list\n");
    return;
    }

qs = pair->queryIx;
ts = pair->targetIx;

/* Print out up to 25 bp of initial non-matching parts. */
if (extraAtEnds)
    {
    int qStart = qs;
    int tStart = ts;
    int i;

    for (i= -25; i < 0; ++i)
        {
        int qIx = qStart + i;
        int tIx = tStart + i;
        qsyms[lineIx] = (qIx >= 0 ? am->query[qIx] : ' ');
        tsyms[lineIx] = (tIx >= 0 ? am->target[tIx] : ' ' );
        asyms[lineIx] = hsyms[lineIx] = ' ';
        if (qIx >= 0 || tIx >= 0)
            {
            ++lineIx;
            if (!gotQs)
                {
                gotQs = TRUE;
                qs = qIx;
                ts = tIx;
                }
            }
        }
    }

/* Print out matching parts */
for (pair = pairList; pair != NULL; pair = pair->next)
    {
    qsyms[lineIx] = pair->querySym;
    tsyms[lineIx] = pair->targetSym;
    asyms[lineIx] = (pair->querySym == pair->targetSym ? '|' : ' ');
    hsyms[lineIx] = pair->hiddenSym;
    if (++lineIx == lineLen)
        {
        qsyms[lineIx] = asyms[lineIx] = tsyms[lineIx] = hsyms[lineIx] = 0;
        fprintf(f, "%5d %s\n", qs, qsyms);
        fprintf(f, "%5s %s\n", "", asyms);
        fprintf(f, "%5d %s\n", ts, tsyms);
        if (showStates)
            fprintf(f, "%5s %s\n", "", hsyms);
        fputc('\n', f);
        lineIx = 0;
        if (pair->next)
            {
            qs = pair->next->queryIx;
            ts = pair->next->targetIx;
            }
        }
    }

/* Print out last bit - first get last pair. */
if (extraAtEnds)
    {
    for (pair = pairList; pair->next != NULL; pair = pair->next)
        ;

        {
        int qIx, tIx, i;
        int residue = lineLen - lineIx;
        for (i=1; i<=residue; i++)
            {
            if ((qIx = pair->queryIx+i) >= am->querySize)
                qsyms[lineIx] = ' ';
            else
                qsyms[lineIx] = am->query[qIx];
            if ((tIx = pair->targetIx+i) >= am->targetSize)
                tsyms[lineIx] = ' ';
            else
                tsyms[lineIx] = am->target[tIx];
            asyms[lineIx] = ' ';
            hsyms[lineIx] = ' ';
            ++lineIx;
            assert(lineIx <= lineLen);
            }
        }
    }

/* Print out anything not printed out yet. */
if (lineIx != 0)
    {
    qsyms[lineIx] = asyms[lineIx] = tsyms[lineIx] = hsyms[lineIx] = 0;
    fprintf(f, "%5d %s\n", qs, qsyms);
    fprintf(f, "%5s %s\n", "", asyms);
    fprintf(f, "%5d %s\n", ts, tsyms);
    if (showStates)
        fprintf(f, "%5s %s\n", "", hsyms);
    fputc('\n', f);
    lineIx = 0;
    }
#undef lineLen
}

struct phmmAliPair* phmmTraceBack	(	struct phmmMatrix *	am,
		struct phmmMommy *	end
	)			[read]

Definition at line 134 of file pairHmm.c.

References AllocVar, phmmMommy::mommy, phmmFindMatrixIx(), phmmFindMommy(), phmmMommyTraceBit, slAddHead, and UBYTE.

{
struct phmmAliPair *pairList = NULL, *pair;
struct phmmMommy *cell, *parent = end;
int parentSix, parentTix, parentQix;
int sIx, tIx, qIx;

phmmFindMatrixIx(am, parent, &parentSix, &parentQix, &parentTix);
for (;;)
    {
    cell = parent;
    sIx = parentSix;
    tIx = parentTix;
    qIx = parentQix;

    if ((parent = phmmFindMommy(am, cell)) == NULL)
        break;
    phmmFindMatrixIx(am, parent, &parentSix, &parentQix, &parentTix);
    
    cell->mommy |= phmmMommyTraceBit;

    AllocVar(pair);
    pair->hiddenIx = (UBYTE)sIx;
    pair->hiddenSym = am->states[sIx].emitLetter; 

    if (parentQix == qIx - 1 && parentTix == tIx - 1 )
        {
        pair->queryIx = qIx-1;
        pair->querySym = am->query[qIx-1];
        pair->targetIx = tIx - 1;
        pair->targetSym = am->target[tIx-1];
        }
    else if (parentQix == qIx) /* && parentTix == tIx-1 */
        {
        pair->queryIx = -1;
        pair->querySym = '-';
        pair->targetIx = tIx-1;
        pair->targetSym = am->target[tIx-1];
        }
    else  /* parentTix == tIx  && parentQix == qIx-1 */
        {
        pair->queryIx = qIx-1;
        pair->querySym = am->query[qIx-1];
        pair->targetIx = -1;
        pair->targetSym = '-';
        }
    slAddHead(&pairList, pair);
    }
return pairList;
}

Here is the call graph for this function:

---

Variable Documentation

UBYTE phmmNullMommy

Definition at line 13 of file pairHmm.c.

Referenced by phmmFindMommy().

---