Ejemplo n.º 1
0
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
Ejemplo n.º 2
0
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
Ejemplo n.º 3
0
def parse_discordant_reads(infh):
    for line in infh:
        fields = line.strip().split('\t')
        frag = DiscordantFragment.from_list(fields)
        yield frag