# 18 December 2010
# Trying v4.0.16-dev.interval

FP=/home/grobertson/FindPeaks/Code/v4.0.16-dev.interval/fp4/FindPeaks.jar
SEPR=/home/grobertson/FindPeaks/Code/v4.0.16-dev.interval/fp4/SeparateReads.jar
SRTF=/home/grobertson/FindPeaks/Code/v4.0.16-dev.interval/fp4/SortFiles.jar

WTBW=/home/grobertson/bin/x86_64/wigToBigWig

# have to run this on xhost06, can't find java16 on xhost08 or 09
JAVA=/usr/local/bin/java16


BEDDIR=/projects/remc_bigdata/Brad/3T3L1_t1_GSM535743_H3K4me1
BED=3T3L1_t1_GSM535743_H3K4me1_ChIP-seq.mm9.bed
NAME=3T3L1_t1_GSM535743_H3K4me1_ChIP-seq.mm9_ht3_q10
FPWIG=3T3L1_t1_GSM535743_H3K4me1_ChIP-seq.mm9_ht3_q10_triangle_standard

#BEDDIR=/projects/remc_bigdata/Brad/Liverf_H3K4me1
#BED=Liverf_H3K4me1_ChIP-seq_mm9.bed
#NAME=Liverf_H3K4me1_mm9_ht3_q10
#FPWIG=Liverf_H3K4me1_mm9_ht3_q10_triangle_standard

#BEDDIR=/projects/remc_bigdata/Brad/Mefs_H3K4me1
#BED=Mefs_H3K4me1_ChIP-seq_mm9.bed
#NAME=Mefs_H3K4me1_ChIP-seq_mm9_ht3_q10
#FPWIG=Mefs_H3K4me1_ChIP-seq_mm9_ht3_q10_triangle_standard


# 'clean' has chrM, not chrMT; no header lines
CLEAN=$FPWIG.clean
CHRM=$FPWIG.chrM.wig
BWFILE=$FPWIG.bigwig


# http://sourceforge.net/apps/mediawiki/vancouvershortr/index.php?title=WorkFlows#Work_Flow_-_BED
#you@machine:/path$ java -Xmx2G -jar /path/to/fp4/SeparateReads.jar bed ./input.bed out/
#you@machine:/path$ java -Xmx2G -jar /path/to/jars/fp4/SortFiles.jar bed out/ out/*.part.bed.gz
#you@machine:/path$ mkdir wig
#you@machine:/path$ java -Xmx2G -jar /home/fejes/workspace/Vancouver_Package/jars/fp4/FindPeaks.jar \
#-input out/*.part.bed.gz -aligner bed -output wig/ -dist_type 1 -name example -one_per

echo "separate reads"
$JAVA -Xmx6G -jar $SEPR interval $BEDDIR/$BED out/
echo "sort reads"
$JAVA -Xmx6G -jar $SRTF interval out/ out/*.part.bed.gz
#mkdir $OUTDIR/wig
echo "findpeaks in 'interval' mode"
echo "$JAVA -Xmx6G -jar $FP -input out/*.part.bed.gz -aligner interval -output wig/ -dist_type 1 200 -name $NAME -qualityfilter 10 -duplicatefilter -minimum 3"
$JAVA -Xmx6G -jar $FP -input out/*.part.bed.gz -aligner interval -output wig/ -dist_type 1 200 -name $NAME -qualityfilter 10 -duplicatefilter -minimum 3


#time java -Xmx2G -jar FindPeaks.jar -name test -input filename.map -output /output_path/ -aligner maq -one_per -dist_type 1 200
#echo "time $JAVA -Xmx6G -jar $FP -name $NAME -input $BAMDIR/$BAMFILE -output $OUTDIR/ -aligner sam -dist_type 1 200 -duplicatefilter -minimum 15 -prepend chr -qualityfilter 10"
#time $JAVA -Xmx6G -jar $FP -name $NAME -input $BAMDIR/$BAMFILE -output $OUTDIR/ -aligner sam -dist_type 1 200 -duplicatefilter -minimum 15 -prepend chr -qualityfilter 10

#=============================================================================
#echo "changing to 'chrMT' from 'chrM', removing header lines..."
#less $OUTDIR/$FPWIG.gz | sed 's/chrom=MT/chrom=chrM/' | sed '1,3d' > $OUTDIR/$CLEAN
#echo "less $OUTDIR/$FPWIG.gz | sed 's/chrom=MT/chrom=chrM/' > $OUTDIR/$CHRM"
#less $OUTDIR/$FPWIG.gz | sed 's/chrom=MT/chrom=chrM/' > $OUTDIR/$CHRM
#echo "gzipping"
#gzip $OUTDIR/$CHRM

#=============================================================================
WIGDIR=/projects/remc_bigdata/Brad/3T3L1_t1_GSM535743_H3K4me1/wig
WIGNAME=3T3L1_t1_GSM535743_H3K4me1_ChIP-seq.mm9_ht3_q10_triangle_standard

#CHROMINFO=/projects/mapp_etc/Mappability/fastq/mm9/mm9.chromInfo.txt
CHROMINFO=/home/grobertson/bin/x86_64/mm9.chrom.sizes
#CHROMINFO=/home/grobertson/bin/x86_64/mm8.chrom.sizes
#GENOMEFA=/projects/mapp_etc/Mappability/fastq/mm9/36mer/bwa_ind/mm9_build37_mouse.fasta
# 'chr1' is '1'
#GENOMEFA=/projects/wtsspipeline/resources/Mus_musculus/bfa_NCBI-37/all_mouse_no_chr.fasta
GENOMEFA=/projects/wtsspipeline/resources/Mus_musculus/bfa_NCBI-37/all_mouse.fasta
echo "sed '1,3d'..."
#less $WIGDIR/$WIGNAME.wig.gz | sed '1,3d' > $WIGDIR/$WIGNAME.sed13.wig
#echo "gzipping..."
#gzip $WIGDIR/$WIGNAME.sed13.wig

#echo "done"; echo "running wigToBigWig... (~/bin/x86_64/wigToBigWig)"
#echo "$WTBW -clip $WIGDIR/$WIGNAME.sed13.wig.gz $CHROMINFO $OUTDIR/$BWFILE" $WIGDIR/$WIGNAME.bigWig"
#$WTBW -clip $WIGDIR/$WIGNAME.sed13.wig.gz $CHROMINFO $OUTDIR/$BWFILE" $WIGDIR/$WIGNAME.bigWig