I am attempting to run a file through an algorithm I have been using, HLA*LA. On running the samtools command within the algorithm, I have unfortunately been getting this error. After trying to debug this following other guides, I am seeking assistance here as I have yet to find a solution.
The script is attempting to run the following command:
samtools view -@ 3 IN.cram '*' | awk '{if ($3 == "*") print $0}' | samtools view -bo UNMAPPED.bam -
The specific error that I am getting is the following:
[main_samview] fail to read the header from "-".
When I print out the header using samtools, there is output and the error is not reproduced. When I solely run the first samtools command and the awk command, there is no error or output meaning the issue is with the last samtools command.
Update: Outside of the algorithm, I was able to produce a bam output without issues when using:
samtools -b IN.cram > OUT.bam
1 answer
Error is because your awk code is stripping out the header lines (generated by samtools view -h). The header lines are needed when you try to write a valid BAM file with second view command.
You can do this.
Grab the header from your file
samtools view -H your.cram > header.file
Filter your file using your original command
samtools view -@ 3 IN.cram '*' | awk '{if ($3 == "*") print $0}' > filtered.sam
Cat the header to top of this filtered file and create your BAM.
cat header.file filtered.sam | samtools view -bo UNMAPPED.bam -
Log in to answer this question.
-
Samtools merge bam issue with header tags
written by sam •
-
SAMtools error closing standard output: -1
written by theadrijasaha •
-
awk: cmd. line:1: fatal: division by zero attempted when calculating average depth when using subpr…
written by Manuel •
-
Unable to convert from sam to bam file.
written by usef •
-
Samtools - [main_samview] fail to read the header
written by kstangline •
-
Samtools reheader Numerical result out of range
written by Medhat
-
[main_samview] fail to read the header from sample.bam
written by sambioinfo2018 •
-
bgzf_read error during samtools run
written by kk.mahsa
-
parse error with samtools
written by a.rex
-
Select concordantly mapped reads from a sam file without header
written by valerie •
You need to add
samtools view -hI edited the command to:
The error occurred again when running the second samtools command. If helpful, I am using samtools version 1.10
Your
awkstep is filtering most of lines out since you are only printing those which match*in third field. You will either need to let the header lines pass through or filter the file with first two steps. Add the header to the filtered file (you can grab it from original by doingsamtools view -Hand then create the BAM you need.The code is meant to identify and snag all the unmatched reads for paired-end data. Would this error just mean that there are no unmatched reads? The file went through a decent amount of QC.