#!/usr/bin/env python

from Bio.GenBank import FeatureParser, Iterator
import sys, os

dir = sys.argv[1]
files = []
for file in os.listdir(dir):
    files.append(os.path.join(dir, file))
pipe = os.popen('./gbkftr ' + ' '.join(files), 'r')

for record in Iterator(pipe, FeatureParser()):
    print record, record.name, '- features:', len(record.features)
