SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
bwa sampe error Mali Salmon Bioinformatics 14 10-27-2014 11:25 AM
BWA error seq_GA Bioinformatics 22 03-22-2013 12:09 AM
stampy error with bwa pankajsvats Bioinformatics 6 01-11-2012 06:39 AM
BWA error jyli Bioinformatics 2 06-28-2010 06:35 PM
bwa sampe 0.5.7 error? rcorbett Bioinformatics 2 04-22-2010 07:13 AM

Reply
 
Thread Tools
Old 07-27-2011, 02:49 AM   #1
Robby
Member
 
Location: Germany

Join Date: Mar 2011
Posts: 68
Default BWA error

Hello everyone,

I used BWA for mapping and received the following error message:
"Parse error at line 34: sequence and quality are inconsistent"
It is a paired end run. The mapping of forward and reverse reads worked without problems. The error message occurs when I use sampe. What could be the mistake? I am gratefull for any ideas.

Thanks Robby
Robby is offline   Reply With Quote
Old 07-28-2011, 12:17 PM   #2
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,543
Default

Have you looked at the input FASTQ file, especially at the area around line 34?

Data can get corrupted in transfer...
maubp is offline   Reply With Quote
Old 07-28-2011, 04:28 PM   #3
Richard Finney
Senior Member
 
Location: bethesda

Join Date: Feb 2009
Posts: 700
Default

What are the commands you use?
Is it in the "aln" program or the subsequent "sampe" program?
Did the aln finish?
Is there junk in the intermediate "sai" file? ( run "strings filename.sai | less" and glance for readable error messages). It's easy to accidentally send stderr to stdout in the aln phase and use this to pump into "sampe". bwa is merciless in expecting correct input.
Richard Finney is offline   Reply With Quote
Old 07-28-2011, 11:17 PM   #4
Robby
Member
 
Location: Germany

Join Date: Mar 2011
Posts: 68
Default

Thanks a lot for your answers.

I checked the fastq-file, but couldn't see any strange data. But I am not really sure, what I have to check exactly. So I had just a look, if the sequence length is the same as the quality scores and if the line breaks are correct. But everything seems to be OK. The onliest difference to other fastq-files is, that in the third line of each sequence is just a "+" instead of "+" and sequence name. But I thought that this is maybe due to the new Illumina format. Does anyone have more information regarding that line?

The error message occurs in the "sampe" program. The "aln program finished without any error message. I had a look into the sai-files, but I don't understand the output. But at least I couldn't find an error message in these files and the files have the expected size.

This was the first run with the new Illumina v3 chemistry and the new fastq format (i.e. Sanger quality scores). Do I have to change in that case anything?

I used the following commands:
bwa aln -I -l 35 ref.fasta reads1.fastq.gz > align1.sai
bwa aln -I -l 35 ref.fasta reads2.fastq.gz > align2.sai
bwa sampe ref.fasta align1.sai align2.sai reads1.fastq.gz reads2.fastq.gz > mapping.sam

I tried the aln-commands without the -I option (for Illumina 1.3+) as well, but received the same error message.
Robby is offline   Reply With Quote
Old 07-29-2011, 12:30 AM   #5
gaffa
Member
 
Location: Gothenburg/Uppsala, Sweden

Join Date: Oct 2010
Posts: 82
Default

Quote:
Originally Posted by Robby View Post
I checked the fastq-file, but couldn't see any strange data. But I am not really sure, what I have to check exactly. So I had just a look, if the sequence length is the same as the quality scores and if the line breaks are correct. But everything seems to be OK. The onliest difference to other fastq-files is, that in the third line of each sequence is just a "+" instead of "+" and sequence name. But I thought that this is maybe due to the new Illumina format. Does anyone have more information regarding that line?
That's fine, the read name is not needed at the third line, just a "+" is valid fastq and recognized by BWA.
gaffa is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:05 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO