BWA mem and Samtools sort in Bash script
Hi, Biostars community.
I'm trying to use both BWA mem and Samtools sort at the same time via passing BWA mem output into samtools.
Command:
bwa mem reference.fa file.fq | samtools sort -o output.bam -
Well, it works fine while I'm using terminal, but the same command in bash script throws errors.
#!/bin/bash
bwa mem reference.fa file.fq | samtools sort -o output.bam -
Error message:
script.sh: line 3: 18: command not found
[main] unrecognized command 'sort'
Can someone show me the correct way to pass bwa output to samtools in bash script?
• 5,367 views
•
link
updated
by
Biostar
•
written
by
Vladislav •
•
written
by
Vladislav •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
BWA alt-aware pipeline troubleshooting
written by Duy •
-
Trying to write bwa mem -> samtools view -> samtools sort loop
written by Lee •
-
How to choose BWA-MEM2 settings
written by mgranada3 •
-
Using shell command in Nextflow DSL=2 for BWA mem aligner
written by rohitsatyam102
-
Samtools indexing problem during bwa mem alignment
written by GlouGlou •
-
Explanation of dimension differences between bam files (BWA aln and BWA mem)
written by qstefano •
-
bwa mem on lowpass wgs
written by igor •
-
samtools arror: [E::hts_idx_push] NO_COOR reads not in a single block at the end 0 -1
written by marongiu.luigi
-
Script to sort bam file not working
written by emth4657 •
-
How to improve the pipeline of call snp in RNA-seq?
written by liuhui •
do you have to load
samtoolsfirst, e.g. via spack or the like? I.e. when you log into your server can you usesamtoolsright away?also, can you post the entire script?
I've posted all script. See "#!/bin/bash"
Well, yes, samtools is alredy installed and in works fine. I guess there is a different way to pass bwa output to samtools is bash script, not via vertical line
|symbol.The pipe
|symbol should work. I can't currently replicate your issue. Can you provide more information about the environment you're working in (workstation/cluster/scheduler/HPC?), and the versions ofbash,bwa,samtoolsyou're using?It works fine in terminal, but not when I try to call the script.
Ubuntu Server 18.04 bwa: 0.7.17-r1198-dirty samtools: 1.9
I think you need to also add
-O BAMargument to write the output to fileNo need for
-O BAM, as recent versions of samtools will recognize the file format based on the extension passed with-o.Thanks WoterDeCoster for this information
what is the output of
?
Do the commands work inside a script if you separate them in to different lines? i.e. output a file with
bwathen read it separately withsamtools sort.There is no reason your pipe syntax shouldn't work, so there is something more 'sinister' going wrong behind the scenes I think.
How did you execute that script?
What makes this question interesting/amusing is that it takes a little bit of effort to write a code that deliberately produces the error. It is like a weird code golf: Write the shortest bash code that produces the error message shown above.
Here is my take:
it prints:
I feel I learned something nontrivial especially when it comes to reading error messages.
What is in my script? Is it possible to make it shorter?
three lines, four words, 16 characters:
As a response to the original poster: I believe that what ends up running is something completely different than what you think/show as being run.
To reduce the suspense:
prints:
Notably this example also helped us discover a bug in the markdown rendering for the javascript preview as well as the backend of the site.
The libraries we use won't properly render leading empty lines (two empty lines to push the error to line
3), hence I can't quite demonstrate my code correctly.