X-Git-Url: https://git.donarmstrong.com/?p=uiuc_igb_scripts.git;a=blobdiff_plain;f=fpkm_from_srr;h=8e207f3a2af315a8a86a937776c1c5db806b418d;hp=a7ab6e0a4ef9c03c818c9737f1c0a810cc7b87f8;hb=HEAD;hpb=13f54c1d3fd7c8ea0294f0b6ae81eedfca6e1ea5

diff --git a/fpkm_from_srr b/fpkm_from_srr
index a7ab6e0..8e207f3 100755
--- a/fpkm_from_srr
+++ b/fpkm_from_srr
@@ -1,17 +1,18 @@
 #!/bin/bash
-# PBS -M id -un
 
 CORES=8
-MAX_MEM=32
+MAX_MEM="18GB"
 
+set -e;
 
 function actually_run_the_alignment {
     module load sratoolkit/2.3.5-2;
     module load bwa/0.7.10;
     module load cufflinks/2.2.1;
-    (cd $FROM_SRR_SRR;
-     fastq-dump ${FROM_SRR_SRR}.sra;
-     bwa mem -t ${CORES} $FROM_SRR_REF ${FROM_SRR_SRR}.fastq.gz | samtools view -b -o ${FROM_SRR_SRR}.bam;
+    module load samtools/1.0;
+    (cd "${FROM_SRR_SRR}";
+     fastq-dump --gzip ${FROM_SRR_SRR}.sra;
+     bwa mem -t ${CORES} $FROM_SRR_REF ${FROM_SRR_SRR}.fastq.gz | samtools view -b -o ${FROM_SRR_SRR}.bam -;
      samtools sort ${FROM_SRR_SRR}.bam -o ${FROM_SRR_SRR}.sorted.bam -T temp
      cufflinks -p ${MAX_MEM} -G $FROM_SRR_GTF ${FROM_SRR_SRR}.sorted.bam;
      rm -f ${FROM_SRR_SRR}.fastq.gz ${FROM_SRR_SRR}.bam
@@ -24,11 +25,10 @@ function set_up_the_run {
     FROM_SRR_GTF="$3"
     export FROM_SRR_SRR FROM_SRR_GTF FROM_SRR_REF;
     mkdir -p "$FROM_SRR_SRR";
-    set -e;
-    (cd $FROM_SRR_SRR;
-     rsync -avP rsync://ftp-trace.ncbi.nlm.nih.gov/sra/sra-instant/reads/ByRun/sra/SRR/${FROM_SRR_SRR:0:6}/${FROM_SRR_SRR}/${FROM_SRR_SRR}.sra;
+    (cd "$FROM_SRR_SRR";
+     rsync -avP rsync://ftp-trace.ncbi.nlm.nih.gov/sra/sra-instant/reads/ByRun/sra/SRR/${FROM_SRR_SRR:0:6}/${FROM_SRR_SRR}/${FROM_SRR_SRR}.sra .;
     );
-    qsub -q budget -v FROM_SRR_SRR,FROM_SRR_GTF,FROM_SRR_REF -S /bin/bash -d $(pwd) -M donarm@illinois.edu -l "nodes=1:ppn=${CORES},mem=${MAX_MEM}" "$0"    
+    qsub -q budget -v FROM_SRR_SRR,FROM_SRR_GTF,FROM_SRR_REF -S /bin/bash -d $(pwd) -m a -M donarm@illinois.edu -l "nodes=1:ppn=${CORES},mem=${MAX_MEM}" "$0"    
 };