From ee4637e5f2ce678421768064f6aead0bcae758b1 Mon Sep 17 00:00:00 2001 From: Till Hartmann Date: Tue, 26 May 2020 14:26:40 +0200 Subject: [PATCH] add dinopy version of fqcnt --- fqcnt/fqcnt_py9x_dinopy.py | 15 +++++++++++++++ 1 file changed, 15 insertions(+) create mode 100755 fqcnt/fqcnt_py9x_dinopy.py diff --git a/fqcnt/fqcnt_py9x_dinopy.py b/fqcnt/fqcnt_py9x_dinopy.py new file mode 100755 index 0000000..9932c31 --- /dev/null +++ b/fqcnt/fqcnt_py9x_dinopy.py @@ -0,0 +1,15 @@ +#!/usr/bin/env python +import sys +import dinopy + +if __name__ == "__main__": + if len(sys.argv) != 2: + print("Usage: {} ".format(sys.argv[0])) + sys.exit(0) + fqr = dinopy.FastqReader(sys.argv[1]) + n, sl, ql = 0, 0, 0 + for (seq, _, qual) in fqr.reads(quality_values=True): + n += 1 + sl += len(seq) + ql += len(qual) + print("{}\t{}\t{}".format(n, sl, ql))