nathanae@quasar-41:/scratch/nathanae/axolotl
$ export assembly=/scratch/nathanae/axolotl/AxoTC2_trinity/Trinity.fasta
$ export oracleset=/tier2/deweylab/data/tanaka_lab/Am_asm_ETs4_RVs.fasta

$ dl-fasta-len < $assembly | python -c "import numpy as np, sys; print np.mean([float(l) for l in sys.stdin.readlines()])"
576.741842532

$ dl-fasta-len < $oracleset | python -c "import numpy as np, sys; print np.mean([float(l) for l in sys.stdin.readlines()])"
623.459056397

$ cat $oracleset | grep '^>' | wc -l
113925

$ cat $assembly | grep '^>' | wc -l
98560