diff --git a/alignment/markdups.sh b/alignment/markdups.sh index 565fac72b6a340d96d75242e0ff2d9eec27b7385..4b897a26be38f1f93717c3f540509014eb81b2c7 100644 --- a/alignment/markdups.sh +++ b/alignment/markdups.sh @@ -56,8 +56,8 @@ elif [ $algo == 'fgbio_umi' ] then module load fgbio bwa/intel/0.7.15 samtools index -@ $SLURM_CPUS_ON_NODE ${sbam} - fgbio GroupReadsByUmi --tmp-dir ./ -s identity -i ${sbam} -o ${pair_id}.group.bam --family-size-histogram ${pair_id}.umihist.txt -e 0 -m 0 - fgbio CallMolecularConsensusReads --tmp-dir ./ -i ${pair_id}.group.bam -p consensus -M 1 -o ${pair_id}.consensus.bam -S ':none:' + fgbio --tmp-dir ./ GroupReadsByUmi -s identity -i ${sbam} -o ${pair_id}.group.bam --family-size-histogram ${pair_id}.umihist.txt -e 0 -m 0 + fgbio --tmp-dir ./ CallMolecularConsensusReads -i ${pair_id}.group.bam -p consensus -M 1 -o ${pair_id}.consensus.bam -S ':none:' samtools index ${pair_id}.consensus.bam samtools fastq -1 ${pair_id}.consensus.R1.fastq -2 ${pair_id}.consensus.R2.fastq ${pair_id}.consensus.bam gzip ${pair_id}.consensus.R1.fastq