INNER CODE UNIT · Python

start

steineggerlab/Metabuli · util/simulate_faa_reads.py:74

            start = rng.randint(0, len(seq) - args.read_len)
            read = seq[start:start + args.read_len]
            out.write(f">sim_{read_id} source={path.name} seq={name} start={start} len={args.read_len}\n")
            for line in wrap(read):
                out.write(line + "\n")

    sys.stderr.write(
        f"Wrote {args.reads} reads from {len(records)} eligible sequences "
        f"across {len(selected_files)} files to {args.output}\n"
    )
    if skipped_short:
        sys.stderr.write(f"Skipped {skipped_short} sequences shorter than {args.read_len} aa\n")


if __name__ == "__main__":
    main()

View source record →