diff --git a/alignment/rnaseqalign.sh b/alignment/rnaseqalign.sh index 4ef42998129846b5c56938ddee0e8207fc6ff420..182a738385e8c8bf3159c151a498df76ec89c2d6 100644 --- a/alignment/rnaseqalign.sh +++ b/alignment/rnaseqalign.sh @@ -69,8 +69,9 @@ else else samtools view -1 --threads $SLURM_CPUS_ON_NODE -o output.bam out.sam fi - samtools sort -@ $SLURM_CPUS_ON_NODE -O BAM -n -o output.nsort.bam output.bam + samtools sort -@ $SLURM_CPUS_ON_NODE -O BAM -o ${pair_id}.bam output.bam java -jar $PICARD/picard.jar FixMateInformation ASSUME_SORTED=TRUE SORT_ORDER=coordinate ADD_MATE_CIGAR=TRUE I=output.nsort.bam O=${pair_id}.bam + fi samtools index -@ $SLURM_CPUS_ON_NODE ${pair_id}.bam