--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/seqlens.py	Fri Feb 26 04:26:20 2021 +0000
@@ -0,0 +1,21 @@
+#!/usr/bin/env python
+from collections import defaultdict
+from operator import itemgetter
+from sys import stdin
+
+
+seqlens = defaultdict(int)
+next_line_seq = False
+count = 0
+for line in stdin:
+    count += 1
+    if line.startswith("@"):
+        count = 0
+        next_line_seq = True
+
+        if next_line_seq and count == 1:
+            next_line_seq = False
+            seqlens[len(line)] += 1
+
+for (length, count) in sorted(seqlens.items(), key=itemgetter(0), reverse=True):
+    print("%d\t%d" % (length, count))