#!/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