Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • PFS
    Member
    • Mar 2010
    • 55

    #1

    How can I tell if an Illumina run is gone bad?

    Sorry for the very basic question. I am dealing with "raw" Illumina data for the first time and I am not sure how to pre-process the reads.

    The *sequence.txt file for each lane contains reads that have tons of 'N' at the end, with quality 'B'. (see below)

    How can I tell if this is normal or if something went wrong with the sequencing run?

    For example:
    @HWI-EAS413_0021:2:1:1159:4065#0/1
    TCGTGCCCGGGTAGCTCTGACTGGGCTGACTGTGGCTGAATACTTNAGNGACNANGAAGGTCANGAGATCGG
    +HWI-EAS413_0021:2:1:1159:4065#0/1
    dddb\dddcd^aaccdcdd\cdd^dadddc^c`ccdcd`d^`b``B\ZB`UYB]BUUU^][]_B_VWWWYT^
    @HWI-EAS413_0021:2:1:1159:11764#0/1
    GTGTGCCTGGTCATGCTGTGGTGGATCACCGTCCCAGGGCATTGGNGANTTNNGNATTTACGANATCGGAAG
    +HWI-EAS413_0021:2:1:1159:11764#0/1
    fffefffffffffeffffdff_ffefffffffeff`eedd`a^b`B\VB`]BB]B_aa\^_`_BaO[[Y]dc
  • dawe
    Senior Member
    • Apr 2009
    • 258

    #2
    GAP >= 1.5 (or OLB...) call bases in a slightly different way, so that if you have a bad quality read, instead of having low scores on each base, you'll see a series of at least 3 consecutive 'B' (or '#' if you converted in sanger) in the "bad zone" (but you'd better check the manual about the number of consecutives). If you see a number of reads mostly called with 'B' qual, you may have an issue. The reads you posted look good.
    This is an example of a bad run:

    @ILLUMINA-F3E58E_10:5:1:2798:3037
    ATGGGCCCTTTATTTCCCATCTTTTCCCAACTTTTT
    +
    A@=>@###############################
    @ILLUMINA-F3E58E_10:5:1:2798:15617
    CATGTATCATGGTCATTTAGCATTATCAAATGTCCT
    +
    EDEABBB@BA:@:??B@:555:@@::5:??@,<B:A
    @ILLUMINA-F3E58E_10:5:1:2798:10376
    TCCCAACAGAACAGCCCCAAGGGGCAAGGGGTTAAA
    +
    A7@#################################
    @ILLUMINA-F3E58E_10:5:1:2798:1024
    TNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
    +
    ####################################
    @ILLUMINA-F3E58E_10:5:1:2798:1640
    ATTGTTCCATCACCTCTGACTTTTTCCCCCACTTTT
    +
    @?##################################

    Comment

    • maubp
      Peter (Biopython etc)
      • Jul 2009
      • 1544

      #3
      This 'B' quality (PHRED 2) is the Illumina 1.5+ Read Segment Quality Control Indicator (RSQCI), see also:
      Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc

      Comment

      Latest Articles

      Collapse

      • SEQadmin2
        16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
        by SEQadmin2



        NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification, detection of antimicrobial...
        Today, 06:46 AM
      • SEQadmin2
        Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
        by SEQadmin2



        Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

        Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
        ...
        10-14-2026, 09:55 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, 10-28-2026, 09:57 AM
      0 responses
      22 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-28-2026, 09:51 AM
      0 responses
      14 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-22-2026, 11:13 AM
      0 responses
      67 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-22-2026, 10:48 AM
      0 responses
      47 views
      0 reactions
      Last Post SEQadmin2  
      Working...