INNER CODE UNIT · Python

limit

Gabaldonlab/redundans · redundans.py:147

        limit = int(readLimit * fastaSize)
        if verbose:
            log.write(" Aligning %s mates per library...\n" % limit)
    return limit


def generate_gfa_file(file, threads="4"):
    """Return a gfa file generated using gfastats"""

    # First process the filename to add the new extension

    abs_path_gfa = os.path.abspath(re.sub(".fa", ".gfa", file))

    with open(abs_path_gfa, "wb") as f:

        args1 = ["gfastats", "-f", file, "--threads", str(threads), "-o", "gfa"]
        proc1 = subprocess.Popen(args1, stderr=subprocess.DEVNULL, stdout=f)
        proc1.communicate()

View source record →