This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Converting and sorting .sam at the same time

Hi

I have a .sam file from STAR; For me picard SortSam gives errors on my machine so I want to use samtools

Is there any command gives a sorted .bam at the same time by samtools without first converting .sam to .bam then sorting .sam?

rna-seq r next-gen
• 4,175 views • link
updated by WouterDeCoster WouterDeCoster • written by zizigolu zizigolu

1 answer 

samtools sort your_alignment.sam -o your_alignment.bam

You could have found that by reading the manual, but instead you decided to have someone else help you.

• 0 views • link
written by WouterDeCoster WouterDeCoster

Alternatively, tell STAR to produce sorted BAM output as in this manual.

• 0 views • link
written by michael.ante michael.ante

I am not a regular STAR user but perhaps STAR needs samtools to be available in your $PATH for that option to work?

• 0 views • link
written by GenoMax GenoMax

What did the log file of STAR have to say about the error? If you want to actually solve issues rather than just trying out new tools you might be unfamiliar with, it might be more useful to try to follow the suggestions by the users here.

• 0 views • link
written by Friederike Friederike

Sorry I faced a problem

[fi1d18@cyan01 310N]$ samtools sort -n ./Aligned.out.sam -o aligned1.sorted.bam
[bam_sort_core] merging from 20 files...

I repeated this command 2 times in 2 different folders but I obtained 2 .bam files with different size one 2046643 KB and another 2036831 KB

How I know which .bam is valid or intact?

• 0 views • link
written by zizigolu zizigolu

How much RAM are you assigning to this job? If you have a large SAM file you may need to assign more RAM.

• 0 views • link
written by GenoMax GenoMax

@ Friederike STAR did not give any error but output was unsorted .sam

• 0 views • link
written by zizigolu zizigolu

You may also need to use --limitBAMsortRAM if you don't have enough RAM available.

• 0 views • link
written by GenoMax GenoMax

Log in to answer this question.

More posts like this