def parse_discordant_by_qname(infh): frags = [] for line in infh: fields = line.strip().split('\t') frag = DiscordantFragment.from_list(fields) qname = frag.qname if len(frags) > 0 and (qname != frags[-1].qname): yield frags frags = [] frags.append(frag) if len(frags) > 0: yield frags
def parse_discordant_reads(infh): for line in infh: fields = line.strip().split('\t') frag = DiscordantFragment.from_list(fields) yield frag