def _parse_preseq_logs(f):
    """ Go through log file looking for preseq output """

    lines = f['f'].splitlines()
    header = lines.pop(0)

    data_is_bases = False
    if header.startswith('TOTAL_READS	EXPECTED_DISTINCT'):
        pass
    elif header.startswith('TOTAL_BASES	EXPECTED_DISTINCT'):
        data_is_bases = True
    elif header.startswith('total_reads	distinct_reads'):
        pass
    else:
        log.debug("First line of preseq file {} did not look right".format(f['fn']))
        return None, None

    data = dict()
    for l in lines:
        s = l.split()
        # Sometimes the Expected_distinct count drops to 0, not helpful
        if float(s[1]) == 0 and float(s[0]) > 0:
            continue
        data[float(s[0])] = float(s[1])

    return data, data_is_bases
def _parse_preseq_logs ( f ) : lines = f [ 'f' ] . splitlines ( ) header = lines . pop ( 0 ) data_is_bases = False if header . startswith ( 'TOTAL_READS	EXPECTED_DISTINCT' ) : pass elif header . startswith ( 'TOTAL_BASES	EXPECTED_DISTINCT' ) : data_is_bases = True elif header . startswith ( 'total_reads	distinct_reads' ) : pass else : log . debug ( "First line of preseq file {} did not look right" . format ( f [ 'fn' ] ) ) return None , None data = dict ( ) for l in lines : s = l . split ( ) if float ( s [ 1 ] ) == 0 and float ( s [ 0 ] ) > 0 : continue data [ float ( s [ 0 ] ) ] = float ( s [ 1 ] ) return data , data_is_bases
Go through log file looking for preseq output