Sorted bam files are empty after sorting them from bam
0
0
Entering edit mode
4 months ago
mgranada3 ▴ 30

Hi,

I have been working with all my DNA analysis files in parallels but I got to a point where I had about 15 files get stuck on one step. Specifically, I notice something is wrong because the files have 0 memory in them when they go from bam to sortedbam. I have compared all the files during previous steps, bam, sam, and fastq. The ones which are stuck are of a comparable size as the rest, up until they get sorted.

enter image description here

This is the code I am using to sort them.

#!/bin/bash
#SBATCH -J bamSORT_all
#SBATCH -A cluster_username
#SBATCH -N 1 --ntasks-per-node=24
#SBATCH --mem-per-cpu=8G
#SBATCH -t 36:00:00
#SBATCH -o Report-%j.out

cd $SLURM_SUBMIT_DIR
ml picard/3.0.0

for i in $(ls *.bam | rev | cut -c 5- | rev | uniq); do java -jar /usr/local/pace-apps/manual/packages/picard/3.0.0/build/libs/picard.jar SortSam -I /scratch/bam_try2/${i}.bam -O /scratch/bam_try2/${i}.sorted.bam --SORT_ORDER coordinate; done

I have tried repeating the bam to sorted bam step, sam to bam, and bwa-mem to sam. By the time I get to the sortedbam step it will show the files are empty. What could be going on?

bam dna • 350 views
ADD COMMENT
0
Entering edit mode

check the logs stdout and stderr should have some information in my opinion

ADD REPLY

Login before adding your answer.

Traffic: 1775 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6