Dear BioStarers,
When I convert my sam file to bam file, I got "truncated file" in prompt, and also "[sam_read1] reference '503821040998320619' is recognized as '*'". It is no such problem, when I converted other sam
files to bam using the same reference, here is "../Alyrref.folded.fa.fai". Now, I do not know what I can do next.
So, I need to get helps from you. What have happened? Is it a serious problem? Can I still use this truncated file for pileup? Thanks in advance.
My command and the prompt:
$ samtools view -bht ../Alyrref.folded.fa.fai map.list.sam > map.list.bam
[sam_header_read2] 1099 sequences loaded.
[sam_read1] reference '503821040998320619' is ecognized as '*'.
[main_samview] truncated file.
•
written
by
Jianfengmao •
1 answer
I remember that I had a similar problem a while ago. In my case it was because the sam file was faulty: bwa was terminated because the qsub time allocation was just too short for the massive amount of data that came from this particular lane (cluster density was way too high). Running your mapping program again to get an intact sam file might solve this problem.
Log in to answer this question.
-
Truncated sam File from Resfinder Database in MGmapper
written by hafiz.talhamalik
-
Parse error at line.....
written by priscilladraj •
-
problem in conversion of sam to bam file
written by priya120195 •
-
How to rectify error while sorting the sam file
written by sunnykevin97
-
Error converting sam to bam
written by valopes •
-
Converting .sam to .bam with samtools and I get a new error with every repeated command
written by ryanzapotocny •
-
failed to convert from cram to bam
written by kygppyou •
-
Select concordantly mapped reads from a sam file without header
written by valerie •
-
Tophat-Fusion Sam To Bam Error
written by Nicolas Rosewick
-
Truncating Bam Files With Samtools
written by User 9996
I suggest you edit your question and include the last few lines of your SAM file - it sounds like the last line in the file is truncated, but maybe something else is going on. Use the tail command to get the final lines, e.g. "tail -3 map.list.sam".