836ef4e2bcac0c455308a6e877c9b6e707a3a436 chinhli Mon Jun 27 11:10:59 2011 -0700 redmine #4036: per data provider's request, skip printing alignment quality for NumtS item bam detail page diff --git src/hg/hgc/bamClick.c src/hg/hgc/bamClick.c index e026ff1..db809a5 100644 --- src/hg/hgc/bamClick.c +++ src/hg/hgc/bamClick.c @@ -25,30 +25,31 @@ /* Maybe make this an option someday -- for now, I find it too confusing to deal with * CIGAR that is anchored to positive strand while showing rc'd sequence. I think * to do it right, we would need to reverse the CIGAR string for display. */ static boolean useStrand = FALSE; static boolean skipQualityScore = FALSE; static void singleBamDetails(const bam1_t *bam) /* Print out the properties of this alignment. */ { const bam1_core_t *core = &bam->core; char *itemName = bam1_qname(bam); int tLength = bamGetTargetLength(bam); int tStart = core->pos, tEnd = tStart+tLength; boolean isRc = useStrand && bamIsRc(bam); printPosOnChrom(seqName, tStart, tEnd, NULL, FALSE, itemName); +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];