#!/bin/sh FILENAME=$(basename "$1") #sample.hist or sample.pe.fq.hist #BASE=${FILENAME/.*} BASE=${FILENAME/.hist} echo "File base: $BASE" echo "set term png" > ${BASE}.gp echo "set output \"${BASE}.png\"" >> ${BASE}.gp echo "set datafile separator \",\"" >> ${BASE}.gp echo "set logscale x" >> ${BASE}.gp echo "set logscale y" >> ${BASE}.gp echo "set grid" >> ${BASE}.gp echo "set xlabel \"kmer frequency\"" >> ${BASE}.gp echo "set ylabel \"number of distinct kmer\"" >> ${BASE}.gp echo "plot \"$1\" using 1:2 with lines" >> ${BASE}.gp gnuplot ${BASE}.gp