def predict_json(json_in, json_out, strategy = 'Houdayer'):
'''
predict the splicing effect
'''
print '***** PREDICTING *****'
db1 = json.loads(open(json_in, 'r').read())
records = list()
for record in db1:
p = SplicePredict(record)
p.set_ref_seq_gene(REFSEQGENE)
p.set_ref_seq(REFSEQ)
p.strategy = strategy
p.predict()
records.append(p)
with open(json_out, 'w') as f:
f.write(json.dumps(records, indent=4, ensure_ascii=False))
