ce07fe9f6695acd811e90fbd646830ca7858b525 markd Mon Sep 28 13:41:50 2015 -0700 fixed various problems detected by making pslCheck more stringent diff --git src/inc/psl.h src/inc/psl.h index ef54630..5d700e4 100644 --- src/inc/psl.h +++ src/inc/psl.h @@ -282,30 +282,33 @@ char *tName, bioSeq *tSeq, int tStart, int tEnd, int exnStarts[], int exnEnds[], int exnCnt, FILE *f); /* Show protein/DNA alignment or translated DNA alignment in HTML format. */ struct psl* pslNew(char *qName, unsigned qSize, int qStart, int qEnd, char *tName, unsigned tSize, int tStart, int tEnd, char *strand, unsigned blockSpace, unsigned opts); /* create a new psl with space for the specified number of blocks allocated. * pslGrow maybe used to expand this space if needed. Valid options are * PSL_XA_FORMAT. */ void pslGrow(struct psl *psl, int *blockSpacePtr); /* Increase memory allocated to a psl to hold more blocks. blockSpacePtr * should point the the current maximum number of blocks and will be * updated to with the new amount of space. */ +void pslComputeInsertCounts(struct psl *psl); +/* compute numInsert and baseInsert fields from the blocks */ + struct psl* pslFromGff3Cigar(char *qName, int qSize, int qStart, int qEnd, char *tName, int tSize, int tStart, int tEnd, char* strand, char *cigar); /* create a PSL from a GFF3-style cigar formatted alignment */ int pslRangeTreeOverlap(struct psl *psl, struct rbTree *rangeTree); /* Return amount that psl overlaps (on target side) with rangeTree. */ float pslIdent(struct psl *psl); /* computer fraction identity */ float pslQueryAligned(struct psl *psl); /* compute fraction of query that was aligned */ INLINE unsigned pslQStart(struct psl *psl, int blkIdx)