#!/bin/sh
#$ -S /bin/sh
# sh chipatlas/sh/makePeakBrowser/splitByChr.sh chipatlas/results/sacCer3/bmap/pmd/Bed/SRX3236110.pmd.bed sacCer3 pmd SRX3236110 bmap
# sh chipatlas/sh/makePeakBrowser/splitByChr.sh chipatlas/results/mm10/Bed10/Bed/SRX999043.10.bed mm10 10 SRX999043 chip

if [ $5 = "chip" ]; then
  mkdir -p tmpDirForPeakBrowser/splitByChr/$2/$3/$JOB_ID
else
  mkdir -p tmpDirForPeakBrowser/splitByChr/$2/bmap/$JOB_ID
fi

cat $1| awk -v genome=$2 -v type=$3 -v JOB_ID=$JOB_ID -v SRX=$4 -v Seq=$5 -v OFS='\t' -F '\t' '{
  if (NR == 1) {
    if (Seq == "chip") fn = "tmpDirForPeakBrowser/splitByChr/"genome"/"type"/"JOB_ID"/"
    if (Seq == "bmap") fn = "tmpDirForPeakBrowser/splitByChr/"genome"/bmap/"JOB_ID"/"
  }
  if (Seq == "chip") print $1, $2, $3, SRX, $5 >> fn $1
  if (Seq == "bmap") print $0, SRX, type >> fn $1
}'