This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Fastq Format Problem When Using Bwa Aln -Q Option (Read Trim)

This is kind of follow-up question of http://biostar.stackexchange.com/questions/14723/trim-the-low-quality-end-of-100bp-read

My paired-end reads are 100bp, and in Illumina 1.5+ format. Now I need to run like:

bwa aln -q 20

I'm just wondering do I need to first convert Illumina 1.5+ format to Sanger format? Or simply set -q 46 ? (Since Sanger format is phred+33; while Ilumina 1.5+ is phred+64)

Actually I'm a little bit confused about Sanger, Illumina 1.3, Illumina 1.5. What are these? Where do they come from?

Thanks

Edit: I just notice bwa aln has an "I" option, which converts Illumina 1.3 to sanger. But what about Illumina 1.5? Can it be applied to 1.5?

bwa trimming quality
• 5,126 views • link
updated by Damian Kao Damian Kao • written by Bioscientist Bioscientist

1 answer

Here is a good thread on seqanswers about converting Illumina to Sanger quality scores.

Here is a good description of what the quality scores are: http://en.wikipedia.org/wiki/FASTQ_format#Quality

• 0 views • link
updated by Ram Ram • written by Damian Kao Damian Kao

Log in to answer this question.

More posts like this
  • Filtering sanger_fastq, illumina_fastq files using Fastx-Toolkit
    written by bioinfo bioinfo
    I am doing some quality filtering of a large bunch of mixed fastq files produced from multiple versions of Illumina platforms. Thus, the quality scores …
  • Bwa Mem: How To Specify The Fastq Phred Format?
    written by Giovanni M Dall'Olio Giovanni M Dall'Olio
    <p>Hello,</p> <p>I have a dataset of short reads in which some fastq files are in the Illumina 1.5 format, and others in the Illumina 1.8. …
  • Convert Illumina 1.3 To Illumina 1.5
    written by Giovanni M Dall'Olio Giovanni M Dall'Olio
    <p>I have some fastq files in the Illumina 1.3 format. What is the best way to convert it to Illumina 1.5?</p> <p>Quoting wikipedia, the differences …
  • Presence Of Series Of Hashes In Illumina Quality Line
    written by deepthithomaskannan deepthithomaskannan
    <p>Hi all,</p> <p>My paired-end illumina data contains series of hash(#) signs in the quality line. Some sequence have # quality for all the bases. When …
  • Fastq Format Confusion
    written by Biomonika (Noolean) Biomonika (Noolean)
    <p>Hi!</p> <p>In the following article <a href='http://file:///home/misa/Sta%C5%BEen%C3%A9/NGS_workshop_1_Intro.pdf'>The Sanger FASTQ file format for sequences with quality scores, and the Solexa/Illumina FASTQ variants.</a> I found this table:</p> …
  • What Is The Default Quality Encoding Expected By Bwa?
    written by Panos Panos
    <p>What is the quality encoding in the input reads that <a href='http://bio-bwa.sourceforge.net/'>BWA</a> expects as default? Is it Sanger, Solexa, Illumina 1.3+, Illumina 1.5+ or Illumina …
  • Best Way To Preprocess Barcoded Illumina Paired-End Data
    written by Ying W Ying W
    <p>I know this topic has come up <a href='http://biostar.stackexchange.com/questions/1675/filtering-paired-end-reads'>many</a> <a href='http://biostar.stackexchange.com/questions/8791/paired-end-sequence-preprocessing-effects-on-downstream-processes'>times</a> <a href='http://biostar.stackexchange.com/questions/3081/multiplexed-pair-end-sequences-barcode-sorting'>before</a> but I wanted to get feedback on what people are doing now …
  • How Mapq In Bwa Is Calculated?
    written by Bioscientist Bioscientist
    <p>I'm just curious about this. I searched a lot and see similar posts but seems the author doesn't give out detailed explanations in their manual, …
  • Removal Of Low Quality Reads With Trailing "B"
    written by Bioscientist Bioscientist
    <p>This is kind of follow-up thread for: <a href="http://biostar.stackexchange.com/questions/14723/trim-the-low-quality-end-of-100bp-read">http://biostar.stackexchange.com/questions/14723/trim-the-low-quality-end-of-100bp-read</a></p> <p>I tried to run <code>bwa -aln -I -q 20</code> to trim off the low-quality sequences at …
  • Trim The Low-Quality End Of 100Bp Read
    written by Bioscientist Bioscientist
    <p>I usually first filtered the fastq read first, tossing out those reads with average phred score&lt; say, 20. But recently I realized long read (100bp) …