def process(sequence, sample_id=None, group_id=None, source=None,
            threshold=None):
    exons = sambamba.depth_output(sequence)
    transcripts = groupby_tx(exons, sambamba=True)
    raw_stats = ((tx_id, tx_stat(tx_id, exons, threshold=threshold))
                 for tx_id, exons in iteritems(transcripts))

    if sample_id is None:
        sample_id = next(itervalues(transcripts))[0]['sampleName']
    sample_obj = Sample(id=sample_id, group_id=group_id, source=source)

    models = (make_model(sample_obj, tx_id, raw_stat) for tx_id, raw_stat
              in raw_stats)
    return Result(models=models, count=len(transcripts), sample=sample_obj)
update separater used here
dataset/ETHPy150Open robinandeer/chanjo/chanjo/load/txload.py/process--crossover with--dataset/ETHPy150Open mission-liao/pyswagger/pyswagger/utils.py/ScopeDict.sep