def _compose_paths(exon_dict, path_dict, number_of_ns):
ns = "N" * number_of_ns
for transcript_id, exon_list in sorted(path_dict.items()):
exon_seqs = [str(exon_dict[exon_id].seq) for exon_id in exon_list]
yield SeqRecord(
id=transcript_id,
description=",".join(exon_list),
seq=Seq(ns.join(exon_seqs))
)
评论列表
文章目录