Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Illumina paired-end names

    I've had these reads in Fastq format from an Illumina machine but noticed the paired end names are not ending in the usual /1 and /2, but instead /1 and /3. Picard doesn't handle these properly. I want to check this doesn't have another meaning I should know about.

    Thanks.

    First end:
    Code:
    @SequenceIdentifier#0/1
    [base calls]
    +SequenceIdentifier#0/1
    [quality scores]
    second end:
    Code:
    @SequenceIdentifier#0/3
    [base calls]
    +SequenceIdentifier#0/3
    [quality scores]

  • #2
    These FASTQ files were probably not produced directly by the Illumina software, perhaps they were output by srf2fastq. srf2fastq records the read segment number after the slash. If the protocol used was a standard illumina index recipe then read 1 is segment 1, the index read is segment 2 and read 2 is segment 3. Alternatively the processing pipeline may be configured to discard the last base(s) of the reads; in this case the saved part of read 1 is segment 1, the discarded part of read 1 is segment 2 and then the saved part of read 2 is segment 3. As you can see in both of these situations the /3 will be added to what you are calling read 2.

    You can fix these names with sed:

    Code:
    # sed -i.bak '/^[@\+]<unique part of your sequence identifier>/ s/\/3$/\/2/' filename.fastq
    This will replace the 3 with a 2 at the end of your sequence identifier lines. It will save a copy of the original file as filename.fastq.bak

    Comment


    • #3
      Hey thanks for that. Thanks for the explanation.

      I wrote something myself earlier to do a similar thing, except I was also checking things were as I expected:
      Code:
      #!/usr/bin/perl -w
      
      use warnings;
      
      =pod
      
      This script fixes the second read pair names in FastQ files, from 
      /3 to /2. 
      
      Usage: 
      cat input-fastq.txt | perl correct-readpairs.pl > output.fq 
      
      =cut
      
      $row = 1;
      while(<STDIN>) { 
              if($row % 2) { # only every first and third row to act on.
                       if(! s/\/3$/\/2/o) {
                              die "We have not considered the possibility of the following line:\n$_";
                      }
              }
              print $_;
              $row++;
      }

      Comment


      • #4
        Originally posted by trickytank View Post
        I've had these reads in Fastq format from an Illumina machine but noticed the paired end names are not ending in the usual /1 and /2, but instead /1 and /3. Picard doesn't handle these properly. I want to check this doesn't have another meaning I should know about.

        Thanks.

        First end:
        Code:
        @SequenceIdentifier#0/1
        [base calls]
        +SequenceIdentifier#0/1
        [quality scores]
        second end:
        Code:
        @SequenceIdentifier#0/3
        [base calls]
        +SequenceIdentifier#0/3
        [quality scores]
        This was probably a run with an index read; in this case the index read gets the number "2", the pairs are "1" and "3".

        hth,
        Sven

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Quality Control Essentials for Next-Generation Sequencing Workflows
          by seqadmin




          Like all molecular biology applications, next-generation sequencing (NGS) workflows require diligent quality control (QC) measures to ensure accurate and reproducible results. Proper QC begins at nucleic acid extraction and continues all the way through to data analysis. This article outlines the key QC steps in an NGS workflow, along with the commonly used tools and techniques.

          Nucleic Acid Quality Control
          Preparing for NGS starts with isolating the...
          02-10-2025, 01:58 PM
        • seqadmin
          An Introduction to the Technologies Transforming Precision Medicine
          by seqadmin


          In recent years, precision medicine has become a major focus for researchers and healthcare professionals. This approach offers personalized treatment and wellness plans by utilizing insights from each person's unique biology and lifestyle to deliver more effective care. Its advancement relies on innovative technologies that enable a deeper understanding of individual variability. In a joint documentary with our colleagues at Biocompare, we examined the foundational principles of precision...
          01-27-2025, 07:46 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 02-07-2025, 09:30 AM
        0 responses
        72 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 02-05-2025, 10:34 AM
        0 responses
        113 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 02-03-2025, 09:07 AM
        0 responses
        90 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 01-31-2025, 08:31 AM
        0 responses
        49 views
        0 likes
        Last Post seqadmin  
        Working...
        X