src/hg/inc/bamFile.h 1.3
1.3 2009/08/03 22:14:39 angie
Oops, silly flag-testing bug (& comment clarification) in bamGetQuerySequence.
Index: src/hg/inc/bamFile.h
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/inc/bamFile.h,v
retrieving revision 1.2
retrieving revision 1.3
diff -b -B -U 4 -r1.2 -r1.3
--- src/hg/inc/bamFile.h 3 Aug 2009 22:00:24 -0000 1.2
+++ src/hg/inc/bamFile.h 3 Aug 2009 22:14:39 -0000 1.3
@@ -32,9 +32,11 @@
return n;
}
char *bamGetQuerySequence(const bam1_t *bam);
-/* Return the nucleotide sequence encoded in bam. */
+/* Return the nucleotide sequence encoded in bam. The BAM format
+ * reverse-complements query sequence when the alignment is on the - strand,
+ * so here we rev-comp it back to restore the original query sequence. */
char *bamGetCigar(const bam1_t *bam);
/* Return a BAM-enhanced CIGAR string, decoded from the packed encoding in bam. */