#!/bin/bash
PATH=/data/results/tools/align/tophat-2.0.10.Linux_x86_64:/home/reczko/ocamlbrew/ocaml-4.01.0/.opam/system/bin:/home/reczko/ocamlbrew/ocaml-4.01.0/.opam/system/bin:/home/reczko/ocamlbrew/ocaml-4.01.0/bin:/data/results/tools/chipseq/macs/bin:/data/results/tools/bedtools/BEDTools-Version-2.14.3/bin:/data/results/tools/chipseq/ngsplot/bin:/data/results/tools/r/R-3.1.2/bin:/usr/local/bin:/usr/bin:/bin:/usr/local/games:/usr/games:/data/results/tools/chipseq/homer/.//bin/; export PATH
echo "Using tophat2 installation at"
which tophat2
tophat --version
PATH=/data/results/tools/align/bowtie2-2.1.0:/home/reczko/ocamlbrew/ocaml-4.01.0/.opam/system/bin:/home/reczko/ocamlbrew/ocaml-4.01.0/.opam/system/bin:/home/reczko/ocamlbrew/ocaml-4.01.0/bin:/data/results/tools/chipseq/macs/bin:/data/results/tools/bedtools/BEDTools-Version-2.14.3/bin:/data/results/tools/chipseq/ngsplot/bin:/data/results/tools/r/R-3.1.2/bin:/usr/local/bin:/usr/bin:/bin:/usr/local/games:/usr/games:/data/results/tools/chipseq/homer/.//bin/; export PATH
echo "Using bowtie installation at"
which bowtie2
bowtie2 --version
mkdir bowtie2_003; cd bowtie2_003;  bowtie2 --local --very-sensitive-local -p 8 --mm -x /data/results/reference/mmu/Mus_musculus/UCSC/mm9/Sequence/Bowtie2Index/genome -U ../R_2015_04_15_12_36_34_user_ION-171-ITlab_YKRNAseq_150415_ITR41-45.ITR41-wt.IonXpressRNA_003.fastq  | /data/results/tools/samtools/samtools-0.1.19/samtools view -uhS -F4 - | /data/results/tools/samtools/samtools-0.1.19/samtools sort - merged ; /data/results/tools/samtools/mysamtools-0.1.19/samtools sort  -@ 8 -n merged.bam sort; rm merged.bam; samtools view sort.bam   | awk -f /data/results/tools/formats/sam-remove-read-multimaps-to-same-location2.awk | samtools view -bt /data/results/reference/mmu/Mus_musculus/UCSC/mm9/Sequence/WholeGenomeFasta/genome.fa.fai - > sort_uniq.bam; cd ..;
mkdir bowtie2_011; cd bowtie2_011;  bowtie2 --local --very-sensitive-local -p 8 --mm -x /data/results/reference/mmu/Mus_musculus/UCSC/mm9/Sequence/Bowtie2Index/genome -U ../R_2015_04_15_12_36_34_user_ION-171-ITlab_YKRNAseq_150415_ITR41-45.ITR45-CEBPaKO.IonXpressRNA_011.fastq  | /data/results/tools/samtools/samtools-0.1.19/samtools view -uhS -F4 - | /data/results/tools/samtools/samtools-0.1.19/samtools sort - merged ; /data/results/tools/samtools/mysamtools-0.1.19/samtools sort  -@ 8 -n merged.bam sort; rm merged.bam; samtools view sort.bam   | awk -f /data/results/tools/formats/sam-remove-read-multimaps-to-same-location2.awk | samtools view -bt /data/results/reference/mmu/Mus_musculus/UCSC/mm9/Sequence/WholeGenomeFasta/genome.fa.fai - > sort_uniq.bam; cd ..;