# HG changeset patch # User petr-novak # Date 1580819972 18000 # Node ID 16150c85fb3afb91239f1325d630890ed5ef8d14 # Parent 7688838470082d1ad930a8cf74fbb555364ef4e9 Uploaded diff -r 768883847008 -r 16150c85fb3a RM_custom_search.py --- a/RM_custom_search.py Mon Feb 03 06:44:58 2020 -0500 +++ b/RM_custom_search.py Tue Feb 04 07:39:32 2020 -0500 @@ -19,7 +19,10 @@ # extract only dirCLXXXX/reads.fas seq_list = [] for filein in z.namelist(): - if filein.lower().startswith("seqclust/clustering/clusters/dir_cl") and filein.endswith("reads.fas"): + c1 = filein.lower().startswith("seqclust/clustering/clusters/dir_cl") + c2 = filein.endswith("reads.fas") + c3 = filein.endswith("reads.fasta") # in newer RE2 versions + if c1 and (c2 or c3): outdir = filein.split("/")[3] outfile = outdir +"/reads.fas" source = z.open(filein) diff -r 768883847008 -r 16150c85fb3a RM_custom_search.xml --- a/RM_custom_search.xml Mon Feb 03 06:44:58 2020 -0500 +++ b/RM_custom_search.xml Tue Feb 04 07:39:32 2020 -0500 @@ -1,4 +1,4 @@ - + Scan clustering results using RepeatMasker against custom database of repeats