Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Split fastq to fasta and qual file?

    Hi all,

    Does anyone have or know about good scripts to split a sanger format fastq file into the corresponding fasta and qual file?? I have a dataset that I'd like to quality trim with LUCY but I can't figure out how to get it split apart! I've tried using the app on the galaxy page -- but its producing weird errors that I don't understand. Any help much appreciated!!

    -Lizzy

  • #2
    This can be done with Biopieces (www.biopieces.org):

    Code:
    read_fastq -i test.fq | write_454 -o test.fna -q test.fna.qual -x

    Cheers,


    Martin

    Comment


    • #3
      In Biopython the simplest way to do it is like this:

      Code:
      from Bio import SeqIO
      SeqIO.convert("example.fastq", "fastq", "example.fasta", "fasta")
      SeqIO.convert("example.fastq", "fastq", "example.qual", "qual")
      You can be more cunning if you want to avoid making two passes through the FASTQ, but the above should be pretty fast anyway.

      See also http://dx.doi.org/10.1093/nar/gkp1137 - I'd have suggested using EMBOSS seqret which can do FASTQ to FASTA, but I don't think it supports the QUAL format.

      Comment


      • #4
        thank you!! The Biopython script did the trick-- even for a python newbie!

        Comment


        • #5
          Do you know how to use Biopython to do the reverse? Fasta +qual = fastq?

          Comment


          • #6
            Well, Biopieces can do that as well:

            Code:
            read_454 -i test.fna -q test.qual | write_fastq -o test.fq -x

            In fact, Biopieces can also trim sequences based on quality scores by using trim_seq:


            Code:
            read_454 -i test.fna -q test.qual | trim_seq | write_fastq -o test.fq -x


            Martin
            Last edited by maasha; 01-06-2011, 10:36 AM.

            Comment


            • #7
              Originally posted by ewilbanks View Post
              Do you know how to use Biopython to do the reverse? Fasta +qual = fastq?
              Since you asked, yes, most easily done with the PairedFastaQualIterator function in the Bio.SeqIO.QualityIO module:

              Code:
              from Bio import SeqIO
              from Bio.SeqIO.QualityIO import PairedFastaQualIterator
              rec_iter = PairedFastaQualIterator(open("Quality/example.fasta"),
                                                 open("Quality/example.qual"))
              SeqIO.write(rec_iter, "Quality/temp.fastq", "fastq")
              This isn't quite as easy as the reverse since we need to take two input files and read over them in sync - and the high level functions in Bio.SeqIO are all intended for just one file. This example is based on the example in the documentation here:

              Comment


              • #8
                Thanks everyone!

                @maasha, I'll have to check it out! Does trim_seq accept sanger format qualities or only Solexa?

                Comment


                • #9
                  trim_seq works on Illumina type qualities.

                  read_fastq and read_454 convert to Illumina type qualities per default. Phred scores are automagically detected and converted. If you have Solexa scores there is a switch.

                  write_fastq output Illumina type qualities.

                  write_454 automagically convertes to decimal scores.



                  Cheers,


                  Martin

                  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
                  67 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X