Update working paths based on /usr/local/share/igc

countSummary

@@ -1,3 +1,3 @@
 #!/bin/bash
-Rscript /gpfs/tools/scripts/countSummary.R $@
+Rscript /usr/local/share/igc/countSummary.R $@
 

fpkmCluster

@@ -1,3 +1,3 @@
 #!/bin/bash
-Rscript /gpfs/tools/scripts/FPKMtoCluster.R $@
+Rscript /usr/local/share/igc/FPKMtoCluster.R $@
 

fpkmPCA

@@ -1,3 +1,3 @@
 #!/bin/bash
-Rscript /gpfs/tools/scripts/pcaFPKM.R $@
+Rscript /usr/local/share/igc/pcaFPKM.R $@
 

justHOMER

@@ -8,19 +8,19 @@ else
 	GENOME=$1
 	shift 1
 	#justSTAR $GENOME $DIR $@
-	for i in *.sam
-	do
-
-		paired=$(head -n 10000 $i | samtools view -S -f 0x1 - | wc -l)
-		if [[ $paired -gt 0 ]]
-		then
-			echo "Making paired end tag dir: tags/${i%Aligned.out.sam}"
-			makeTagDirectory tags/${i%Aligned.out.sam} $i -sspe -flip 2> ${i%Aligned.out.sam}.tagDirectory.log
-		else
-			echo "Making single end tag dir: tags/${i%Aligned.out.sam}"
-			makeTagDirectory tags/${i%Aligned.out.sam} $i -flip 2> ${i%Aligned.out.sam}.tagDirectory.log
-		fi
-	done 
+#	for i in *.sam
+#	do
+#
+#		paired=$(head -n 10000 $i | samtools view -S -f 0x1 - | wc -l)
+#		if [[ $paired -gt 0 ]]
+#		then
+#			echo "Making paired end tag dir: tags/${i%Aligned.out.sam}"
+#			makeTagDirectory tags/${i%Aligned.out.sam} $i -sspe -flip 2> ${i%Aligned.out.sam}.tagDirectory.log
+#		else
+#			echo "Making single end tag dir: tags/${i%Aligned.out.sam}"
+#			makeTagDirectory tags/${i%Aligned.out.sam} $i -flip 2> ${i%Aligned.out.sam}.tagDirectory.log
+#		fi
+#	done 
 	analyzeRepeats.pl rna ${GENOME} -d tags/* -raw -count exons -condenseGenes -strand + -normMatrix > rawPlus.txt
 	analyzeRepeats.pl rna ${GENOME} -d tags/* -raw -count exons -condenseGenes -strand - -normMatrix > rawMinus.txt
 	plus=$(cat rawPlus.txt | cut -f 9 | tail -n +2 | paste -sd+ | bc)
@@ -41,7 +41,7 @@ else
 	countSummary
 	fpkmCluster
 	fpkmPCA
-	cleanSAMs .
+#	cleanSAMs .
 	echo "done aligning, building tag dirs, quantifying, and making sorted indexed bams"
 	rm -rf *_tmp *_STARtmp
 fi

justSTAR

@@ -5,7 +5,7 @@ then
 	echo "Usage: runSTAR genomeDir dirWith_R1.fastq.gz [star commands]"
 	echo "Genomes available:"
 
-	genomes=$(ls -d /gpfs/genomes/*.star)
+	genomes=$(ls -d /usr/local/share/igc/genomes/*.star)
 	for i in $genomes
 	do
 		echo ${i##*/}
@@ -14,7 +14,7 @@ else
 	GENOME=$1
 	DIR=$2
 	shift 2
-	STAR --genomeDir /gpfs/genomes/$GENOME --genomeLoad Remove > /dev/null 2> /dev/null
+	STAR --genomeDir /usr/local/share/igc/genomes/$GENOME --genomeLoad Remove > /dev/null 2> /dev/null
 	# first lets combine all of the R1 and R2 reads across lanes and multiple files
 	#for L in $(seq 99)
 	#do
@@ -48,16 +48,16 @@ else
 		if [ -e $R2 ]
 		then
 			echo "Processing paired reads: $R1 and $R2 using genome $GENOME and outputing to $PREFIX"
-			STAR --readFilesCommand zcat --genomeDir /gpfs/genomes/$GENOME --runThreadN 24 --outFilterIntronMotifs RemoveNoncanonical --readFilesIn $R1 $R2 --outFileNamePrefix $PREFIX $@ --genomeLoad LoadAndKeep --outSAMstrandField intronMotif $@ > ${PREFIX}.log
+			STAR --readFilesCommand zcat --genomeDir /usr/local/share/igc/genomes/$GENOME --runThreadN 24 --outFilterIntronMotifs RemoveNoncanonical --readFilesIn $R1 $R2 --outFileNamePrefix $PREFIX $@ --genomeLoad LoadAndKeep --outSAMstrandField intronMotif $@ > ${PREFIX}.log
 		else
 			echo "Processing single-end reads: $R1 using genome $GENOME and outputing to $PREFIX"
-			STAR --readFilesCommand zcat --outFilterIntronMotifs RemoveNoncanonical  --genomeDir /gpfs/genomes/$GENOME --runThreadN 24 --readFilesIn $R1 --outFileNamePrefix $PREFIX $@ --genomeLoad LoadAndKeep --outSAMstrandField intronMotif $@ > ${PREFIX}.log
+			STAR --readFilesCommand zcat --outFilterIntronMotifs RemoveNoncanonical  --genomeDir /usr/local/share/igc/genomes/$GENOME --runThreadN 24 --readFilesIn $R1 --outFileNamePrefix $PREFIX $@ --genomeLoad LoadAndKeep --outSAMstrandField intronMotif $@ > ${PREFIX}.log
 		fi
 		end_time="$(date -u +%M)"
 		elapsed="$((10#$end_time-10#$start_time))"
 		echo "Took $elapsed minutes to align $PREFIX"
 	done
-	STAR --genomeDir /gpfs/genomes/$GENOME --genomeLoad Remove > /dev/null
+	STAR --genomeDir /usr/local/share/igc/genomes/$GENOME --genomeLoad Remove > /dev/null
 	rm Aligned.out.sam
 	rm -rf *_tmp *_STARtmp
 fi

runSTAR

@@ -7,7 +7,7 @@ then
 	echo "set 4th arg if you don't want to run sam to bam conversion at the end (e.g. if you want sams only)"
 	echo "Genomes available:"
 
-	genomes=$(ls -d /gpfs/genomes/*.star)
+	genomes=$(ls -d /usr/local/share/igc/genomes/*.star*)
 	for i in $genomes
 	do
 		echo ${i##*/}