Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • Giles
    Member
    • Feb 2010
    • 39

    #1

    Fastq groomer from command line

    Does anyone know of a script that allows one to convert Qual scores from a fastq file in order to map w/ BWA or Bowtie using the command line? I'm familar with the Galaxy tools function, but Galaxy has really been unreliable lately, I'd love to move away from it whenever possible.
  • nilshomer
    Nils Homer
    • Nov 2008
    • 1283

    #2
    I may be confused by your question, but BWA/bowtie take fastq files as input.

    Comment

    • Giles
      Member
      • Feb 2010
      • 39

      #3
      The qual portion of fastq comes in different varieties. Bwa and bowtie only accept one format. The galaxy fastq groomer changes the qual score format to allow bwa or bowtie compatibility. I would like to run that script from the command line, and I wondered if anyone knew where to find it.

      Comment

      • kmcarr
        Senior Member
        • May 2008
        • 1181

        #4
        Originally posted by Giles View Post
        The qual portion of fastq comes in different varieties. Bwa and bowtie only accept one format...
        Not true. Both bwa and bowtie can use directly FASTQ files with Q-scores encoded as either Phred+33 (Sanger standard) or Phred+64 (Illumina prior to CASAVA 1.8). You just need to set the proper option on the command line.

        Code:
        bowtie: --phred-33 (default) or --phred-64
        
        bwa aln: (default is Phred+33, no option required) or -I (Phred+64).
        Check the manuals for both of these programs for specifics.

        Comment

        • maubp
          Peter (Biopython etc)
          • Jul 2009
          • 1544

          #5
          There are several other options for doing the conversion if you need to, e.g. EMBOSS seqret
          http://dx.doi.org/10.1093/nar/gkp1137

          Comment

          Latest Articles

          Collapse

          • SEQadmin2
            Beyond CRISPR/Cas9: Understand, Choose, and Use the Right Genome Editing Tool
            by SEQadmin2



            CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).

            Despite this, “CRISPR helped turn genome editing from a specialized technique into
            ...
            07-31-2026, 11:01 AM
          • SEQadmin2
            Proteomic Platforms: How to Choose the Right Analytical Strategy to Improve Detection and Clinical Applications
            by SEQadmin2


            Proteomics platforms are evolving rapidly, with advances in mass spectrometry and affinity-based approaches expanding what researchers can detect and at what scale. As the field moves toward deeper proteome coverage and clinical applications, scientists face an increasingly complex landscape of tools. This article will explore how researchers are navigating these choices to find the right platform for their work.

            The systematic characterization of the human proteome has
            ...
            07-20-2026, 11:48 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by SEQadmin2, 08-06-2026, 07:41 AM
          0 responses
          16 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 08-03-2026, 10:13 AM
          0 responses
          32 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 07-31-2026, 02:55 AM
          0 responses
          42 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 07-24-2026, 12:17 PM
          0 responses
          26 views
          0 reactions
          Last Post SEQadmin2  
          Working...