diff --git a/alignment/rnaseqalign.sh b/alignment/rnaseqalign.sh index 0ed40366290a8abea1bb3bc98d9515f8ead7b806..b4cb8851ece1c99f502b735d7db9b11e560993db 100644 --- a/alignment/rnaseqalign.sh +++ b/alignment/rnaseqalign.sh @@ -40,9 +40,12 @@ if [[ -z $SLURM_CPUS_ON_NODE ]] then SLURM_CPUS_ON_NODE=1 fi + +diff $fq1 $fq2 > difffile.out + if [ $algo == 'star' ] then - if ($fq1 != $fq2) + if (-s difffile) then module load star/2.4.2a STAR --genomeDir ${index_path}/star_index/ --readFilesIn $fq1 $fq2 --readFilesCommand zcat --genomeLoad NoSharedMemory --outFilterMismatchNmax 999 --outFilterMismatchNoverReadLmax 0.04 --outFilterMultimapNmax 20 --alignSJoverhangMin 8 --alignSJDBoverhangMin 1 --alignIntronMin 20 --alignIntronMax 1000000 --alignMatesGapMax 1000000 --outSAMheaderCommentFile COfile.txt --outSAMheaderHD @HD VN:1.4 SO:coordinate --outSAMunmapped Within --outFilterType BySJout --outSAMattributes NH HI AS NM MD --outSAMstrandField intronMotif --outSAMtype BAM SortedByCoordinate --quantMode TranscriptomeSAM --sjdbScore 1 --limitBAMsortRAM 60000000000 --outFileNamePrefix out @@ -53,7 +56,7 @@ then mv outAligned.sortedByCoord.out.bam output.bam else module load hisat2/2.1.0-intel - if [ $fq1 == $fq2 ] + if [ -s difffile ] then hisat2 -p $SLURM_CPUS_ON_NODE --rg-id ${pair_id} --rg LB:tx --rg PL:illumina --rg PU:barcode --rg SM:${pair_id} --no-unal --dta -x ${index_path}/hisat_index/genome -U $fq1 -S out.sam --summary-file ${pair_id}.alignerout.txt else diff --git a/preproc_fastq/trimgalore.sh b/preproc_fastq/trimgalore.sh index b941ef1b9710831f9f8fbcf2af34949b84e6733c..1ae8ba12606bc808337fcea4e3fbe6fe7b28fea1 100644 --- a/preproc_fastq/trimgalore.sh +++ b/preproc_fastq/trimgalore.sh @@ -32,13 +32,13 @@ r2base="${fq2%.fastq*}" source /etc/profile.d/modules.sh module load trimgalore/0.4.1 cutadapt/1.9.1 -if [ $fq1 == $fq2 ] +if [ -s $fq2 ] then - trim_galore -q 25 --illumina --gzip --length 35 ${fq1} - mv ${r1base}_trimmed.fq.gz ${pair_id}.trim.R1.fastq.gz - cp ${pair_id}.trim.R1.fastq.gz ${pair_id}.trim.R2.fastq.gz -else trim_galore --paired -q 25 --illumina --gzip --length 35 ${fq1} ${fq2} mv ${r1base}_val_1.fq.gz ${pair_id}.trim.R1.fastq.gz mv ${r2base}_val_2.fq.gz ${pair_id}.trim.R2.fastq.gz +else + trim_galore -q 25 --illumina --gzip --length 35 ${fq1} + mv ${r1base}_trimmed.fq.gz ${pair_id}.trim.R1.fastq.gz + cp ${pair_id}.trim.R1.fastq.gz ${pair_id}.trim.R2.fastq.gz fi