Advertisement
RonWeber

Untitled

Feb 22nd, 2020
161
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 1.39 KB | None | 0 0
  1. print((''.join(["{}\t{}\t{}\t{}\n".format(su["ls"][1], su["ls"][0], su["count"], su["count"] / totalnumseqs) for su in [vsu for vsu in sorted([[(lsdict.__setitem__("count", lsdict["count"] + tcval) or lsdict if lsdict["ls"][5] in tckey else lsdict) for (tckey, tcval) in functools.reduce(lambda count, parts: ((globals().__setitem__("totalnumseqs", totalnumseqs + int(parts[1])) or count.__setitem__(parts[0], count.get(parts[0], 0) + int(parts[1]))) or count), [e.split() for e in open(sys.argv[1]).readlines() + (open(sys.argv[2]).readlines())], globals().__setitem__("totalnumseqs", 0) or {}).items()][0] for lsdict in [dict(ls=ls, count=0) for ls in [a.split() for a in open(sys.argv[3]).readlines()[1:]]]], key=lambda entry: entry["count"], reverse=True) if (vsu["count"] > 0)]]) if os.path.basename(sys.argv[0]) == "mirna.py" else ''.join(["{} {}\n".format(onecount["string"], onecount["count"]) for onecount in [a for a in [dict(string=s, count=len(re.findall("(?=" + s + ")", open(sys.argv[1]).read()))) for s in [''.join(s) for s in itertools.product('ACGT', repeat=int(sys.argv[2]))]] if (a["count"] > int(sys.argv[3]))]])) if ([globals().__setitem__(lib, __import__(lib)) for lib in ["sys","functools","itertools","re","os"]] or (os.path.basename(sys.argv[0]) in ["kmer.py","mirna.py"])) else "Please name this file either kmer.py or mirna.py and invoke with the correct arguments.") #Just a simple one-liner
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement