Hi guys,
I am aligning RNA-seq data using STAR and will need in later steps to use sorted BAM files. I was wondering what people suggest to use for BAM sorting? The options that I have considerer/used include:
STAR --outSAMtype BAM SortedByCoordinate - but this crashes because of memory (we only have 128GB of RAM our server). That is when I'm merging lanes and have paired-end reads using STARSTAR --ReadFilesIn read1_lane1_r1,fq,read1_lane2_r1 read1_lane1_r2,read1_lane2_r2 This is fixable by limiting RAM STAR --limitBAMsortRAMOR
samtools sort -b -o outSorted.bam inFromSTAR.bamOR
Are there performance differences between STAR and samtools for example? If anyone knows. And Is sorting algorithms the same similar between different tools?
Cheers
•
written
by
Kirill Tsyganov
Sorting tool is your choice of preference. Either STAR, samtools sort, picard sortSam or sambamba-sort they all do the same thing. (samtools sort and sambamba are multithreaded and works much faster. Plus sambamba indexing works at lightning speed). But unsorted bam file works much faster when you are using count tools such as featureCounts. (It takes like ~5 mins for assigning 60 million reads to genes for unsorted bam file against 60 mins for position sorted, of which most of the time is spent on format conversion).
May be you do second option, output as unsorted bam files, get counts using featureCounts for reference gtf (if this is what you are up to) and sort it using any of the above tool and remove unsorted one if you want to save the disk space.
•
written
by
poisonAlien
Great, thanks for answer. I know that featureCounts is faster than htseq-count. Would you say the same for htseq-count that it works faster on the unsorted bam? thanks
Log in to answer this question.
