bowtie2 - (ERR): bowtie2-align exited with value 13

Question

[saifulislam99121]

I am trying to run bowtie2. but following error are occuring everytime

bowtie2 --very-fast-local -x bowtie -q -1 R1.fastq -2 R2.fastq  -s aligned.sam
Saw ASCII character 10 but expected 33-based Phred qual.
terminate called after throwing an instance of 'int'
Aborted (core dumped)
(ERR): bowtie2-align exited with value 134

[Kevin Blighe]

Why have you got -x bowtie?

What is the output of:

head R1.fastq
zcat R1.fastq | head

[saifulislam99121]

-x bowtie, here bowtie is my index name.

head R1.fastq
@NB501155:17:HFVMTAFX3:1:11101:23243:1106 1:N:0:GATTCTGAAT+GTCCACCGCT
TTGCTAAAGTTCAGATCTTGATAAAAGACATTAGAAATCACACCAGCAATAATGTCTGAAATACGGAAGAACCCAATCAAGAATAGAAGTACTAAGGCAAGTTTGACACCATAACGTTTAAAGAAATCCAGAATTGGTGAAATCCAAGT
+
AAAAA/EEEEE/EEEEEEEEEE/EEEAEEEEEEEEEEA/E/EA/EE/EEE/AEEEEEEEEAEAEE/E<AAAAE6EAEAEEAEEEE/EEAEE/AAEE/<EEEA/E<E/E//<EAEA<EEEA//<EEEEA/</EEEE/E///E/</6AE<E
@NB501155:17:HFVMTAFX3:1:11101:4115:1116 1:N:0:GAATCTGAAT+GTCCACCGCT
AATCAATGCCGTATTTGGTTTTAAATGCATGAAAAACAGCAGCTTTTAAAATTGGATTTTCACTTGCAGCAACTTTACCCACTACACGACTTAACTGATGCTGAGGAACAAGATTTTGTGCCTTAATAAATAATTCTTTTTTTAAGCGA
+
AAAAAEEEEEEAEEEEEE/EEEEEEE/EEAEEE/EEEEEEEE<E//E/EEEE/E/EE<E/EEE</EAE/AEEEAE/AA<<AEE<AAAEAE6EAEE//EEAEE//E/EEEEEEE/EE/<A<<AA/6EAAEEA<E//AEEEAEE/AAA/A<
@NB501155:17:HFVMTAFX3:1:11101:9798:1117 1:N:0:GATTCTGAAT+GTCCACCGCT
GTTTTACCTAAACCTACTTCATCAGCAAGTAAAACACGTGGTGCGAAACGTTTACCGACTTCATGTGCAATATAAAGCTGATGCGGAATTAAACCTACACGAGCCCCAACCAAGCCGCGTAATGGGCTGGTTTGCATTTGCGCCTGCAT
zcat R1.fastq | head

gzip: R1.fastq: not in gzip format

Here is the output.

Answer 1

[Mensur Dlakic]

It sounds like some of the lines are of unequal length between the sequence and quality records. That is why ASCII character 10 (end-of-line) is encountered while the expectation was it would be a 33-based Phred quality value. I suggest you try fixing the files with repair.sh which is part of the BBtools package.

[saifulislam99121]

Thank you so much..

Works fine

[ATpoint]

Also note that -s (lowercase) is not the option for saving a SAM file, that would be uppercase -S, bet better would be to save as BAM right away, bowtie2 -x (...) | samtools view -o your.bam