Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Unmatched read name error when converting csfasta to fastq with solid2fastq.pl

    I'm encountered with "unmatched read name" error when I converting paired SOLiD csfasta files to fastq by using solid2fastq.pl from bwa-0.5.9

    The error message:
    ================================
    ** unmatched read name: '>4_26_52_F3
    ' != '>4_26_52_F3
    '
    ================================

    I read there was suggestion to modify the script as mentioned by

    1. drio's post in PE SOLiD reads alignment by bwa thread.

    2. javijevi's post in bwa samse segmentation fault thread.

    I tired the modification and the "unmatched read name" error remained.

    Could anyone shares your solution for this?

    Thank you very much.

  • #2
    It's dieing in this routine:
    Code:
    [FONT="Courier New"]sub read1 {
      my $i = shift(@_);
      my $j = ($i-1)<<1;
      my ($key, $seq);
      my ($fhs, $fhq) = ($fhr[$j], $fhr[$j|1]);
      while (<$fhs>) {
            my $t = <$fhq>;
            if (/^>(\d+)_(\d+)_(\d+)_[FR]3/) {
              $key = sprintf("%.4d_%.4d_%.4d", $1, $2, $3); # this line could be improved on 64-bit machines
              die(qq/** unmatched read name: '$_' != '$_'\n/) unless ($_ eq $t);
              my $name = "$pre:$1_$2_$3/$i";
              $_ = substr(<$fhs>, 2);
              tr/0123./ACGTN/;
              my $s = $_;
              $_ = <$fhq>;
              s/-1\b/0/eg;
              s/^(\d+)\s*//;
              s/(\d+)\s*/chr($1+33)/eg;
              $seq = qq/\@$name\n$s+\n$_\n/;
              last;
            }
      }
      return defined($seq)? ($key, $seq) : ();
    }[/FONT]
    Things you might try are modifying this to print out $t and the line number. Please do it on a copy of solid2fastq.pl so you keep the original. You can look at that line number in both files to see what's up. I'm not really sure what's up with that extra carriage return(s) in the error message.
    Last edited by Richard Finney; 11-06-2011, 08:51 AM. Reason: format

    Comment


    • #3
      Thanks Richard Finney.

      I came across to this piece of modified solid2fastq.pl script by bgulko in bwa color-space index thread.

      I did minor changed as suggested by drio's post in PE SOLiD reads alignment by bwa thread for pairend reads on this modified solid2fasta.pl script.

      Without that minor changed on that script, it only converted the F3 csfasta to single end reads fastq.

      So far, the script worked with no complaint. I got read1.fastq and read2.fastq after the converting task done. However, I'm not yet carry out the mapping task. I will do it later.

      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
      7 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, Yesterday, 06:07 PM
      0 responses
      7 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
      66 views
      0 likes
      Last Post seqadmin  
      Working...
      X