Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Solid2FASTQ

    Hello

    I just got a question. I wanna align SOLID data with BWA and therefore transformed it in FASTQ files with the solid2fastaq.pl script. But the script does not take the whole sequence of the csfasta but excludes the first two bp. In my opinion it would make sense to exclude the first base, as it is not in color space but not two.

    Here an example:
    an entry of my csfasta files has following sequence
    >1_3_123
    T000111222333

    and the output in the FASTQ file is a sequence with:
    AACCCGGGTTT

    By comparing the sequences it can be seen that two bp are missing:
    T000111222333

    AACCCGGGTTT

    this is because of line 100 in the script delivered by BWA 0.5.7 in the method read1:
    $_ = substr(<$fhs>, 2);
    This funktions takes the sequence beginning from the 3rd position.

    Is there any reason for this or is this an error ?


    Thanks for answering
    martin

  • #2
    Originally posted by Martin R View Post
    Hello

    I just got a question. I wanna align SOLID data with BWA and therefore transformed it in FASTQ files with the solid2fastaq.pl script. But the script does not take the whole sequence of the csfasta but excludes the first two bp. In my opinion it would make sense to exclude the first base, as it is not in color space but not two.

    Here an example:
    an entry of my csfasta files has following sequence
    >1_3_123
    T000111222333

    and the output in the FASTQ file is a sequence with:
    AACCCGGGTTT

    By comparing the sequences it can be seen that two bp are missing:
    T000111222333

    AACCCGGGTTT

    this is because of line 100 in the script delivered by BWA 0.5.7 in the method read1:
    $_ = substr(<$fhs>, 2);
    This funktions takes the sequence beginning from the 3rd position.

    Is there any reason for this or is this an error ?


    Thanks for answering
    martin
    This is very subtle but important property of BWA: 2bp are lost from each read for SOLiD data.

    The first color encodes the adapter and the first base, and therefore is partly dependent on DNA not present in the underlying genome (the adapter). BWA trims the first color and adapter so you lose the first two bases (a 50bp reads is now 48bp). This is performed so when the color read is compared to the encoded reference (in color space), the colors from the read will match the color reference exactly (when there is no error). Including the first color (or adapter) would violate those rules. Nonetheless, you lose two bases from each read, which BWA should be able to include/recover during local alignment (Smith-Waterman/Needleman-Wunsch), but does not. Other software, like BFAST (I am the author), can use all the bases in the final alignment, so it is possible. You should bug the BWA author (lh3) to remove the 2bp base trimming as it is not going to be fixed by changing the solid2fastq.pl, but only through changing much of the BWA color space code.

    <rant>
    It also seems like SOLiD support is an afterthought of many alignment software authors. Such a shame, since it is a much more fun and challenging problem to try and properly align in color space.
    </rant>

    Comment


    • #3
      Originally posted by nilshomer View Post

      <rant>
      It also seems like SOLiD support is an afterthought of many alignment software authors. Such a shame, since it is a much more fun and challenging problem to try and properly align in color space.
      </rant>
      hahah agreed.. is <rant> a new BB code? ROFL
      http://kevin-gattaca.blogspot.com/

      Comment


      • #4
        Thank you very much for that fast answer.
        Now I understand the handling of neglecting the first 2bp.

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Strategies for Sequencing Challenging Samples
          by seqadmin


          Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
          03-22-2024, 06:39 AM
        • seqadmin
          Techniques and Challenges in Conservation Genomics
          by seqadmin



          The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

          Avian Conservation
          Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
          03-08-2024, 10:41 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, Yesterday, 06:37 PM
        0 responses
        10 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, Yesterday, 06:07 PM
        0 responses
        9 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-22-2024, 10:03 AM
        0 responses
        49 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-21-2024, 07:32 AM
        0 responses
        67 views
        0 likes
        Last Post seqadmin  
        Working...
        X