#!/bin/sh
#$ -S /bin/sh
# qsub -e /dev/null -o /dev/null chipatlas/sh/makePeakBrowser/allPeaks_light.sh hg38 05
genome=$1
qval=$2

inBed=tmpDirForPeakBrowser/sortedBed/$genome/$genome.$qval.bed
outFn=ALL.ALL.$qval.AllAg.AllCell.bed
litFn=tmpDirForPeakBrowser/sortedBed/$genome/allPeaks_light.$genome.$qval.bed
linNum=tmpDirForPeakBrowser/linNum/$genome.$outFn

if [ $qval != "bs" ]; then
  cat $inBed| awk -v linNum=$linNum -v genome=$genome -v qval=$qval '{
    print
    x[$4]++
  } END {
    cmd = "sort"
    if (qval == "05") {
      fn = "chipatlas/results/"genome"/summary/allLineNum.txt"
      for (srx in x) print srx "\t" x[srx] + 0 |& cmd
      close(cmd, "to")
      while ((cmd |& getline) > 0) print > fn
      close(cmd)
    }
    print NR > linNum
  }'| split -l 50000000 - tmpDirForPeakBrowser/inSilicoChIP/$genome/$outFn"."
else
  cat $inBed| awk -F '\t' '$5 == "hypermr" {
    x[$4]++
  } END {
    for (srx in x) print srx "\t" x[srx] + 0
  }'| sort > chipatlas/results/$genome/bmap/summary/allLineNum.txt
fi

gzip -c $inBed > tmpDirForPeakBrowser/allPeaks_light/$genome/allPeaks_light.$genome.$qval.bed.gz

exit