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
              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
            • seqadmin
              The Impact of AI in Genomic Medicine
              by seqadmin



              Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
              02-26-2024, 02:07 PM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 03-14-2024, 06:13 AM
            0 responses
            34 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-08-2024, 08:03 AM
            0 responses
            72 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-07-2024, 08:13 AM
            0 responses
            82 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-06-2024, 09:51 AM
            0 responses
            68 views
            0 likes
            Last Post seqadmin  
            Working...
            X