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
                  Best Practices for Single-Cell Sequencing Analysis
                  by seqadmin



                  While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
                  06-06-2024, 07:15 AM
                • seqadmin
                  Latest Developments in Precision Medicine
                  by seqadmin



                  Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

                  Somatic Genomics
                  “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
                  05-24-2024, 01:16 PM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, Yesterday, 06:54 AM
                0 responses
                10 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 06-14-2024, 07:24 AM
                0 responses
                17 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 06-13-2024, 08:58 AM
                0 responses
                16 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 06-12-2024, 02:20 PM
                0 responses
                17 views
                0 likes
                Last Post seqadmin  
                Working...
                X