Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Aligning only unique reads in Bowtie

    Hi everyone,

    I am aligning some ChIP-seq data using Bowtie. I have been using the -m option to throw out any reads with > 1 reportable alignment, but I would also like to try omitting non-unique reads. Is there a command line option to throw out any reads that are identical?

    Thanks!

  • #2
    If you just have the raw reads, you can use the "uniq" command in Linux to extract the unique reads (after sorting).

    SpliceMap: De novo detection of splice junctions from RNA-seq
    Download SpliceMap Comment here

    Comment


    • #3
      Thanks John!

      That sounds like it should be a useful command. I was just wondering if you could give me a little more detail.

      I have the ChIP-seq data as FASTQ files which I align using bowtie. Would I use the uniq command on the FASTQ prior file to alignment to generate another FASTQ containing only unique reads?

      i.e., prior to alignment, run uniq -u on the FASTQ?

      Thanks!

      Comment


      • #4
        No worries, but the method I suggested is a bit of a hack... It will require you to fiddle with the data a bit.

        Firstly, do you need to preserve the read-quality information? If so then it is probably best to write your own python or perl script to do it. I'm pretty sure there are existing tools to do this though... I just can't re-call off the top of my head.

        -----

        The method I suggested is to firstly extract the raw-reads from the FASTQ file by using
        instructions here



        Then sort the reads with http://en.wikipedia.org/wiki/Sort_(Unix)

        sort input_file > output_file

        Finally use "uniq"

        uniq -u input_file > output_file

        After you do this, you can align your reads using bowtie with the "-r" option for raw reads.
        SpliceMap: De novo detection of splice junctions from RNA-seq
        Download SpliceMap Comment here

        Comment


        • #5
          You can try fastx_collapser from http://hannonlab.cshl.edu/fastx_toolkit/

          Comment


          • #6
            Re: Aligning only unique reads in Bowtie

            I have few questions regarding the best practices that are adopted, in dealing with multiple alignments from a single read and presence of identical reads in the data (from Biology stand point) :

            I am curious, how important it is to deal with identical reads.
            Having many identical reads in data means something wrong with the
            experiment?

            What could be considered as max. cutoff value for the number of identical reads in the data, so as to not consider those reads?

            In the other case of a single read aligning at multiple places in a genome, what should be the cutoff value for number of multiple alignments, so as to not consider those reads?

            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
            10 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, Yesterday, 06:07 PM
            0 responses
            9 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-22-2024, 10:03 AM
            0 responses
            50 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