Skip to main content
https://github.com/arjunrajlaboratory/RajLabSeqTools
10 December 2020, 08:13:14 UTC
Revision d64c7115b852e3a269ed9e3c069a3485b34bbea5 authored by Eric Sanford on 10 October 2019, 18:04:36 UTC, committed by Eric Sanford on 10 October 2019, 18:04:36 UTC
added location of hg19 reference files as comments to setEnvironmentVariables.sh script
1 parent c5b88eb
Raw File Download Save again
Tip revision: d64c7115b852e3a269ed9e3c069a3485b34bbea5 authored by Eric Sanford on 10 October 2019, 18:04:36 UTC
added location of hg19 reference files as comments to setEnvironmentVariables.sh script
Tip revision: d64c711
runHTSeq.sh 
 #!
# run from within "repo" directory

EXPERIMENT=$1
SAMPLEID=$2
gtfFile=$3

commandNAME=runHTSeq

if [ ! -d $EXPERIMENT/analyzed/$SAMPLEID/log ]; then
    mkdir $EXPERIMENT/analyzed/$SAMPLEID/log
fi

JOURNAL=$EXPERIMENT/analyzed/$SAMPLEID/log/$(date +%Y-%m-%d_%H-%M).$commandNAME.log


if [ ! -d $EXPERIMENT/analyzed/$SAMPLEID/htseq ]; then
    mkdir $EXPERIMENT/analyzed/$SAMPLEID/htseq
fi

inputFile="$EXPERIMENT/analyzed/$SAMPLEID/htseq/$SAMPLEID.nameSorted.mateFixed.sam"

countsOutFile="$EXPERIMENT/analyzed/$SAMPLEID/htseq/$SAMPLEID.htseq.stdout"
logOutFile="$EXPERIMENT/analyzed/$SAMPLEID/htseq/$SAMPLEID.htseq.stderr"

htseqCommand="python -m HTSeq.scripts.count \
	--idattr=gene_id \
	--mode=union \
	--stranded=no \
	--type=exon \
	--order=name \
	$inputFile \
	$gtfFile \
	1> $countsOutFile \
	2> $logOutFile"

echo "Starting..." >> $JOURNAL
date >> $JOURNAL
echo "$htseqCommand" >> $JOURNAL
eval "$htseqCommand"

date >> $JOURNAL
echo "Done" >> $JOURNAL