Skip to content
Snippets Groups Projects
Commit 82f61996 authored by Rachel  LEGENDRE's avatar Rachel LEGENDRE
Browse files

fix some code issues

parent f99e0f4e
Branches
No related tags found
No related merge requests found
......@@ -35,7 +35,7 @@ rule kallisto_index:
log:
kallisto_index_log
threads:
config['pseudomapping']['threads']
config['kallisto']['threads']
envmodules:
"kallisto/0.46.2"
shell:
......
......@@ -37,7 +37,7 @@ rule kallisto_quant:
log:
kallisto_pseudo_log
threads:
config['pseudomapping']['threads']
config['kallisto']['threads']
envmodules:
"kallisto/0.46.2"
shell:
......
......@@ -74,13 +74,13 @@ rule sortmerna:
mv out_paired_rev.fq${{ext}} {{outnorRNA[1]}}
mv out_paired_ribo_fwd.fq${{ext}} {{outrRNA[0]}}
mv out_paired_ribo_rev.fq${{ext}} {{outrRNA[1]}}
mv output.rRNA.log {log.out}
#mv output.rRNA.log {log.out}
else
sortmerna --ref ${{fasta}} --idx-dir ${{idxdir}} --workdir {params.tmpdir}/{wildcards.SAMPLE} -threads {threads} --reads {input.fastq} --aligned output.rRNA --fastx --other output.no_rRNA -v > {log.out} 2> {log.err}
mv output.no_rRNA.fq${{ext}} {output.no_rRNA}
mv output.rRNA.fq${{ext}} {output.rRNA}
mv output.rRNA.log {log.out}
#mv output.rRNA.log {log.out}
fi
......
......@@ -43,18 +43,19 @@ rule star_index:
shell:
"""
set +o pipefail
GenomeLength=`grep -v ">" {input} | tr -d '\n' | wc -c`
GenomeLength=`grep -v ">" {input.fasta} | tr -d '\n' | wc -c`
SAindex=$(echo $GenomeLength | awk '{{a=log($1)/2 - 1; if (a>14) a=14; printf("%d\\n",a)}}')
NumberOfReferences=`grep -c ">" {input.fasta}`
ChrBin=$(echo | awk -v g=$GenomeLength -v n=$NumberOfReferences '{{a=log(g/n); if (a>18) a=18; printf("%d\\n",a)}}' )
if [[ -s {input.gff_file} && {input.gff_file} == "*.gtf" ]]
then
STAR --runMode genomeGenerate --genomeFastaFiles {input.fasta} --genomeDir {params.wkdir} --runThreadN {threads} --sjdbGTFfile {input.gff_file} --genomeSAindexNbases $SAindex
STAR --runMode genomeGenerate --genomeFastaFiles {input.fasta} --genomeDir {params.wkdir} --runThreadN {threads} --sjdbGTFfile {input.gff_file} --genomeSAindexNbases $SAindex --genomeChrBinNbits $ChrBin
elif [[ -s {input.gff_file} && {input.gff_file} == "*.gff" ]]
then
STAR --runMode genomeGenerate --genomeFastaFiles {input.fasta} --genomeDir {params.wkdir} --runThreadN {threads} --sjdbGTFfile {input.gff_file} --sjdbGTFtagExonParentTranscript Parent --genomeSAindexNbases $SAindex
STAR --runMode genomeGenerate --genomeFastaFiles {input.fasta} --genomeDir {params.wkdir} --runThreadN {threads} --sjdbGTFfile {input.gff_file} --sjdbGTFtagExonParentTranscript Parent --genomeSAindexNbases $SAindex --genomeChrBinNbits $ChrBin
else
STAR --runMode genomeGenerate --genomeFastaFiles {input.fasta} --genomeDir {params.wkdir} --runThreadN {threads} --genomeSAindexNbases $SAindex
STAR --runMode genomeGenerate --genomeFastaFiles {input.fasta} --genomeDir {params.wkdir} --runThreadN {threads} --genomeSAindexNbases $SAindex --genomeChrBinNbits $ChrBin
fi
samtools faidx {input.fasta}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment