Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • inesdesantiago
    Member
    • Jan 2009
    • 44

    #1

    What is the difference between sol2std and sol2sanger

    Hello, does any one know what is the difference between these two maq scripts (both available in the latest maq version 0.7.1):

    sol2sanger
    sol2sanger convert Solexa FASTQ to standard/Sanger FASTQ
    Usage: maq sol2sanger in.fastq out.fastq


    sol2std
    sol2std Convert Solexa/Illumina FastQ to the standard/Sanger FASTQ
    Usage: pearl fq_all2std.pl sol2std in.txt > out.fastq


    I tried both scripts with the file "_sequence.txt" that I got from the Illumina. Here are the results for the first sequence:

    s_sequence.txt
    @GAII01_3:8:1:136:1560
    GAGTTCAATTTTGTTAATTCCTTTCTGAAAATAGAT
    +GAII01_3:8:1:136:1560
    WWWWWWWWWWWWWWWWWWWWWWWWWWLWWWRRRIRO


    sol2sanger_s_sequence.fastq
    @GAII01_3:8:1:136:1560
    GAGTTCAATTTTGTTAATTCCTTTCTGAAAATAGAT
    +
    88888888888888888888888888-888333+30


    sol2std_s_sequence.fastq
    @GAII01_3:8:1:136:1560
    GAGTTCAATTTTGTTAATTCCTTTCTGAAAATAGAT
    +
    88888888888888888888888888-888333+30!


    The only difference between the fastq file generated by sol2sanger and sol2std scripts is that the sol2std writes a "!" in the end of the 4th line..
    I checked, it does this in every 4th line of the file..
    What does this means?
    What script is the best?

    Thanks!
    Ines de Santiago
  • regyre
    Member
    • Jan 2009
    • 11

    #2
    Hi,

    If you run the following command:
    Code:
    perl fq_all2std.pl instruction
    You get this output:
    Code:
    FASTQ format is first used in the Sanger Institute, and therefore
    we take the Sanger specification as the standard FASTQ. Although
    Solexa/Illumina reads file looks pretty much like the standard
    FASTQ, they are different in that the qualities are scaled
    differently. In the quality string, if you can see a character
    with its ASCII code higher than 90, probably your file is in the
    Solexa/Illumina format.
    
    Sometimes we also use an integer, instead of a single character,
    to explicitly show the qualities. In that case, negative
    qualities indicates that Solexa/Illumina qualities are used.
    You might find these two threads informative as well:

    Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc


    Any topic/question that does not fit into the subcategories below. If you're unsure of where to put something, ask in here!


    Well, to answer your question, both should provide the same output. However as far as I can see, the ! is not introduced on purpose, looks more like a bug to me.

    I've "annotated" the perl sol2std function of the fq_all2std.pl script (maq 0.7.1) below:

    Code:
    sub sol2std {
      my $max = 0;
      while (<>) {
    	if (/^@/) {
    	  
              # print the first entry line
              print;
    	 
              # print the seq line, skip the second ID line and 
              # get the content of the quality line
    
    	  # the actual problem is there, the carriage return 
              # char is captured as well	  
              $_ = <>; print; $_ = <>; $_ = <>;
              
              # my solution is to add this line:
              #  chomp;	  
              my @t = split('', $_);
    
              # here the at array contains all the quality 
              # chars + a carriage return 
     	  my $qual = '';
    
              # here the additional carriage return will be converted 
              # as well; therefore the ! which is introduced
    	  $qual .= $conv_table[ord($_)] for (@t);
    	  print "+\n$qual\n";
    	}
      }
    }

    I hope lh3 will read this thread and can comment on it.

    Therefore, I would recommand to use this for the time being
    Code:
    maq sol2sanger
    Best,

    Comment

    • inesdesantiago
      Member
      • Jan 2009
      • 44

      #3
      thanks very much!

      Comment

      • malachig
        Senior Member
        • Aug 2010
        • 117

        #4
        I have also encountered the trailing '!' produced by the perl sol2std function of the fq_all2std.pl script (maq 0.7.1). This causes an exception if the subsequent fastq file is fed into hmmSplicer. Thanks for the advice on using sol2sanger instead.

        Comment

        Latest Articles

        Collapse

        • 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
        • SEQadmin2
          New Genomics Technologies Take Aim at Long-Standing Limits
          by SEQadmin2


          Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

          We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
          ...
          09-28-2026, 10:25 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by SEQadmin2, 10-15-2026, 10:52 AM
        0 responses
        23 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-14-2026, 05:39 AM
        0 responses
        27 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-13-2026, 09:56 AM
        0 responses
        23 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-09-2026, 11:14 AM
        0 responses
        28 views
        0 reactions
        Last Post SEQadmin2  
        Working...