https://github.com/arjunrajlaboratory/RajLabSeqTools
Revision 2311bee1010fb2075a1866ce75c1b5956ec3859c authored by Sydney Shaffer on 11 November 2014, 19:43:50 UTC, committed by Sydney Shaffer on 11 November 2014, 19:43:50 UTC
1 parent 6313ee2
Raw File
Tip revision: 2311bee1010fb2075a1866ce75c1b5956ec3859c authored by Sydney Shaffer on 11 November 2014, 19:43:50 UTC
I forgot a variable here - added it back
Tip revision: 2311bee
unzipAndConcatenateZippedFastq.sh
#!/bin/bash

ZIPFILEDIRECTORY=$1

OUTFASTADIRECTORY=$2

for dirname in $ZIPFIELDIRECTORY/* ; do
    cd $dirname

    INPUT=`ls *001_R1*`
    SAMPLE=${INPUT%%_*}

    if [ ! -d $OUTFASTADIRECTORY/raw ]; then
        mkdir $OUTFASTADIRECTORY/raw
    fi

    if [ ! -d $OUTFASTADIRECTORY/raw/$SAMPLE ]; then
        mkdir $OUTFASTADIRECTORY/raw/$SAMPLE
    fi

    for i in *.gz; do
        gunzip -c $i > ${i%.*}
    done

    FASTQR1=${SAMPLE}_R1.fastq
    FASTQR2=${SAMPLE}_R2.fastq

    cat *R1*fastq > $OUTFASTADIRECTORY/raw/$SAMPLE/$FASTQR1
    cat *R2*fastq > $OUTFASTADIRECTORY/raw/$SAMPLE/$FASTQR2

done

back to top