SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
error in picard (SAM validation error) dongshenglulv Bioinformatics 3 07-05-2013 12:58 PM
SAM to FASTQ converter - Picard sdm Bioinformatics 14 03-19-2013 10:19 AM
Picard error JohnK Bioinformatics 14 12-14-2012 05:41 AM
picard error seq_GA Bioinformatics 1 11-25-2010 09:47 PM
N character in Illumina fastq format balamudiam Illumina/Solexa 1 12-13-2009 11:53 PM

Reply
 
Thread Tools
Old 02-22-2012, 08:19 PM   #1
sehrrot
Member
 
Location: Australia

Join Date: Jul 2010
Posts: 55
Default Picard Error: IllegalArgumentException: Invalid fastq character:

hi all

I'm having trouble in using Picard to covert sam to bam. it seems to be halted by weird characters (I cannot type it here)

Quote:
Ignoring SAM validation error due to lenient parsing:
Error parsing text SAM file. length(QUAL) != length(SEQ); Line 28
Line: DFCDZDN1:130:B01D5ACXX:4:1101:1099:2046 99 chr3 134076497 60 100M = 134076563 166 AGGGCAGAGGAGGGGAGAGAATGGCTCAGAGTCCTGAAGCACCACCTCCTTCAGATCAGTATCTCCACCATGTCTGACAAGTCCTGGANTCTGGATGTAG $$$''''%)))))++*+((**+*+++++**+'))++**++**+**+(*)++++))))()'''''''%&%%$$"$&%$$$$%"$%$%% (weird character here)
[Thu Feb 23 15:13:19 EST 2012] net.sf.picard.sam.SortSam done. Elapsed time: 0.00 minutes.
Runtime.totalMemory()=123666432
Exception in thread "main" java.lang.IllegalArgumentException: Invalid fastq character: (weird character here)
at net.sf.samtools.SAMUtils.fastqToPhred(SAMUtils.java:408)
at net.sf.samtools.SAMUtils.fastqToPhred(SAMUtils.java:387)
at net.sf.samtools.SAMRecord.setBaseQualityString(SAMRecord.java:254)
at net.sf.samtools.SAMTextReader$RecordIterator.parseLine(SAMTextReader.java:422)
at net.sf.samtools.SAMTextReader$RecordIterator.next(SAMTextReader.java:278)
at net.sf.samtools.SAMTextReader$RecordIterator.next(SAMTextReader.java:250)
at net.sf.samtools.SAMFileReader$AssertableIterator.next(SAMFileReader.java:641)
at net.sf.samtools.SAMFileReader$AssertableIterator.next(SAMFileReader.java:619)
at net.sf.picard.sam.SortSam.doWork(SortSam.java:67)
at net.sf.picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:177)
at net.sf.picard.cmdline.CommandLineProgram.instanceMainWithExit(CommandLineProgram.java:119)
at net.sf.picard.sam.SortSam.main(SortSam.java:79)
sehrrot is offline   Reply With Quote
Old 02-22-2012, 09:15 PM   #2
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 868
Default

Either the fastq was made wrongly, or the .sam was made wrongly. Those are not normal quality characters.

Last edited by swbarnes2; 02-22-2012 at 09:18 PM.
swbarnes2 is offline   Reply With Quote
Old 02-22-2012, 11:57 PM   #3
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,303
Default

What is the weird characters? Try using hexdump.
maubp is offline   Reply With Quote
Old 02-23-2012, 05:07 PM   #4
sehrrot
Member
 
Location: Australia

Join Date: Jul 2010
Posts: 55
Default

Hi all

now it's working fine after I removed "-I" from bwa alignment process ("bwa aln -t 4 -f input.sai -I hg19 input.fastq")
sehrrot is offline   Reply With Quote
Old 02-24-2012, 12:30 AM   #5
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,303
Default

So you probably have Sanger style FASTQ, which is what the latest Illumina pipeline produces - but you told bwa it was the legacy Illumina specific FASTQ encoding.

See:

http://en.wikipedia.org/wiki/FASTQ_format
http://dx.doi.org/10.1093/nar/gkp1137
maubp 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 06:38 AM.


Powered by vBulletin® Version 3.8.6
Copyright ©2000 - 2014, Jelsoft Enterprises Ltd.