INNER CODE UNIT · Python

stdfrac

Gabaldonlab/redundans · redundans.py:124

        stdfrac = isstd / ismean
        # capture large stdev
        if stdfrac > stdfracTh:
            log.write(
                "[WARNING] Highly variable insert size (%.f +- %.2f) in %s - %s!\n"
                % (ismean, isstd, fq1, fq2)
            )
        # SSSPACE accepts stdfrac 0-1.0
        if stdfrac > 1:
            stdfrac = 1.0
        libraries[-1][5].append(stdfrac)
        # update counter
        i += 1
    return libraries


def get_read_limit(fasta, readLimit, verbose, log=sys.stderr):
    """Return read limit and libraries."""

View source record →