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

new adapater set for trimming

parent 70f09c08
No related branches found
No related tags found
No related merge requests found
...@@ -71,7 +71,7 @@ export -f dge_fastqc ...@@ -71,7 +71,7 @@ export -f dge_fastqc
dge_cutadapt(){ dge_cutadapt(){
export Ill_ADAPTERS=/projects/bioinfo/common/illumina_universal_index.fasta export Ill_ADAPTERS=/projects/bioinfo/common/cutadapt_patterns.fasta
for fastqFile in $* ; do for fastqFile in $* ; do
# DEBUG fastqFile=/projects/bioinfo/holger/projects/helin/mouse/treps_pooled/mouse_liver_polar_stage1_rep4.fastq.gz # DEBUG fastqFile=/projects/bioinfo/holger/projects/helin/mouse/treps_pooled/mouse_liver_polar_stage1_rep4.fastq.gz
...@@ -288,7 +288,7 @@ tmpDbDir=$(mktemp -d) ...@@ -288,7 +288,7 @@ tmpDbDir=$(mktemp -d)
cp -r . $tmpDbDir cp -r . $tmpDbDir
## todo remove this hack ## todo remove this hack
genome=$(echo $gtfFile | cut -f7 -d'/'); echo "genome is $genome" genome=$(echo $gtfFile | cut -f8 -d'/' | tr '[:upper:]' '[:lower:]'); echo "genome is $genome"
## make sure to use temp-r to avoid file locking problems ## make sure to use temp-r to avoid file locking problems
export R_LIBS=/tmp/r_index export R_LIBS=/tmp/r_index
......
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