Script <myIllu_02_jellyfish.sh>

#!/bin/sh

export LD_LIBRARY_PATH=/usr/local/apps/jellyfish-2.2.3/lib
KMER=21
THREADS=20

BASE=${1%%.*}
echo "File name base: $BASE"

echo "Running jellyfish..."
jellyfish count -s 1G -m ${KMER} -t ${THREADS} -o kmer -C $1
jellyfish dump -c -t kmer > ${BASE}.kmer${KMER}.txt

echo "Preparing histogram..."
awk '{print $2}' ${BASE}.kmer${KMER}.txt | sort -n | uniq -c | awk '{print $2 "," $1}' > ${BASE}.jf.hist