dge_master_template.sh

export baseDir=<<PATH_TO_BASEDIR>>
export project=<<PROJECTNAME>>
# screen -R $project


## intially
DGE_HOME=/projects/bioinfo/holger/bioinfo_templates/dge_workflow
## when customization is needed
#export DGE_HOME=/sw/users/brandl/projects/<<PROJECTNAME>>

source $DGE_HOME/dge_utils.sh
export PATH=$DGE_HOME:$DGE_HOME/../misc/:$PATH


########################################################################################################################
### QC

dge_fastqc $(ls *fastq.gz) &


########################################################################################################################
### Trim low-quality bases and remove left-over adapters

mcdir $baseDir/trimmed

dge_cutadapt $(ls $baseDir/treps_pooled/*fastq.gz) 2>&1 | tee cutadapt.log

dge_fastqc $(ls *fastq.gz) &


########################################################################################################################
### Align the reads

mcdir $baseDir/mapped

fastqFiles=$(ls $baseDir/trimmed/*.fastq.gz)
igenome=<<<<TBD>>>>
## Example:
## igenome=/projects/bioinfo/igenomes/Canis_familiaris/Ensembl/CanFam3.1


dge_tophat_se -i $igenome $fastqFiles 2>&1 | tee mapped.log

mailme "$project: mapping done"


########################################################################################################################
### Basic Alginment QC and technical replicate grouping

mcdir $baseDir/trep_pooled


bio_reps=<<<biological replicates labels>>>
## Examples
# bio_reps=$(csvcut  -tc bio_sample  ../lanereps_pooled/renaming_scheme.txt |  tail -n+2 | sort -u | xargs echo | tr " " ",")
## bio_reps="ctrl,isnm1"

dge_merge_treps $baseDir/mapped/ $bio_reps


########################################################################################################################
### Do the differential expression analysis

mcdir $baseDir/cuffdiff

gtfFile=$igenome/Annotation/Genes/genes.gtf

## define labels to split bam files into replicate groups
#labels=$(csvcut -t -c1 $baseDir/lanereps_pooled/sample2replicates.txt | tail -n+2 | uniq | xargs echo | sed 's/ /,/g')
labels=<<<<TBD>>>>

dge_cuffdiff $gtfFile $baseDir/mapped $labels

mailme "$project: cuffdiff done"


mcdir $baseDir/cuffdiff/dge_report
#export DGE_PARAMS="-S"
spin.R $DGE_HOME/dge_analysis.R $baseDir/cuffdiff


########################################################################################################################
### Sync back to project space

# ... project specific stuff
screen -R rsync_$project

## main sync
~/bin/unison $baseDir ssh://bioinfo///home/brandl/mnt/<<MOUNT_PATH>> -fastcheck true -times -perms 0


rsync -avsn --delete  $baseDir brandl@fileserver:/projects//file/server/path

mailme "$project: rsync done"