This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Error in sam to bam conversion after bowtie alignment

I used this command line in bowtie:

bowtie -v 0 -m 1 ref file1.fastq file2.sam

and now I'm trying to convert the sam file I obtained:

samtools view -b -S file2.sam > file2.bam

and I get this error (even if I specify the -t ref.fa.fai option):

[W::sam_read1] parse error at line 1
[main_samview] truncated file.

Here's the header of my sam file:

HWI-1KL110:149:HAUWDADXX:1:1101:1419:2097 1:N:0:ATCACG  -       MT      7644    TATCACCTTTCATGATCACGCC  IEIFFIIIIIEEDDDDD8DDD?  0
HWI-1KL110:149:HAUWDADXX:1:1101:2046:2200 1:N:0:ATCACG  +       17      57918633        TAGCTTATCAGACTGATGTTGACT        CCCFFFFFHHDFHIIIIJJIJJJI        0
HWI-1KL110:149:HAUWDADXX:1:1101:3118:2245 1:N:0:ATCACG  +       17      28444112        TGAGGGGCAGAGAGCGAGACTTT @@?DDD1D8CCDFE<FG6??CDE 0

Can someone help me?

sam bowtie samtools bam
• 8,118 views • link
updated by Ram Ram • written by Maguelonne • Maguelonne

Hi. Have you tried the following command ? :

samtools view -bhS file2.sam > file2.bam

(with -h option)

• 1 views • link
written by Thibault D. Thibault D.

1 answer

Your output is not SAM. I think it's bowtie legacy output (another tsv format). You need to use bowtie --sam

• 1 views • link
updated by Ram Ram • written by thackl thackl

There are 2 possibilities in Bowtie to get an output in SAM format:

  • Add the tag -S <output.sam> if one wants directly the SAM file (space consuming)
  • In case one wants to pipe the sam output directly in samtools to get a sorted bam file (smaller), add flag --sam such as bowtie --sam | samtools view -bS - > <output.sorted.bam>

NB: Most parameters are omitted in the bowtie command.

• 0 views • link
written by Johan Zicola Johan Zicola

Log in to answer this question.

More posts like this
  • Low alignment score for microRNA
    written by ntsopoul • ntsopoul
    Hi, I am trying to align reads to miRbase however, I get very low alignment scores and I do not know why. I read in …
  • Picard ERROR:INVALID_VERSION_NUMBER of sorted BAM file
    written by marongiu.luigi marongiu.luigi
    Dear all, I have checked some sorted BAM files with [Picard 2.18.13][1] and I got this error: ## HISTOGRAM java.lang.String Error Type Count ERROR:INVALID_VERSION_NUMBER 1 …
  • Fastq File Without Index?
    written by lsvijfhuizen • lsvijfhuizen
    <p>Dear All,</p> <p>I have a fastq file which failed to align with bowtie. The error message is: reads file does not look like a FASTQ …
  • Parsing Fastq Files
    written by deepthitheresa • deepthitheresa
    <p>Hi all,</p> <p>I have Fastq reads something like </p> <pre><code>@HWI-ST1162:73:C0KEFACXX:6:1101:1816:1918 1:N:0:CGATGT NACCCTAGAAATTATAAATCTCTTCAAGTGAGATTGTAAGGAGAAGGAGAAACTTGGTCTGGAATTTGTTATAAAAGCACTT + #1=DDFFFHHGGHIJJJJJIJJJJJJJJCHGHIIJJEFHIJIJJIIJIIIIJHHIJJFHIIJJJJJJJIJIJIJIIJHEHHHHFFFFFFEEEDEEEDCDDC </code></pre> <p>I aligned this fastq file with a reference genome using …