script for creating R commands for printing plots of genome coverage and SNP peaks per chromosome