src/blat/blat.c 1.113
1.113 2009/02/10 21:58:55 kent
Makin blat work with U's instead of T's if q=rna is set.
Index: src/blat/blat.c
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/blat/blat.c,v
retrieving revision 1.112
retrieving revision 1.113
diff -b -B -U 4 -r1.112 -r1.113
--- src/blat/blat.c 25 Nov 2007 17:17:55 -0000 1.112
+++ src/blat/blat.c 10 Feb 2009 21:58:55 -0000 1.113
@@ -306,8 +306,10 @@
Bits *qMaskBits = maskQuerySeq(seq, isProt, maskQuery, lcMask);
struct dnaSeq trimmedSeq;
ZeroVar(&trimmedSeq);
trimSeq(seq, &trimmedSeq);
+if (qType == gftRna || qType == gftRnaX)
+ memSwapChar(trimmedSeq.dna, trimmedSeq.size, 'u', 't');
searchOne(&trimmedSeq, gf, outFile, isProt, maskHash, qMaskBits);
*retTotalSize += seq->size;
*retCount += 1;
bitFree(&qMaskBits);