Converting and sorting .sam at the same time
1
0
Entering edit mode
5.3 years ago
zizigolu ★ 4.3k

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 • 3.3k views
ADD COMMENT
5
Entering edit mode
5.3 years ago
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.

ADD COMMENT
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

Thank you I tried telling STAR but giving nothing

C: How I know a file is sam or bam, sorted or unsorted

ADD REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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.

ADD REPLY
0
Entering edit mode

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?

ADD REPLY
1
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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

ADD REPLY
1
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 2441 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