Advertisement
CAROJASQ

Untitled

Apr 27th, 2017
95
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 0.46 KB | None | 0 0
  1. def extract_superkmers(reads_matrix, minimizer_matrix):
  2.     superkmers = [] # Composed of tuples (superkmer, minimizer)
  3.     for read_number, row in minimizer_matrix:
  4.         for v in row:
  5.             minimizer = v & 0b11111111111100000000000000000000
  6.             start = v & 0b00000000000000000000001111111111
  7.             end = v & 0b00000000000011111111110000000000
  8.             superkmers.append((reads_matrix[read_number][start:end], minimizer)
  9.     return superkmers
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement