Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • different Illumina convention in fastq files?

    We have a few different Illumina data sets here, and I'm wondering if there has been a change in the characters used in fastq files, or that we have just been given bad data.

    we have some sequences for which the translation works fine:
    10*log(1+10**(c - 64/10.0))/log(10)

    @SLXA-EAS1_89:3:1:715:750/1
    GTCTTGAAAGCTATGATGTCAAGATTAATTTAATC
    +SLXA-EAS1_89:3:1:715:750/1
    bbbbbbbbbbbbbbbbb^bbbbbbbbbbbbbbbba
    @SLXA-EAS1_89:3:1:747:756/1
    GTGTATTGCTCAATCTTCGAACGGGGGGAGGATTG
    +SLXA-EAS1_89:3:1:747:756/1
    bbbbb^bbbbbbbbbbbbbbbbbbbbbb\bbObbb
    @SLXA-EAS1_89:3:1:859:343/1
    GTTAATAGATTTAATTGCCACCGCAATACCAGCAA
    +SLXA-EAS1_89:3:1:859:343/1
    ccccccccccccccccbcccccccbbbbbbbYbb^

    However we have other sequences where the quality is 10 or below using the same translation. Is one using Sanger center convention and the other Illumina?

    @ILunknown_unknown_1_1_107_394
    TATTCCCCGCCCCCCCGTCGTGCCCGGTCTTTGTCC
    +
    IIIIIIGI/<I8IIIIII.II//I);II)%<>III&
    @ILunknown_unknown_1_1_118_376
    TTGGGAAGCGCACCCGGCCCGTGTTGGCTTTCGCCT
    +
    I%I!I(,@+I"%",!#&)"5$'&$&%$.*%$"#"""
    @ILunknown_unknown_1_1_342_174
    GCGGCGGGTTATCGGGGCACTCCCCCCCTCCCGCAC
    +
    IIIIII2I<#?II+@CI1/CI0III5%(5'I%0-%$
    @ILunknown_unknown_1_1_121_440
    TGGCGTCCAGGCCGGCCTGGCTGCCGTCCGGCCCCC
    +
    II50I5I.&-*+9-&>2$&&3#&#"#%%"&"#&#&&
    @ILunknown_unknown_1_1_109_530
    TTTGCGTTAGGTAAAATGCTAGAAGCAGGTGAGCAG
    +
    IIIIIIIIIIIIIIIIIIIIIII&II0IIII%II0A

    thanks,
    -mark

  • #2
    The top is illumina/solexa format (what you would see in s_[1-8]_sequence.txt) and the bottom is the fastq Sanger format.

    To convert solexa values
    perl -e 'my $qual = 10 * log(1 + 10 ** ((ord($myChar) - 64) / 10.0))/log(10); print int($qual);'

    To convert fastq values
    perl -e 'my $qual = ord($myChar) -33; print $qual;'

    hope that helps

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Essential Discoveries and Tools in Epitranscriptomics
      by seqadmin


      The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist on Modified Bases...
      Yesterday, 07:01 AM
    • seqadmin
      Current Approaches to Protein Sequencing
      by seqadmin


      Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
      04-04-2024, 04:25 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 04-11-2024, 12:08 PM
    0 responses
    39 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 10:19 PM
    0 responses
    41 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 09:21 AM
    0 responses
    35 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-04-2024, 09:00 AM
    0 responses
    55 views
    0 likes
    Last Post seqadmin  
    Working...
    X