Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Bowtie2 cannot read gzip format files

    Hey, everyone!

    The help document of bowtie2 have said that the files of reads could be gzip'ed.

    When I use the paired reads and try the argument ( -1 <m1> -2 <m2>), the bowtie2 can read the gzip'ed files.

    However, if I use the gzip'ed unpaired reads and apply the argument (-U <r>), the bowtie2 suggest that "Error: reads file does not look like a FASTQ file". The gzip'ed file can not be recognized.

    Does anyone have come up with this problem?

    Thank You!

  • #2
    # HOW to add GUNZIP support to ANY program that doesn't have it:

    mkfifo file1.fifo
    mkfifo file2.fifo
    gunzip file1.gz > file.fifo &
    gunzip file2.gz > file.fifo &
    bowtie -1 file1.fifo -2 file2.fifo

    Comment


    • #3
      Originally posted by earonesty View Post
      # HOW to add GUNZIP support to ANY program that doesn't have it:

      mkfifo file1.fifo
      mkfifo file2.fifo
      gunzip file1.gz > file.fifo &
      gunzip file2.gz > file.fifo &
      bowtie -1 file1.fifo -2 file2.fifo
      Thank You Very Much!

      This is very useful!

      Comment


      • #4
        Originally posted by earonesty View Post
        # HOW to add GUNZIP support to ANY program that doesn't have it:

        mkfifo file1.fifo
        mkfifo file2.fifo
        gunzip file1.gz > file.fifo &
        gunzip file2.gz > file.fifo &
        bowtie -1 file1.fifo -2 file2.fifo
        In many cases this will work, but not all cases. This doesn't work for programs that attempt to seek through files (low-memory fastq randomisation being one example that I've recently seen).

        Comment


        • #5
          I needed to add a "-c" to that command, otherwise gunzip will delete the original compressed file, i.e.
          Code:
          mkfifo file1.fifo
          mkfifo file2.fifo
          gunzip -c file1.gz > file1.fifo &
          gunzip -c file2.gz > file2.fifo &

          Comment


          • #6
            Originally posted by dsenalik View Post
            I needed to add a "-c" to that command, otherwise gunzip will delete the original compressed file, i.e.
            Code:
            mkfifo file1.fifo
            mkfifo file2.fifo
            gunzip -c file1.gz > file1.fifo &
            gunzip -c file2.gz > file2.fifo &


            this patch fixes fifo support, and also adds built-in gzip support to bowtie 1.0 (which remains the fastest, fastq->sam program out there)

            An ultrafast memory-efficient short read aligner. Contribute to earonesty/bowtie development by creating an account on GitHub.

            Comment


            • #7
              FYI, zcat is often a shortcut to 'gunzip -c'

              Also under BASH you should be able do something like:

              bowtie -1 <( zcat file1.gz ) -2 <( zcat file2.gz )

              or in the original posters question

              bowtie2 -U <( zcat file.gz )

              Comment


              • #8
                actually there's a bug in bowtie ... a race condition that can make it unreliable. you should not use fifo's unless you used the patched version that supports them

                An ultrafast memory-efficient short read aligner. Contribute to earonesty/bowtie development by creating an account on GitHub.

                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
                11 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, Yesterday, 06:07 PM
                0 responses
                10 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-22-2024, 10:03 AM
                0 responses
                51 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-21-2024, 07:32 AM
                0 responses
                68 views
                0 likes
                Last Post seqadmin  
                Working...
                X