https://github.com/jgurtowski/ectools
Tip revision: 031eb0300c82392915d8393a5fedb4d3452b15bf authored by James Gurtowski on 01 December 2014, 14:07:00 UTC
Update LICENSE.md
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