obitools/obitools
2
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

test a new obitools script

Browse Source
This commit is contained in:
Eric Coissac
2012-06-11 17:38:26 +00:00
parent d5747e89e5
commit c20b9d33a9
1 changed files with 4 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
src/extractreads.py
View File

@ -179,7 +179,7 @@ if __name__ == '__main__':
writer = sequenceWriterGenerator(options)
nbseq=0
for seq in sequences:
for s in seq:
@ -194,5 +194,7 @@ if __name__ == '__main__':
for s in seq:
words,wordlist=addSeq2Pattern(s, wordlist, options)
st = lookforseq(s, word2seq, seq2word, options)
print >>sys.stderr,"\nWrote %d new sequences %s" % (len(st),[x.id for x in st])
nbseq+=len(st)
ldict=len(wordlist)
print >>sys.stderr,"\nWrote %d/%d [wordlist=%d] new sequences (%s)" % (len(st),nbseq,ldict,",".join([x.id for x in st]))
writer(st)