SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
BWA sam and Samtools sam->bam conversion problem maasha Bioinformatics 6 06-05-2013 08:39 AM
Sam / Bam to Ace conversion empyrean Bioinformatics 4 05-13-2013 05:32 AM
sam 2 bam conversion error DavyK Bioinformatics 2 01-14-2013 07:13 PM
Issue with Sam-Bam conversion samtools - how to remove last line of Sam file? TabeaK Bioinformatics 3 11-19-2012 11:05 AM
SAM to BAM conversion issue rndouglas Bioinformatics 2 08-27-2012 04:10 PM

Reply
 
Thread Tools
Old 01-25-2017, 11:29 PM   #1
finswimmer
Member
 
Location: Europe

Join Date: Oct 2016
Posts: 60
Default Some clarification about sam to bam conversion

Hello,
for getting a sorted bam file out of an unsorted sam file i often found these commands:

Code:
samtools view -bS file.sam | samtools sort - > file_sorted.bam
samtools view -bT reference.fa file.sam | samtools sort - > file_sorted.bam
My questions are:
1. When do I need to give the reference file as an argument?
2. Are there any disadvantages if I do this:
Code:
samtools sort -O BAM file.sam > file_sorted.bam
fin swimmer
finswimmer is offline   Reply With Quote
Old 01-26-2017, 12:22 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

  1. You only need to do that with CRAM files (and usually not even then)
  2. Samtools sort doesn't accept BAM input.

As an aside, you can use the following for a slight speed improvement:

Code:
samtools view -uS file.sam | samtools sort - > file_sorted.bam
Note that the syntax for "samtools sort" has changed in more recent versions.
dpryan is offline   Reply With Quote
Old 01-26-2017, 12:48 AM   #3
finswimmer
Member
 
Location: Europe

Join Date: Oct 2016
Posts: 60
Default

Hello
Quote:
Originally Posted by dpryan View Post
  1. You only need to do that with CRAM files (and usually not even then)
  2. Samtools sort doesn't accept BAM input.
you mean doesn't accept SAM input?

Code:
samtools 1.3.1
Using htslib 1.3.1
Copyright (C) 2016 Genome Research Ltd.

Usage: samtools sort [options...] [in.bam]
Options:
  -l INT     Set compression level, from 0 (uncompressed) to 9 (best)
  -m INT     Set maximum memory per thread; suffix K/M/G recognized [768M]
  -n         Sort by read name
  -o FILE    Write final output to FILE rather than standard output
  -T PREFIX  Write temporary files to PREFIX.nnnn.bam
  -@, --threads INT
             Set number of sorting and compression threads [1]
      --input-fmt-option OPT[=VAL]
               Specify a single input file format option in the form
               of OPTION or OPTION=VALUE
  -O, --output-fmt FORMAT[,OPT[=VAL]]...
               Specify output format (SAM, BAM, CRAM)
      --output-fmt-option OPT[=VAL]
               Specify a single output file format option in the form
               of OPTION or OPTION=VALUE
      --reference FILE
               Reference sequence FASTA FILE [null]
The strange thing is, that I didn't see any problems until now by sorting the way I described above...

fin swimmer
finswimmer is offline   Reply With Quote
Old 01-26-2017, 12:57 AM   #4
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Oops, yes I meant SAM input. That at least used to not work, but perhaps it's been silently changed.
dpryan is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 05:49 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO