Skip to content
Snippets Groups Projects
Commit 7965e302 authored by Holger Brandl's avatar Holger Brandl
Browse files

started new star analysis

parent 0edeacfc
No related branches found
No related tags found
No related merge requests found
......@@ -107,15 +107,21 @@ for fastqFile in $fastqFiles ; do
# --quantMode GeneCounts see https://groups.google.com/forum/#!searchin/rna-star/GeneCounts/rna-star/gZRJx3ElRNo/p5FjBYKuY00J
# --outSJfilterCountUniqueMin see https://groups.google.com/forum/#!topic/rna-star/_1BeAlGUmpA
mysub "${project}__star__${fastqBaseName}" "
jl submit -j .starjobs -q medium -t 5 -n "${project}__star__${fastqBaseName}" "
STAR --genomeDir $star_index --readFilesIn $fastqFile --runThreadN 6 --readFilesCommand zcat --outFileNamePrefix ${fastqBaseName}. --outSAMtype BAM SortedByCoordinate --outSAMstrandField intronMotif --sjdbGTFfile $gtfFile --outFilterIntronMotifs RemoveNoncanonicalUnannotated --outFilterType BySJout --quantMode GeneCounts --outFilterMultimapNmax 1 --outSJfilterCountUniqueMin 8 3 3 3
mv ${fastqBaseName}.Aligned.sortedByCoord.out.bam ${fastqBaseName}.bam
samtools index ${fastqBaseName}.bam
" -n 5 -R span[hosts=1] -q medium | joblist .starjobs
"
done
wait4jobs .starjobs
jl wait .starjobs
if [ -n "$(jl failed .starjobs)" ]; then
echo "some jobs failed"; exit 1
fi
jl status --report .starjobs
rm -rf *STARgenome *.Log.progress.out _STARtmp *.SJ.out.tab *Log.out
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment