# HG changeset patch
# User Steve Losh <steve@stevelosh.com>
# Date 1582502369 18000
# Node ID dfec45c7e50024b7ee0ee3926f7c5feb0e7bf1f3
# Parent  2039f25256fedf25fca1657627fc24dae4a74ade
Update

diff -r 2039f25256fe -r dfec45c7e500 README.markdown
--- a/README.markdown	Sun Feb 23 17:41:04 2020 -0500
+++ b/README.markdown	Sun Feb 23 18:59:29 2020 -0500
@@ -935,3 +935,49 @@
 just… did nothing at all.  Looks like there's a *second* Hershey Text plugin,
 this time under the AxiDraw extension menu.  I guess I'm supposed to use that
 one now?  Who knows?
+
+Played around with how to make a dot plot in gnuplot while the AxiDraw drew
+Rochester.  First I grabbed the read counts for our various FASTQ files, pulled
+from the FastQC output files:
+
+    ffind _data.txt \
+        | xargs grep 'Total Sequences' \
+        | sed -Ee 's|./([^/]+)_fastqc/.*\t|\1 |' \
+        | sort --stable -k1,1 -r \
+        | sort -k2,2 -n \
+        > seqs.txt
+
+    cat seqs.txt | head -n2
+    C3_1 21081615
+    C3_2 21081615
+
+Then I can plot a dot plot with something like:
+
+    # https://twitter.com/stevelosh/status/1231434196213796864
+    set termoption font "Routed Gothic,12"
+
+    # I wonder if there's a way to capitalize stylistically
+    set title "READ COUNTS OF INDIVIDUAL FASTQ FILES"
+    set xlabel "READS (MILLIONS)"
+
+
+    # major x tics every 2 million, with 2 minor divisions per major (i.e. minor tics are every 1 million)
+    set xtics 2
+    set mxtics 2
+
+    set grid y nox
+
+    # left justify
+    set ytics left offset -2.2, 0.0
+
+    unset key
+
+    #                                     scale to millions by hand
+    #                                                       read FASTQ names from file
+    plot [20:35][-1:26] "seqs.txt" using ($2 / 1000000.0):0:ytic(1) lc rgb 'black'
+
+And I get something like:
+
+![plot](https://i.imgur.com/YAPXHaQ.png)
+
+Neat!