#!/usr/bin/env python3 import pykmertools as kt def main() -> None: seq = ( "ATGCGATATCGTAGGCGTCGATGGAGAGCTAGATCGATCGATCTAAATCCCGATCGATTCCGAGCGCGATCAAAGCGCGATAGGCTAGCTAAAGCTAGCA" ) wsize = 31 msize = 7 print(f"sequence length = {len(seq)}") print(f"window size = {wsize}") print(f"minimiser size = {msize}") print() print("Minimiser segments:") mg = kt.MinimiserGenerator(seq, wsize, msize) for idx, (mmer, start, end) in enumerate(mg, start=1): mmer_str = mg.to_acgt(mmer) print( f"{idx:2d}. " f"min={mmer_str} ({mmer}) " f"range=[{start}, {end}) " f"segment={seq[start:end]}" ) if __name__ == "__main__": main()