INNER CODE UNIT · Python
run
MariaNattestad/assemblytics · assemblytics/index.py:9
def run(args):
coords = args.coords
output_dir = args.out
if output_dir:
os.makedirs(output_dir, exist_ok=True)
f = open(coords)
f.readline() # ignore header
fields_by_query = {}
existing_query_names = set()
existing_reference_names = set()
reference_lengths = []
query_lengths = {}
for line in f:
fields = line.strip().split(",")
query_name = fields[7]