INNER CODE UNIT · Python

parse_files

torognes/swarm · scripts/graph_plot.py:87

def parse_files(swarms, internal_structure, cluster, drop):
    """
    """
    # List amplicon ids and abundances
    amplicons = list()
    with open(swarms, "r") as swarms:
        for i, swarm in enumerate(swarms):
            if i == cluster - 1:
                # Deal with ";size=" in a rather clumsy way... but it works
                print("Reading target cluster", file=sys.stdout)
                amplicons = [
                    tuple(
                        item.replace(";size=", "_").rstrip(";").rsplit("_", 1))
                    for item in swarm.strip().split(" ")]
                break

    # Drop amplicons with a low abundance (remove connections too)
    if drop:

View source record →