https://github.com/jgurtowski/ectools
Raw File
Tip revision: 031eb0300c82392915d8393a5fedb4d3452b15bf authored by James Gurtowski on 01 December 2014, 14:07:00 UTC
Update LICENSE.md
Tip revision: 031eb03
kmer.py
#!/usr/bin/env python

import sys

from itertools import imap

from seqio import iteratorFromExtension
from nucio import fileIterator


##Create Kmers

if not len(sys.argv) == 3:
    sys.exit("Usage: kmer.py k-size in.fa\n")

fn = sys.argv[2]
ksize = int(sys.argv[1])

for record in fileIterator(fn, iteratorFromExtension(fn)):
    seq = record.seq
    starts = range(len(seq)-ksize+1)
    kmers = imap(lambda start: seq[start:start+ksize], starts)
    for kmer in kmers:
        print kmer
back to top