#!/bin/bash

PATH=/data/results/tools/align/cufflinks-2.2.1.Linux_x86_64:$PATH; export PATH

LABELS[1]="GKR14_0h"
LABELS[2]="GKR20_0h"
LABELS[3]="GKR26_0h"
LABELS[4]="GKR15_1h"
LABELS[5]="GKR21_1h"
LABELS[6]="GKR27_1h"
LABELS[7]="GKR16_3h"
LABELS[8]="GKR22_3h"
LABELS[9]="GKR28_3h"
LABELS[10]="GKR17_6h"
LABELS[11]="GKR23_6h"
LABELS[12]="GKR29_6h"
LABELS[13]="GKR18_24h"
LABELS[14]="GKR24_24h"
LABELS[15]="GKR30_24h"
LABELS[16]="GKR19_7d"
LABELS[17]="GKR25_7d"
LABELS[18]="GKR31_7d"

#ls G*/accepted_hits/transcripts.gtf > assemblies.txt
#/data/results/tools/align/cufflinks-2.2.1.Linux_x86_64/cuffmerge   -p 8 -s /data/results/reference/mmu/Mus_musculus/UCSC/mm9/Sequence/WholeGenomeFasta/genome.fa -o merged_asm -g /data/results/reference/mmu/Mus_musculus/UCSC/mm9/Annotation/Archives/archive-2013-03-06-15-01-24/Genes/genes.gtf  assemblies.txt



#BASE_DIR=/home/dimopoulos/TopHat_test/

#rm time_Cuff-mm9.txt
CONCUR=1
THREADS=5

i=1
j=2
k=3
l=7
m=8
n=9

echo ${LABELS[i]},${LABELS[j]},${LABELS[k]}
echo ${LABELS[l]},${LABELS[m]},${LABELS[n]}
cuffdiff -L ${LABELS[i]},${LABELS[l]} --no-effective-length-correction --library-type fr-secondstrand --no-update-check --num-threads $THREADS -q -o cuffdiff_${LABELS[i]}_vs_${LABELS[l]}  merged_asm/merged.gtf www/${LABELS[$i]}.bam,www/${LABELS[$j]}.bam,www/${LABELS[$k]}.bam www/${LABELS[$l]}.bam,www/${LABELS[$m]}.bam,www/${LABELS[$n]}.bam  &> cuffdiff_${LABELS[i]}_vs_${LABELS[l]}.log

i=1
j=2
k=3
l=10
m=11
n=12

echo ${LABELS[i]},${LABELS[j]},${LABELS[k]}
echo ${LABELS[l]},${LABELS[m]},${LABELS[n]}
cuffdiff -L ${LABELS[i]},${LABELS[l]} --no-effective-length-correction --library-type fr-secondstrand --no-update-check --num-threads $THREADS -q -o cuffdiff_${LABELS[i]}_vs_${LABELS[l]}  merged_asm/merged.gtf www/${LABELS[$i]}.bam,www/${LABELS[$j]}.bam,www/${LABELS[$k]}.bam www/${LABELS[$l]}.bam,www/${LABELS[$m]}.bam,www/${LABELS[$n]}.bam  &> cuffdiff_${LABELS[i]}_vs_${LABELS[l]}.log

i=1
j=2
k=3
l=13
m=14
n=15

echo ${LABELS[i]},${LABELS[j]},${LABELS[k]}
echo ${LABELS[l]},${LABELS[m]},${LABELS[n]}
cuffdiff -L ${LABELS[i]},${LABELS[l]} --no-effective-length-correction --library-type fr-secondstrand --no-update-check --num-threads $THREADS -q -o cuffdiff_${LABELS[i]}_vs_${LABELS[l]}  merged_asm/merged.gtf www/${LABELS[$i]}.bam,www/${LABELS[$j]}.bam,www/${LABELS[$k]}.bam www/${LABELS[$l]}.bam,www/${LABELS[$m]}.bam,www/${LABELS[$n]}.bam  &> cuffdiff_${LABELS[i]}_vs_${LABELS[l]}.log

i=1
j=2
k=3
l=16
m=17
n=18

echo ${LABELS[i]},${LABELS[j]},${LABELS[k]}
echo ${LABELS[l]},${LABELS[m]},${LABELS[n]}
cuffdiff -L ${LABELS[i]},${LABELS[l]} --no-effective-length-correction --library-type fr-secondstrand --no-update-check --num-threads $THREADS -q -o cuffdiff_${LABELS[i]}_vs_${LABELS[l]}  merged_asm/merged.gtf www/${LABELS[$i]}.bam,www/${LABELS[$j]}.bam,www/${LABELS[$k]}.bam www/${LABELS[$l]}.bam,www/${LABELS[$m]}.bam,www/${LABELS[$n]}.bam  &> cuffdiff_${LABELS[i]}_vs_${LABELS[l]}.log

exit


i=1
j=2
k=3
l=4
m=5
n=6

echo ${LABELS[i]},${LABELS[j]},${LABELS[k]}
echo ${LABELS[l]},${LABELS[m]},${LABELS[n]}
cuffdiff -L ${LABELS[i]},${LABELS[l]} --no-effective-length-correction --library-type fr-secondstrand --no-update-check --num-threads $THREADS -q -o cuffdiff_${LABELS[i]}_vs_${LABELS[l]}  merged_asm/merged.gtf www/${LABELS[$i]}.bam,www/${LABELS[$j]}.bam,www/${LABELS[$k]}.bam www/${LABELS[$l]}.bam,www/${LABELS[$m]}.bam,www/${LABELS[$n]}.bam  &> cuffdiff_${LABELS[i]}_vs_${LABELS[l]}.log
exit