lib/chainToAxt.c File Reference

#include "common.h"
#include "chain.h"
#include "dnautil.h"
#include "dnaseq.h"
#include "axt.h"
#include "chainToAxt.h"

Include dependency graph for chainToAxt.c:

Go to the source code of this file.

Functions
static struct axt *	axtFromBlocks (struct chain *chain, struct cBlock *startB, struct cBlock *endB, struct dnaSeq *qSeq, int qOffset, struct dnaSeq *tSeq, int tOffset)
axt *	chainToAxt (struct chain *chain, struct dnaSeq *qSeq, int qOffset, struct dnaSeq *tSeq, int tOffset, int maxGap, int maxChain)
Variables
static char const	rcsid [] = "$Id: chainToAxt.c,v 1.4 2005/04/10 14:41:21 markd Exp $"

---

Function Documentation

static struct axt* axtFromBlocks	(	struct chain *	chain,
		struct cBlock *	startB,
		struct cBlock *	endB,
		struct dnaSeq *	qSeq,
		int	qOffset,
		struct dnaSeq *	tSeq,
		int	tOffset
	)			[static, read]

Definition at line 12 of file chainToAxt.c.

References AllocVar, axtScoreDnaDefault(), blockSize, cloneString(), dnaSeq::dna, needLargeMem(), cBlock::next, axt::qEnd, cBlock::qEnd, chain::qName, axt::qName, axt::qStart, cBlock::qStart, chain::qStrand, axt::qStrand, axt::qSym, axt::score, axt::symCount, axt::tEnd, cBlock::tEnd, chain::tName, axt::tName, axt::tStart, cBlock::tStart, axt::tStrand, and axt::tSym.

Referenced by chainToAxt().

{
int symCount = 0;
int dq, dt, blockSize = 0, symIx = 0;
struct cBlock *b, *a = NULL;
struct axt *axt;
char *qSym, *tSym;

/* Make a pass through figuring out how big output will be. */
for (b = startB; b != endB; b = b->next)
    {
    if (a != NULL)
        {
        dq = b->qStart - a->qEnd;
        dt = b->tStart - a->tEnd;
        symCount += dq + dt;
        }
    blockSize = b->qEnd - b->qStart;
    symCount += blockSize;
    a = b;
    }

/* Allocate axt and fill in most fields. */
AllocVar(axt);
axt->qName = cloneString(chain->qName);
axt->qStart = startB->qStart;
axt->qEnd = a->qEnd;
axt->qStrand = chain->qStrand;
axt->tName = cloneString(chain->tName);
axt->tStart = startB->tStart;
axt->tEnd = a->tEnd;
axt->tStrand = '+';
axt->symCount = symCount;
axt->qSym = qSym = needLargeMem(symCount+1);
qSym[symCount] = 0;
axt->tSym = tSym = needLargeMem(symCount+1);
tSym[symCount] = 0;

/* Fill in symbols. */
a = NULL;
for (b = startB; b != endB; b = b->next)
    {
    if (a != NULL)
        {
        dq = b->qStart - a->qEnd;
        dt = b->tStart - a->tEnd;
        if (dq == 0)
            {
            memset(qSym+symIx, '-', dt);
            memcpy(tSym+symIx, tSeq->dna + a->tEnd - tOffset, dt);
            symIx += dt;
            }
        else
            {
            assert(dt == 0);
            memset(tSym+symIx, '-', dq);
            memcpy(qSym+symIx, qSeq->dna + a->qEnd - qOffset, dq);
            symIx += dq;
            }
        }
    blockSize = b->qEnd - b->qStart;
    memcpy(qSym+symIx, qSeq->dna + b->qStart - qOffset, blockSize);
    memcpy(tSym+symIx, tSeq->dna + b->tStart - tOffset, blockSize);
    symIx += blockSize;
    a = b;
    }
assert(symIx == symCount);

/* Fill in score and return. */
axt->score = axtScoreDnaDefault(axt);
return axt;
}

Here is the call graph for this function:

Here is the caller graph for this function:

struct axt* chainToAxt	(	struct chain *	chain,
		struct dnaSeq *	qSeq,
		int	qOffset,
		struct dnaSeq *	tSeq,
		int	tOffset,
		int	maxGap,
		int	maxChain
	)			[read]

Definition at line 91 of file chainToAxt.c.

References axtFromBlocks(), chain::blockList, cBlock::next, cBlock::qStart, slAddHead, slReverse(), and cBlock::tStart.

{
struct cBlock *startB = chain->blockList, *a = NULL, *b;
struct axt *axtList = NULL, *axt;

for (b = chain->blockList; b != NULL; b = b->next)
    {
    if (a != NULL)
        {
        int dq = b->qStart - a->qEnd;
        int dt = b->tStart - a->tEnd;
        if ((dq > 0 && dt > 0) || dt > maxGap || dq > maxGap || (b->tEnd - startB->tStart) > maxChain)
            {
            axt = axtFromBlocks(chain, startB, b, qSeq, qOffset, tSeq, tOffset);
            slAddHead(&axtList, axt);
            startB = b;
            }
        }
    a = b;
    }
axt = axtFromBlocks(chain, startB, NULL, qSeq, qOffset, tSeq, tOffset);
slAddHead(&axtList, axt);
slReverse(&axtList);
return axtList;
}

Here is the call graph for this function:

---

Variable Documentation

char const rcsid[] = "$Id: chainToAxt.c,v 1.4 2005/04/10 14:41:21 markd Exp $" [static]

Definition at line 10 of file chainToAxt.c.

---