Tweet
Is it possible to output bwa output directly into samtools view to convert into bam format without using intermediate files? I am running low on hard disk space and I am attempting to automate my script. I would like to convert my bwa output to bam, sort it, and index it. I know the sam-bam conversion can be piped into the sort command, but is it possible for the samtools view to take its input from STDIN?
-
-
figured out on linux at least, if I give it the filename /dev/stdin that will satisfy it since the OS does the work, even though samtools does not appear to be expressly programmed to switch to stdin.
It appears to be working too, which is going to save me some headaches.
--edit
apparently I could also have done this with named pipesLast edited by Kotoro; 07-18-2011, 04:06 PM. -
bwa + samtools have been developed with pipes in mind:Originally posted by Kotoro View Post
HTHCode:$ bwa aln [OPTIONS] [DB] [FASTQ] | bwa samse [OPTIONS] [DB] - [FASTQ] | samtools view -Su - | samtools sort - [PREFIX]
dComment
-
-
Almost, this should work (essentially it is a double process substitution).Originally posted by gavin.oliver View Post
Code:$ bwa sampe [OPTIONS] [DB] <(bwa aln [OPTIONS] [DB] [FASTQ1]) <(bwa aln [OPTIONS] [DB] [FASTQ2]) [FASTQ1] [FASTQ2] | samtools view -Su - | samtools sort - [PREFIX]
Comment
-
Thanks!
The wonders of the command line, eh?
Comment
-
Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.
The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
...06-02-2026, 10:05 AM -
With the launch of new single-cell sequencing platforms in 2026, the field stands at an exciting inflection point. This article surveys the most impactful advances in the field and discusses how they’re reshaping research in cancer, immunology, and beyond.
Introduction
Single-cell sequencing technologies have undergone remarkable advances over the past decade, transitioning from low-throughput experimental approaches to highly scalable platforms capable of...05-22-2026, 06:42 AM
Comment