def get_transcript_dict(gp, filter_set):
transcripts = seq_lib.get_gene_pred_transcripts(gp)
d = seq_lib.transcript_list_to_dict(transcripts, noDuplicates=True)
r = defaultdict(list)
for aln_id, rec in d.iteritems():
tx_id = strip_alignment_numbers(aln_id)
if tx_id in filter_set:
r[tx_id].append(rec)
return r
def get_transcript_dict(gp, filter_set):
transcripts = seq_lib.get_gene_pred_transcripts(gp)
d = seq_lib.transcript_list_to_dict(transcripts, noDuplicates=True)
r = defaultdict(list)
for aln_id, rec in d.iteritems():
tx_id = strip_alignment_numbers(aln_id)
if tx_id in filter_set:
r[tx_id].append(rec)
return r
def load_gp(path):
tm_recs = seq_lib.get_gene_pred_transcripts(path)
tm_dict = seq_lib.transcript_list_to_dict(tm_recs, noDuplicates=True)
return tm_dict
def load_gp(path):
tm_recs = seq_lib.get_gene_pred_transcripts(path)
tm_dict = seq_lib.transcript_list_to_dict(tm_recs, noDuplicates=True)
return tm_dict
