Mercurial > repos > petr-novak > re_utils
comparison RM_custom_search.py @ 11:16150c85fb3a draft
Uploaded
author | petr-novak |
---|---|
date | Tue, 04 Feb 2020 07:39:32 -0500 |
parents | a4cd8608ef6b |
children | 62fefa284036 |
comparison
equal
deleted
inserted
replaced
10:768883847008 | 11:16150c85fb3a |
---|---|
17 try: | 17 try: |
18 z=zipfile.ZipFile(f) | 18 z=zipfile.ZipFile(f) |
19 # extract only dirCLXXXX/reads.fas | 19 # extract only dirCLXXXX/reads.fas |
20 seq_list = [] | 20 seq_list = [] |
21 for filein in z.namelist(): | 21 for filein in z.namelist(): |
22 if filein.lower().startswith("seqclust/clustering/clusters/dir_cl") and filein.endswith("reads.fas"): | 22 c1 = filein.lower().startswith("seqclust/clustering/clusters/dir_cl") |
23 c2 = filein.endswith("reads.fas") | |
24 c3 = filein.endswith("reads.fasta") # in newer RE2 versions | |
25 if c1 and (c2 or c3): | |
23 outdir = filein.split("/")[3] | 26 outdir = filein.split("/")[3] |
24 outfile = outdir +"/reads.fas" | 27 outfile = outdir +"/reads.fas" |
25 source = z.open(filein) | 28 source = z.open(filein) |
26 os.mkdir(outdir) | 29 os.mkdir(outdir) |
27 target = open(outfile, "wb") | 30 target = open(outfile, "wb") |