#!/usr/bin/env bash
##Author: guochangquan
##Copyright 20210929
##company: NuProbe

set -e

function usage() {
    echo "Usage: $0 [ -s <sample> ] [ -i <fq_dir> ] [ -o <analysis_dir> ] [ -x <script_dir> ] [ -c <FeBY2_CNV_outdir> ]"
    echo "-s sample name"
    echo "-i input fastq dir"
    echo "-x script_dir"
    echo "-o analysis_dir"
    echo "-c FeBY2 CNV outdir"
}

if [ $# -eq 0  ]
then
    usage
    exit
fi

while getopts ":s:c:x:o:i:h" opt; do
    case $opt in
        i)
            fq_dir=$OPTARG
            ;;
        o)
            out_dir=$OPTARG
            ;;
        s)
            sample=$OPTARG
            ;;
        c)
            cnv_dir=$OPTARG
            ;;
        x)
            feby2_script=$OPTARG
            ;;
        h)
            usage
            exit
            ;;
        \?)
            echo "Invalid option: -$OPTARG" > /dev/stderr
            usage
            exit
            ;;
    esac
done
echo "Input FASTQ dir:" $fq_dir
echo "Output BAM dir:" $out_dir
echo "Sample Name:" $sample
echo "CNV Analysis dir:" $cnv_dir
echo $feby2_script

. $feby2_script/conf/feby2.cfg

test -d $out_dir/${sample} || mkdir ${outdir}/${sample}
cd $out_dir/${sample}
module=$feby2_script/module

### FASTQ TRIM and Mapping
test -s $sample.bam || sh $module/run_FQ2BAM_Speedseq.sh -s $sample -i $fq_dir -o $out_dir -x $feby2_script
### SNV/INDEL CAll and annotation
test -s $sample.gatk.vcf.gz || sh $module/run_BAM2VCF_HC.sh -s $sample -o $out_dir -x $feby2_script
test -s $sample.sv.vcf.gz || sh $module/run_BAM2SV_Manta.sh -s $sample -o $out_dir/${sample} -x $feby2_script
### QC
sh $module/run_QC.sh -s $sample -o $out_dir -x $feby2_script
### yesuan
sh $module/run_yesuan.sh -s $sample -o $out_dir -x $feby2_script
### prepare for CNV
sh $module/run_CNV_prep.sh -s $sample -o $out_dir -x $feby2_script -c $cnv_dir