# HG changeset patch # User davidvanzessen # Date 1415109802 18000 # Node ID cd56ce86c4bf07be38940f9c78238935941d1650 # Parent 18ede27dc6b00ea44896677d3f61d92e801e8ec8 Uploaded diff -r 18ede27dc6b0 -r cd56ce86c4bf wrapper.sh --- a/wrapper.sh Tue Nov 04 08:59:34 2014 -0500 +++ b/wrapper.sh Tue Nov 04 09:03:22 2014 -0500 @@ -64,14 +64,3 @@ echo "$barcode$count$file.fastq$file.fasta${file}_trimmed.fastaReport" >> $output done < output.txt echo "" >> $output - -:' -mids = "mid\tforward\treverse" -fal = "file with sequenses on single line (no multiline)" -for mid in mids: - forwardmids = len([x for x in fal if re.match("^" + mid[1], x)]) - reversemids = len([x for x in fal if re.match(mid[2] + "$", x)]) - both = len([x for x in fal if re.match("^" + mid[1] + ".*" + mid[2] + "$", x)]) - if (forwardmids + reversemids + both) > 0: - print "%s: %d - %d - %d" % (mid[0], forwardmids, reversemids, both) -'