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