Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • how to parse reads containing a particular sequence in any orientation

    I'm pretty new to bioinformatics and sorry if it's not worthy of asking here.

    I have a FLASH stitched fastq file from my paired end data, from which I want to sort the reads containing a particular sequence or part of that sequence in any orientation, and make a new fastq file with them. Is there any easy tool/code to do that?
    Last edited by PoorSeq; 10-30-2013, 10:25 PM.

  • #2
    I have a FLASH stitched fastq file from my paired end data, from which I want to sort the reads containing a particular sequence
    Okay, that's doable with grep, and very quick.

    or part of that sequence
    Wait, what? You want any subsequence? Will one base pair do? What are your limits?

    Comment


    • #3
      If you are trying to get reads that align to a particular sequence, try bowtie2. Not sure if that's your purpose though?

      Comment


      • #4
        Originally posted by PoorSeq View Post
        I'm pretty new to bioinformatics and sorry if it's not worthy of asking here.

        I have a FLASH stitched fastq file from my paired end data, from which I want to sort the reads containing a particular sequence or part of that sequence in any orientation, and make a new fastq file with them. Is there any easy tool/code to do that?
        Hi- As gringer pointed out, you should clarify what you are trying to do. Anyway... This sequence of unix commands reads a fastq file and outputs the reads in fastq format matching a regular expression. See if it helps:

        Code:
        ## Get reads containing substring AAA or its revcomp TTT
        gunzip -c fastq.fq.gz \
        | paste - - - - \
        | grep -P '^@.*?\t(.*?AAA.*?)|(.*?TTT.*?)\t\+' \
        | tr '\t' '\n' \
        | gzip > sub.fq.gz
        
        ## Example input fastq:
        @seq1
        ACTGAAACTG
        +comment
        IIIIIIIIII
        @seq2
        ACTGNNNCTGTTT
        +comment
        BBBBBBBBBBBBB
        @seq3
        CCCCCCCCCCCCC
        +comment
        BBBBBBBBBBTTT
        @seq4
        AAACCCCCCCCCC
        +comment
        BBBBBBBBBBTTT
        
        ## Output sub.fq.gz
        @seq1
        ACTGAAACTG
        +comment
        IIIIIIIIII
        @seq2
        ACTGNNNCTGTTT
        +comment
        BBBBBBBBBBBBB
        @seq4
        AAACCCCCCCCCC
        +comment
        BBBBBBBBBBTTT
        If your input is unzipped use "paste - - - - < fastq.fq" instead of "gunzip -c fastq.fq.gz \
        | paste - - - -"

        Comment


        • #5
          Thanks all for the answers, specially dariober for the code. I have also developed a bioawk code later which I was able to use. The code is:

          bioawk -c fastx '/SEQUENCE/ {print "@"$name; print $seq; print "+"; print $qual }' inut.fq > output.fq

          Yes, my aim is to collect all the reads that contains a sequence (or subsequence, at least 14nt) and make a new file.

          Comment


          • #6
            Originally posted by PoorSeq View Post
            bioawk -c fastx '/SEQUENCE/ {print "@"$name; print $seq; print "+"; print $qual }' inut.fq > output.fq
            Hmm, bioawk seems pretty neat.

            Note that what you've got there won't work for reverse complement orientation, so you'll need to have both forward and reverse included. Also, picking any 14+nt substring of SEQUENCE (or its reverse-complement) will be a bit trickier to implement.

            Comment


            • #7
              Good point gringer! however, because I have already flash-stitched the reads, I expect that the sequence I'm looking for will be in the same orientation in all reads. Still, I would see if I have any potholes. Also, I chose 14nt because my samples are from a bacteria with a genome size of 4.2 million, so I expect anything equal or more than 12nt will be unique.

              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, 02:20 PM
              0 responses
              9 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 06-07-2024, 06:58 AM
              0 responses
              181 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 06-06-2024, 08:18 AM
              0 responses
              228 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 06-06-2024, 08:04 AM
              0 responses
              184 views
              0 likes
              Last Post seqadmin  
              Working...
              X