Header Leaderboard Ad

Collapse

Bowtie: Faulty Combination of Paired Ends

Collapse

Announcement

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

  • Bowtie: Faulty Combination of Paired Ends

    Hi All,

    My problem is that Bowtie is aligning paired ends together that shouldn't be.

    The alignment uses the options, -m1 -v2 --best --strata -X800

    The beginning of the alignment output looks fine...

    Ex Output at top of file:

    Id0001 + chr1 10000 sequence quality
    Id0001 - chr1 10200 sequence quality
    ...

    The end of the file seems to align reads together that don't go together (the input is highly repetitive, so this will still align)

    Ex Output at bottom of file:

    ...
    Id0052 + chr1 10000 sequence quality
    Id0124 - chr1 10200 sequence quality

    So it doesn't look like it's just outputting singletons for these reads, but is actually trying to align them together.

    Importantly I've double checked that the two paired end inputs are matched. EVERY READ in the input is in the same location in the first paired end file as it's counter part in second paired end file!

    Lastly, I'm still using Bowtie 1 so this isn't something weird with Bowtie 2.

    Thanks ahead of time !!!!

  • #2
    I haven't solved this yet, but here's an update...

    When I break the input fastq into smaller subfiles and run those with the same options, I don't come across the problem (as often).

    Does anyone know if there's some sort of buffer or memory storage quirk with Bowtie that could be resulting in this? Again, I have very repetitive sequences, and I'm running them with tight restrictions, which is what I think might be causing the problem.

    Again, thanks ahead of time. I hope someone out there can help

    Comment


    • #3
      slink

      Can you post the full command as well as the bowtie version ?

      -A

      Comment


      • #4
        Sure thing!

        The version is 0.12.7 (64-bit)

        and the full command is:

        bowtie --best --strata -m1 -v2 -X800 [hg19 reference] -1[input.fastq] -2[input.fastq] >[output]

        Comment


        • #5
          Okay maybe I'm going overboard with this... but after I split the reads into subfiles and realigned with the same options as previously, I graphed the number of alignments / subfile on one axis and the existence of mixed reads on an alternate y-axis.

          As expected, the number of total alignments dropped when the file also contained these mixed reads.

          What was odd, though, is that there seemed to be apprehension? The number of alignments in a subfile would drop but not have mixed reads. Then, the next subfile would have the drastic drop in alignment as well as an occurrence of mixed reads. Almost as if the program was struggling in the previous subfile's alignment, but not quite dead yet. Then...all would be good and everything would go back to normal until the next occurrence. Again, this kind of activity sounds to me like a buffer system, but I'm not sure and would love any other ideas anyone has. For now I'm going to run the subfiles separately, I'm getting enough data to work with, but would love a resolution to the problem.

          Thanks again,

          Sara

          -- Attached is a picture of the graph
          Attached Files

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Targeted Sequencing: Choosing Between Hybridization Capture and Amplicon Sequencing
            by seqadmin




            Targeted sequencing is an effective way to sequence and analyze specific genomic regions of interest. This method enables researchers to focus their efforts on their desired targets, as opposed to other methods like whole genome sequencing that involve the sequencing of total DNA. Utilizing targeted sequencing is an attractive option for many researchers because it is often faster, more cost-effective, and only generates applicable data. While there are many approaches...
            03-10-2023, 05:31 AM
          • seqadmin
            Expert Advice on Automating Your Library Preparations
            by seqadmin



            Using automation to prepare sequencing libraries isn’t a new concept, and most researchers are aware that there are numerous benefits to automating this process. However, many labs are still hesitant to switch to automation and often believe that it’s not suitable for their lab. To combat these concerns, we’ll cover some of the key advantages, review the most important considerations, and get real-world advice from automation experts to remove any lingering anxieties....
            02-21-2023, 02:14 PM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 03-17-2023, 12:32 PM
          0 responses
          8 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-15-2023, 12:42 PM
          0 responses
          17 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-09-2023, 10:17 AM
          0 responses
          66 views
          1 like
          Last Post seqadmin  
          Started by seqadmin, 03-03-2023, 12:03 PM
          0 responses
          64 views
          0 likes
          Last Post seqadmin  
          Working...
          X