This is a test version of Biostars. For the public version, visit https://www.biostars.org.
samtools Mpileup quality scores?

Dear all,

What is the default quality scores expected by samtools mpileup?

I see you can specify the option

-6, --illumina1.3+ Assume the quality is in the Illumina 1.3+ encoding.

Is the default therefore Sanger?

And what do you specify if you have Illumina 1.9 where the quality scores are different to Illumina 1.3??

Thanks!

mpileup quality-scores samtools
• 4,465 views • link
updated by Ram Ram • written by AW AW

Sorry I should have been clearer.

Does it expect phred33 or phred64?

• 0 views • link
written by AW AW

1 answer

All modern tools will expect Sanger encoding. It is best to convert all data to that.

See also:

Convert Illumina Reads To Sanger Score

• 0 views • link
written by Istvan Albert Istvan Albert

Thanks I meant does it expect Sanger phred33, illumina 1.3 phred 64 or illumina 1.8 phred 33?

• 1 views • link
written by AW AW

Thanks!

So how do I specify illumina phred33?

• 0 views • link
updated by Ram Ram • written by AW AW

actually most of the time you may not even need to do anything for illumina phred 33 data

the offset is the same, 33, that is what throws tools off,

the encoding differs only in the largest value 40 vs 41 that on its own matters only if someone wrote a very sloppy tool that can't think outside of 40

• 1 views • link
updated by Ram Ram • written by Istvan Albert Istvan Albert

Log in to answer this question.

More posts like this
  • Why does Phred+ encoding changes after trimming fastq files??
    written by Bikal • Bikal
    When I did fastqc analysis on my original fastq files, the encoding that I can see in the biostatistics was sanger/illumina 1.9. But when I …
  • mpilup output quality
    written by mluypaert • mluypaert
    Hi All I was look at the `samtools mpileup` output format and found something odd in the reported quality scores. 1 633887 G 28 >><<<<CAAA.AA.AAAcaAAAC.A..a …
  • Scafolding and gap filling
    written by bhagyathimmappa • bhagyathimmappa
    Hello all, for one of the novel fungi, we did illumina mi-seq paired end sequencing and used velvet to assemble, I have 15500 contigs. some …
  • RSEM ignores quality scores?
    written by AW AW
    Dear all, I want to quantify gene expression for my Trinity transcriptome assembly using RSEM. I am getting confusing results depending on how I specify …
  • Can Bowtie 0.12.7 Handle Sanger / Illumina 1.9 ?
    written by fire_water fire_water
    Hello. I have a user that needs to use Bowtie 0.12.7 (Yeah, old version. Don't ask.) but the quality score format of our FASTQ files …
  • Fastx Invalid Quality Score Value
    written by sckinta sckinta
    <p>I have multiple RNAseq libraries to parse, so I wrote a pipeline using bash and submitted them in batch. Most of libraries ran well and …
  • 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> …
  • Changing Quality Scores Within Bam Files
    written by matted matted
    <p>I have a set of <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> files that include reads with quality scores in a mix of <a href='http://en.wikipedia.org/wiki/FASTQ_format'>incompatible formats</a>. This includes Solexa, Illumina …
  • 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 …
  • Fastq Format Problem When Using Bwa Aln -Q Option (Read Trim)
    written by Bioscientist Bioscientist
    <p>This is kind of follow-up question of <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>My paired-end reads are 100bp, and in Illumina 1.5+ format. Now I need to run like:</p> …