Convert SAM to BAM error
Hi all,
I have mapped a pair-end SAMPLE_R1.fastq and SAMPLE_R2.fastq to hg.38 in Bowtie2. The mapping is fine as I could see the summary stats from the stederr file.But when I tried to convert the SAM file to BAM file.
samtools view -b -o SAMPLE.bam -S SAMPLE.sam
I have this error message:
[W::sam_read1] Parse error at line 198
samtools view: error reading file "SAMPLE.sam"
This is line 198 in the sam file
@SRR891270.1
HWI-ST281:266:C1LTTACXX:2:1101:1339:1962/1%0ANTATCACACCCCATCCTAAAGTAAGGTCAGCTAAATAAGCTATCGGGCCC%0A+%0A#41BABDDDHHFHCEEEAAFAFCCCGFF<@?E:*1?DDFDGGFFHGIH<F%0A
What problem of this line would be?
Thanks
1 answer
That line is the four lines of a FASTQ record, with the four newlines percent-encoded as %0A. It is very far from being a valid SAM line.
You would need to look at the rest of your “SAM” file and at the commands that produced it to figure out what has gone wrong here.
• 0 views
•
link
written
by
John Marshall
Log in to answer this question.
More posts like this
-
samtools convert SAM to BAM error code
written by gbmygbm •
-
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
-
[main_samview] fail to read the header from sample.bam
written by sambioinfo2018 •
-
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 •
-
Select concordantly mapped reads from a sam file without header
written by valerie •
-
Samtools Adding Head reduces file size
written by martijnkelder •
-
Preprocessing for GATK
written by anailg92 •
I'm currently learning on the SAM to BAM conversion. From my understanding the code line should be;
instead of samtools view -b -o SAMPLE.bam -S SAMPLE.sam
Have a try on it and let's see if this works (mine works well). Hope this helps!
Any time you see
samtools view -bS, it is an indication that the tutorial you are reading is outdated.Your suggested command line uses the
-b,-S, and-ooptions, which is exactly the same set of options that BeeWork is already using. So it is effectively identical to the OP's version (except that it puts an option after the filename argumentSAMPLE.sam, which is a bad habit especially in scripts as it does not work on all Unixes).Samtools has automatically detected the format of its input files since 2014. The
-Soption (forsamtools view), which previously specified that the input file was SAM rather than BAM, has been ignored since samtools 1.0. There is no reason to use it.Thanks John,
Suchakan - John is correct. these flags are the same and don't need -S to specific input SAM format, samtools view should be able to detect. in samtools-view manual, it said "ignored" in -S flag.
Oh thank you so much for telling me this invaluable information; I have learned so much from you two John and BeeWork. Sorry for the inconvenience caused.