Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Has anyone used FastX Toolkit with IonTorrent data?

    Hi y'all,
    Has anyone used FastX Toolkit with IonTorrent data? I'm trying to analyze my fastq files and I get errors like:
    fastq_to_fasta: Invalid quality score value (char '0' ord 48 quality value -16) on line 4
    and something similar with fastx_quality_stats

    Am I right at thinking that the problem is that the quality values need to be expressed in the Solexa system and thus FastX won't work with Ion Torrent? That would be a shame...

    Thanks,
    Gustavo

  • #2
    I think it is -Q 33 or -Q 64 that specifies which PHRED quality scores it will use. By default it uses the Solexa quals.

    Comment


    • #3
      Dear mnkboy,
      Thanks for your prompt response.
      I do not see the -Q 33 or -Q 64 option either at fastq_to_fasta or at fastx_quality_stats:

      usage: fastx_quality_stats [-h] [-N] [-i INFILE] [-o OUTFILE]
      Part of FASTX Toolkit 0.0.13 by A. Gordon ([email protected])

      [-h] = This helpful help screen.
      [-i INFILE] = FASTQ input file. default is STDIN.
      [-o OUTFILE] = TEXT output file. default is STDOUT.
      [-N] = New output format (with more information per nucleotide/cycle)

      usage: fastq_to_fasta [-h] [-r] [-n] [-v] [-z] [-i INFILE] [-o OUTFILE]
      Part of FASTX Toolkit 0.0.13 by A. Gordon ([email protected])

      [-h] = This helpful help screen.
      [-r] = Rename sequence identifiers to numbers.
      [-n] = keep sequences with unknown (N) nucleotides.
      Default is to discard such sequences.
      [-v] = Verbose - report number of sequences.
      If [-o] is specified, report will be printed to STDOUT.
      If [-o] is not specified (and output goes to STDOUT),
      report will be printed to STDERR.
      [-z] = Compress output with GZIP.
      [-i INFILE] = FASTA/Q input file. default is STDIN.
      [-o OUTFILE] = FASTA output file. default is STDOUT.

      is there another way to tell the program I'm working on a -Q 64 environment?

      Thanks,
      Gustavo

      Comment


      • #4
        It is not in the documentation, which is why it is confusing. I think the default is -Q 64 and if you use the other you have to do -Q 33. I do not think it supports anything else.

        Comment


        • #5
          thanks, mnkyboy! It worked like a charm!

          Comment


          • #6
            I agree with mnkboy... it is not doccumented ... it even works with recent illumina pipeline...(CASSAVA 1.9 I guess)

            Comment


            • #7
              tools for Ion proton quality?

              Did you obtain an answer because I have the same problem with Ion Proton reads? or did you use an other tools that fastx?
              Thank you vb

              Originally posted by gkijak View Post
              Hi y'all,
              Has anyone used FastX Toolkit with IonTorrent data? I'm trying to analyze my fastq files and I get errors like:
              fastq_to_fasta: Invalid quality score value (char '0' ord 48 quality value -16) on line 4
              and something similar with fastx_quality_stats

              Am I right at thinking that the problem is that the quality values need to be expressed in the Solexa system and thus FastX won't work with Ion Torrent? That would be a shame...

              Thanks,
              Gustavo

              Comment


              • #8
                Originally posted by vbiaudet View Post
                Did you obtain an answer because I have the same problem with Ion Proton reads? or did you use an other tools that fastx?
                Thank you vb
                when you use "fastx_quality_stats", you should add "-Q 33", for its default is Illumina's coding format.

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