diff --git a/alignment/dnaseqalign.sh b/alignment/dnaseqalign.sh index 0e60add806be8ac7c29243512146683a63dce889..12955d6f64a32627029ebed0d0b544602e06d899 100644 --- a/alignment/dnaseqalign.sh +++ b/alignment/dnaseqalign.sh @@ -47,7 +47,7 @@ else bwa mem -M -t $SLURM_CPUS_ON_NODE -R "@RG\tID:${pair_id}\tLB:tx\tPL:illumina\tPU:barcode\tSM:${pair_id}" ${index_path}/genome.fa ${fq1} ${fq2} > out.sam fi -if [[ $umi=='umi' ]] +if [[ $umi == 'umi' ]] then k8 /cm/shared/apps/bwakit/0.7.15/bwa-postalt.js -p ${pair_id}.hla ${index_path}/genome.fa.alt out.sam | python ${baseDir}/add_umi_sam.py -s - -o output.unsort.bam elif [[ $index_path == '/project/shared/bicf_workflow_ref/GRCh38' ]]