3846f517009c43abc65d227a4695645c9b5f3e8a braney Fri Feb 15 18:31:21 2013 -0800 changes necessary to support assembly hubs (#8072) diff --git src/hg/hgc/bamClick.c src/hg/hgc/bamClick.c index 0430243..37bb181 100644 --- src/hg/hgc/bamClick.c +++ src/hg/hgc/bamClick.c @@ -39,33 +39,31 @@ if (!skipQualityScore) printf("Alignment Quality: %d
\n", core->qual); printf("CIGAR string: %s (", bamGetCigar(bam)); bamShowCigarEnglish(bam); printf(")
\n"); printf("Tags:"); bamShowTags(bam); puts("
"); printf("Flags: 0x%02x:
\n   ", core->flag); bamShowFlagsEnglish(bam); puts("
"); if (bamIsRc(bam)) printf("Note: although the read was mapped to the reverse strand of the genome, " "the sequence and CIGAR in BAM are relative to the forward strand.
\n"); puts("
"); -char nibName[HDB_MAX_PATH_STRING]; -hNibForChrom(database, seqName, nibName); -struct dnaSeq *genoSeq = hFetchSeq(nibName, seqName, tStart, tEnd); +struct dnaSeq *genoSeq = hChromSeq(database, seqName, tStart, tEnd); char *qSeq = bamGetQuerySequence(bam, FALSE); if (isNotEmpty(qSeq) && !sameString(qSeq, "*")) { char *qSeq = NULL; struct ffAli *ffa = bamToFfAli(bam, genoSeq, tStart, useStrand, &qSeq); printf("Alignment of %s to %s:%d-%d%s:
\n", itemName, seqName, tStart+1, tEnd, (isRc ? " (reverse complemented)" : "")); ffShowSideBySide(stdout, ffa, qSeq, 0, genoSeq->dna, tStart, tLength, 0, tLength, 8, isRc, FALSE); } if (!skipQualityScore && core->l_qseq > 0) { printf("Sequence quality scores:
\n\n"); UBYTE *quals = bamGetQueryQuals(bam, useStrand); int i;