Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Aligning paired end Illumina data with Bowtie

    I have 2 fastq files read1.fastq and read2.fastq.
    If I align read1.fastq to the reference genome, I get a set of readids that align. Lets say a1, a2, a3, a4, a5, a6, a7, a8, a9 for simplicity.
    If I align read2.fastq to the reference genome, I get another set of readids: a1, a2, a4, a5, a6, a8, a9, a10.
    But when I align both together as paired end data using -1 and -2 Bowtie command line options, I get only a1, a9 (I expected a1, a2, a4, a5, a6, a8, a9, i.e. common to read1 and read2) along with a bunch of memory warnings that look like this:

    Warning: Exhausted best-first chunk memory for read HWI-XXXXX:103046AACXX:5:1101:4313:3969 1:N:0:CAGATC/1 (patid 7142); skipping read

    Can someone explain this?
    Thanks,
    Vishal

  • #2
    I don't know bowtie well, but what do the single end .sam entries look like for A2-8? Are they repetative? Are they way far apart?

    Comment


    • #3
      It seems that bowtie does not guarantee to find the best match, after several tries, it may give up.

      Comment


      • #4
        A couple of the default bowtie parameters for paired end data can cause you to lose hits you might want to keep.

        The 'exhausted memory' error can be fixed by adding --chunkmbs 512 to the bowtie options.

        Bowtie also has a strict limit (250bp) on the distance between paired sequences which is quite low. We tend to use --maxins 1000 to allow hits which are a bit further apart.

        Comment


        • #5
          I'm having similar problem. Currently using Bowtie for Illumina on Galaxy

          Single-end mapping with each dataset yield high number (similar number) of reads and coverage. However when I used paired-end mapping number of stacks were dramatically reduced (after applying filter). Regions of genome where you expect to have stacks aren't there?

          Tried increasing the insert size, it did yield more coverage.
          Is it because it needs to be in a proper pair?

          Also does seed length affect the outcome? E.g. sequences trimmed to 40bp and 25bp, for read1 and read2, respectively. When -l = 28, would read2 be considered as valid sequences?

          Comment


          • #6
            Yes, proper pairing does indeed matter. When both paired-ends are aligned together, the results will include only the alignments in which the paired-end are within a certain distance apart, and in the proper orientation. With SOLiD reads, properly paired-end reads must be directed toward each other, and on opposite strands. I'm not surprised that increasing the insertion size improves the results slightly. But when the reads of a pair are mapped to different chromosomes, then increasing the insertion size won't help.

            If your results are in BAM format, try running "flagstat" from samtools on it. This will give you pairing stats. In my experience, when the percentage of properly paired reads is very low (e.g. 5% - 10%), that indicates that the sequences may not be good.

            Adjusting the seed length probably won't improve the paired-end results. If you increase the seed length from 25 to 28, then read2 will probably have fewer alignments -- but they will be more confident alignments. If you shorten the seed length, then the reads will align in more places, and the number of paired-reads may increase slightly. But I believe that this approach simply creates false positives.

            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, 03-27-2024, 06:37 PM
            0 responses
            12 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-27-2024, 06:07 PM
            0 responses
            11 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-22-2024, 10:03 AM
            0 responses
            53 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