Skip to content

Instantly share code, notes, and snippets.

@sjdv1982
Last active July 31, 2019 13:25
Show Gist options
  • Save sjdv1982/d91fef1225263fbd1a8884c5f6f2f186 to your computer and use it in GitHub Desktop.
Save sjdv1982/d91fef1225263fbd1a8884c5f6f2f186 to your computer and use it in GitHub Desktop.
import numpy as np
fragments = np.load("fragments_clust.npy")
redun = np.load("redundancy-masks/1.0-seqclust-50.npy").astype(bool)
fragments = fragments[redun]
fragments_aa = np.array([f for f in fragments if f["motif"] == b"AAA"])
clus = [f for f in fragments_aa if f["clust1.0"] == 700]
print(clus)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment