Debuging read indexer

This commit is contained in:
2012-06-26 21:00:57 +00:00
parent 8e193d3206
commit a2ecbace5e

View File

@ -543,11 +543,11 @@ cdef class ReadIndex:
start1=self._buffer+ i * self._seqsize
start2=start1 + self._seqsize / 2
if (GET_SEQUSED(start1[0])==1 or GET_SEQUSED(start2[0])==1) \
and not (GET_SEQUSED(start1[0])==1 or GET_SEQUSED(start2[0])==1):
and not (GET_SEQUSED(start1[0])==1 and GET_SEQUSED(start2[0])==1):
if GET_SEQUSED(start1[0])==1:
yield self.getSeqPairAt(start1,False)
yield self.getSeqAt(start1,False)
else:
yield self.getSeqPairAt(start2,False)
yield self.getSeqAt(start2,False)
def iterreads(self,bytes word):
cdef obinuc nword[1024]