#!/usr/bin/env python
"""
This script needs work! Or will we never use it again?
There should be a --find option, or something.
"""
from __future__ import print_function
from dark.blast import BlastRecords
from dark.graphics import evalueGraph
if __name__ == '__main__':
import sys
if len(sys.argv) < 2:
print('Usage: %s BLAST-hitfile...' % sys.argv[0], file=sys.stderr)
else:
blastRecords = BlastRecords(sys.argv[1:])
evalueGraph(blastRecords.records(), 10, 10,
find=lambda title: title.find('HKU4') > -1 and
title.find('complete genome') > -1, titles=False)
#!/usr/bin/env python
from __future__ import print_function
from dark.blast import BlastRecords
from dark.blast.records import printBlastRecordForDerek
if __name__ == '__main__':
import sys
if len(sys.argv) < 2:
print('Usage: %s hitfiles...' % sys.argv[0], file=sys.stderr)
else:
blastRecords = BlastRecords(sys.argv[1:])
hits = 0
for i, result in enumerate(blastRecords.records(), start=1):
oldHits = hits
hits += printBlastRecordForDerek(result)
