Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Bowtie does not align my both ends as a pair

    Dear All

    I am facing a problem with bowtie. I generated a library using MetaSim software for two viruses (Serotypes). I used an empirical error model, and --><-- as read pairs layout, the number of reads is 5000.

    I am trying to paired-end align this data in Bowtie using the following command:

    bowtie -S ebwt_file -f -1 reads_1.fasta -2 reads_2.fasta file .sam

    The problem is that Bowtie is giving very poor results when the command above would executed:

    # reads with at least one reported alignment: 16 (0.32%)
    # reads that failed to align: 4984 (99.68%)

    While i am expecting to see the opposite figure. However; both ends are aligned perfectly when i ask Bowtie to align them individually.

    I need to align both ends as a pair using Bowtie, i can't work out what might be the reason behind these poor results. I tried to find out a pattern by examining the pair-algined reads ( 0.32% which is 16 reads), but i failed to find any interesting figures.

    Many thanks

  • #2
    What insert size did you simulate for your paired reads? bowtie's default for a maximum valid insert (-X/--maxins) is 250 bp. If the read pair does not meet that criteria bowtie will not call it a valid alignment.

    Comment


    • #3
      I am using 330 as mean insert size, i will try to scale it down and give it a try! I am wondering then, bowtie wont be the best choice to align long inserts. What would you suggest?

      thanks very much kmcarr.

      Comment


      • #4
        You don't have to scale your insert size down as bowtie can align reads from larger inserts, just enter appropriate values for your insert size with bowtie's minimum and maximum insert size parameters (-I and -X). There are a number of parameters which bowtie uses to determine what is and is not a valid alignment and these need to be tuned to your particular experiment.

        Comment


        • #5
          1) Try setting your min and max insert size options to a huge range
          2) Try bwa, which will align each on individually, and will return them, no matter what the insert size is.

          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
          69 views
          0 likes
          Last Post seqadmin  
          Working...
          X